diff(diff(y(x),x),x)*y(x)+a*(diff(y(x),x)^2+1)