step3;
h=0.0001;
fdx=(f(1+h,1)-f(1,1))/h;
quiver3(1,1,f(1,1),1,0,fdx,2,'color','b','linewidth',3,'maxheadsize',2);
