function out=fce03(x);
out=atan(x^2-4*x+1);        


