DSolve[((-(WeierstrassP[x, {a, b}]*(-a/2 + 6*WeierstrassP[x, {a, b}]^2)) - WeierstrassP[x, {a, b}]^2*WeierstrassPPrime[x, {a, b}] + WeierstrassPPrime[x, {a, b}]^2)*y[x])/(WeierstrassP[x, {a, b}]^2 + WeierstrassPPrime[x, {a, b}]) + ((a/2 - 6*WeierstrassP[x, {a, b}]^2 + WeierstrassP[x, {a, b}]^3 - WeierstrassP[x, {a, b}]*WeierstrassPPrime[x, {a, b}])*Derivative[1][y][x])/(-WeierstrassP[x, {a, b}]^2 + WeierstrassPPrime[x, {a, b}]) + Derivative[2][y][x] == 0,y[x],x]