diff(diff(y(x),x),x)*y(x)+y(x)^2-a*x-b