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