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