function y = quad(x,p)
%
y = x .* x + p;