function fu=fun(th)
fu=2+.5*sin(th)-.1*cos(th)+cos(th).^2+cos(4*th);
