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