% Taylor-utvikling av (1-y^2)*cos(x) rundt origo
% over kvadratisk område i xy-planet
maks=5;
n=20;
[x,y]=meshgrid(linspace(-maks,maks,n));
hold off
% funksjonen i solid flate
surf(x,y,cos(x).*(1-y.^2));
hold on
% Taylor-utvikling til potens 2 i nett
mesh(x,y,1-0.5*x.^2 - y.^2);
xlabel('X')
ylabel('Y')
% bruk musa på figuren for å dreie den rundt slik at du får et godt
% inntrykk av hvor godt Taylor-utviklingen tilnærmer funksjonen i x- og
% y-retning
