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