function fu=fun(th)
fu=(cos(th)+4*sin(th).^3);
