function y=p(x)
ep=1.e-2;
y=x^2/ep;

