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