(b + a*JacobiSN[x, k]^2)*y[x] + Derivative[2][y][x] == 0