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