diff(diff(diff(y(x),x),x),x)-(4*n*(n+1)*WeierstrassP(x,g2,g3)+a)*diff(y(x),x)-2*n*(n+1)*WeierstrassPPrime(x,g2,g3)*y(x)