function y=flagr3(x)
y = 1./(1+x.^2);
