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