h=0.001;
x=0.8;

dfdx=(sin(x+h)-sin(x-h))/(2*h)