dsolve(diff(diff(y(x),x),x)+(diff(y(x),x)-y(x)/x)^a*h(x,y(x))=0,y(x))