function y=fu(x)
%y=cos(x);
%y=sin(2*x).*cos(x).^2;
y=(2+.5*sin(x)-.1*cos(x)+cos(x).^2+cos(4*x));
