%DerivativePlotter

a=0;
b=1;
n=100;

x=linspace(a,b,n);
dx=x(2)-x(1);

f=x.^2;

upperdiag=ones(1,n-1);
lowerdiag=-1*ones(1,n-1);

D=diag(upperdiag,1)+diag(lowerdiag,-1);
D(1,1)=-3;
D(1,2)=4;
D(1,3)=-1;
D(n,n-2)=1;
D(n,n-1)=-4;
D(n,n)=3;

D=1/(2*dx)*D;

fp=D*f';

plot(x,fp);
