DSolve[(beta + alpha*JacobiSN[z, x]^2)*y[x] + b*Derivative[1][y][x] - (a + 12*k^2*JacobiSN[z, x]^2)*Derivative[2][y][x] + Derivative[4][y][x] == 0,y[x],x]