1 2 3 4 0 0 -500 -1e3 -1.5e3 -2e3 jacobi_ns(k, u) k k=0 k=0.5 k=0.75 k=0.95 k=1 Jacobi Elliptic ns