function out=fun1(x)
out=x^2-sin(x);