dsolve(diff(diff(y(x),x),x)-(1/30*diff(diff(diff(diff(p(x),x),x),x),x)+7/3*diff(diff(p(x),x),x)+a*p(x)+b)*y(x)=0,y(x))