diff(diff(y(x),x),x)+d+b*y(x)^2+c*y(x)+a*y(x)^3