function [f df]=Newt_sec_1(x0,x1)
 fx = @(x)x^3-sin(x);
 f  = fx(x1);
 if nargout<2, return, end
 df = (fx(x1)-fx(x0))/(x1-x0);
end