dsolve(diff(diff(diff(y(x),x),x),x)+(-n^2+1)*WeierstrassP(x,g2,g3)*diff(y(x),x)+1/2*((-n^2+1)*WeierstrassPPrime(x,g2,g3)-a)*y(x)=0,y(x))