plotting Jacobian elliptic function

Question

How can I plot bellow function in MATHEMATICA, MATLAB, or MAPLE?

Here $cn$ is the Jacobi elliptic function:

$$ f(x)=0.7951556995+0.0048443005cn^2(0.4831042546x,0.1118800243)$$

Answer 1

In Mathematica:

Plot[0.7951556995 + 0.0048443005 JacobiCN[0.4831042546 x,0.1118800243]^2,{x,0,5}]

Answer 2

In Maple the function is JacobiCN, but the parameter is different: Maple uses the modulus $k$ where Mathematica uses the amplitude $k^2$. So to get caverac's plot in Maple, you would use

plot(0.7951556995 + 0.0048443005*JacobiCN(0.4831042546*x, sqrt(0.1118800243))^2, x=0..5);