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