function fu=fun(th)
fu=1+sin(th);
