plot (sn(x,0.3))^2 -(sin(sqrt(1-0.15666)*x)/sqrt(1-0.15666*cos^2(sqrt(1-0.15666)*x)))^2 from 0 to 5