dsolve(diff(diff(y(x),x),x)+k^2*JacobiSN(x,k)*JacobiCN(x,k)/JacobiDN(x,k)*diff(y(x),x)+n^2*y(x)*JacobiDN(x,k)^2=0,y(x))