dsolve(y(x)^m*x^n*(a*x*diff(y(x),x)+b*y(x))+alpha*x*diff(y(x),x)+beta*y(x) = 0,y(x))