function y=flagr1(x)
y=sin(x)-x.*cos(x);
