N=30;
disp(' ')
disp('Geplottet werden die Taylor-Polynome der Ordnung n = 0,...,30')
disp('an die Funktion  f(x) = sin(x) + cos(2x),  0 <= x <= 2*pi :')
disp(' ')
disp('Bitte Taste drücken ...')
disp(' ')
pause
x=linspace(0,2*pi,400);
f=sin(x)+cos(2*x);
clf;
h1=plot(x,f,'k');
set(h1,'LineWidth',2);
axis([0 7 -2 1.5]);
hold on;
pause
T=(2i).^[0:N];
T(2:2:N)=(-1).^[0:N/2-1];
T=T./[1 cumprod(1:N)];
co={'b' 'r' 'm' 'y' 'g' 'c'};
for n=1:N+1
  t=fliplr(T(1:n));
  p=polyval(t,x);
  plot(x,p,co{rem(n,6)+1});
  pause
end
disp('Testprogramm beendet')
disp(' ')