function y=quadf(x,s), y= -(x.^2)/(2*s); end
