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