dsolve(x*y(x)*diff(y(x),x)-y(x)^2+a*x^3*cos(x) = 0,y(x))