function v=dfun(x)
v=cos(x);
%v=2*x;
return