function y=fu1(x,a,b)
y=(x-1).^5-a+b*sin(x);