n=10:10:1000; k=length(n); s=zeros(1,k); err=zeros(1,k); f=@sin; a=0; b=pi; y=2; %=-cos(pi)-(-cos(0)) wahre Lösung for i=1:k s(i)=integriere(f,a,b,n(i)); err(i)=abs(y-s(i)); end plot(n,log(err)) title('Fehler gegen Anzahl der Stützstellen') xlabel('n') ylabel('|y-err|')