%ex integral clear all x=linspace(0,pi,10); y=sin(x); inty1=trapz(x,y) x=linspace(0,pi,100); y=sin(x); inty2=trapz(x,y) x=linspace(0,pi,1000); y=sin(x); inty3=trapz(x,y) inty4=integral(@sin,0,pi) %_________________________________ %fkn och dess derivata x=[-10:0.01:10]; y=x.^2; plot(x,y) title('y=x^2') dx=diff(x); dy=diff(y); dydx=dy./dx; figure(2) plot(x(1:length(x)-1),dydx) title('dy/dx')