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