dsolve(diff(diff(diff(y(x),x),x),x)-(3*k^2*JacobiSN(z,x)^2+a)*diff(y(x),x)+(b+c*JacobiSN(z,x)^2-3*k^2*JacobiSN(z,x)*JacobiCN(z,x)*JacobiDN(z,x))*y(x)=0,y(x))