step5;
h=0.0001;
fdx=(f(1+h,1)-f(1,1))/h;
quiver(1,f(1,1),1,fdx,'b','linewidth',2,'Maxheadsize',1);
text('string','g''(1)=f_x(1,1)','position',[2 3]);