elliptic (version 1.4-0)

pari: Wrappers for PARI functions

Description

Wrappers for the three elliptic functions of PARI

Usage

P.pari(z,Omega,pari.fun="ellwp",numerical=TRUE)

Arguments

z

Complex argument

Omega

Half periods

pari.fun

String giving the name of the function passed to PARI. Values of ellwp, ellsigma, and ellzeta, are acceptable here for the Weierstrass \(\wp\) function, the \(\sigma\) function, and the \(\zeta\) function respectively

numerical

Boolean with default TRUE meaning to return the complex value returned by PARI, and FALSE meaning to return the ascii string returned by PARI

Value

Returns an object with the same attributes as z.

Details

This function calls PARI via an R system() call.

References

http://www.parigp-home.de/

Examples

# NOT RUN {
 #this in a dontrun environment because it requires pari/gp 
z  <- seq(from=1,to=3+2i,len=34)
p <- c(1,1i)
plot(abs(P.pari(z=z,Omega=p) - P(z=z,Omega=p)))
plot(zeta(z=z,params=parameters(Omega=p))- P.pari(z=z,Omega=c(p),pari.fun="ellzeta"))

# }

Run the code above in your browser using DataCamp Workspace