function y = fun1(x)
y = sin(10*x)+cos(3*x);