dsolve(diff(diff(y(x),x),x)-(n*(n+1)*WeierstrassP(x,g2,g3)+B)*y(x)=0,y(x))