%tegner en sinusfunksjon sin(x)
% a= er amplituden
close all
a=1;
x=0:pi/200:2*pi;                % velger omr?de og oppdeling i x
y = a*sin(x);                   % regner ut funksjonen
plot(x,y,'b','LineWidth',3)     % kurve farges blaa 'b' med linjetykkelse 3
legend('sin(x)')                % setter p? tekst
av = [0 2*pi -2 2];             % velger omr?det som skal plottes i x og y
axis equal
axis(av)
xlabel('x')                     % setter navn p? x-aksen                               
ylabel('y')                     % setter navn p? y-aksen
hold on                         % holder p? plottet til neste plott            
waitforbuttonpress              % venter p? tastetrykk (return) for ? g? videre
                                % NB! kursor m? st? i plottet
% f?rste leddet i Taylorrekken 
x1=0:pi/100:pi;
y1 = a*x1;
plot(x1,y1,'r','LineWidth',3)
hold on
axis(av)
legend('sin(x)','x')
waitforbuttonpress
% andre leddet i Taylorrekken
x2=0:pi/100:pi;
y2 = a*(x2-(x2.*x2.*x2)/6);
plot(x2,y2,'g','LineWidth',3)
hold on
axis(av)
waitforbuttonpress
% tredje leddet i Taylorrekken
x3=0:pi/100:pi;
%y3 = a*(x3-(x3.*x3.*x3)/6+(x3.*x3.*x3.*x3.*x3)/120);
y3 = a*(x3-(x3.^3)/6+(x3.^5)/120);
plot(x3,y3,'m','LineWidth',3)
hold off