dsolve(diff(diff(y(x),x),x)-(n*(n+1)*k^2*JacobiSN(x,k)^2+b)*y(x)=0,y(x))