function y=p(x)
  y=cos(x);
%y=0;
