x=linspace(-2,2,30);y=x; [X,Y]=meshgrid(x,y); f=@(x,y) (x.*y.^2)./(x.^2+y.^4) surfc(x,y,f(X,Y)); colorbar; shg figure contour(x,y,f(X,Y)) contour(x,y,f(X,Y),linspace(-1,1,20)) % tihennetään hilaa. x=linspace(-2,2,60);y=x; [X,Y]=meshgrid(x,y);contour(x,y,f(X,Y),linspace(-1,1,20));shg [c,h] = contour(x,y,f(X,Y)); clabel(c,h), colorbar
Tietokonearitmetiikkaa Katsotaan vielä polynomien evaluointiasiaa. Vrt. myös Lv7.mws Kirjoitimme Matlab-funktion horner.m (Otettin suoraan [BB]-kirjasta.) Tämä on matlab-hakemistossamme. Ellei, niin cut/paste ja UNIX:ssa: function p=horner(c,x) l=length(c);p=c(1); for k=2:l p=p.*x+c(k); end;Toimii myös vektorilla x.
# Luento viikolla 7, kiintopisteiteraatio # V2 kevät 2000 HA # Määritellään kiintopisteiteraatiota havainnollistavia # apufunktioita. # Määrittelemme ensin funktion porras, joka palauttaa # grafiikkaobjektin, joka puolestaan edustaa "porrasta" # [x,x],[x,f(x)],[f(x),f(x)].