pari | R Documentation |
Wrappers for the three elliptic functions of PARI
P.pari(z,Omega,pari.fun="ellwp",numerical=TRUE)
z |
Complex argument |
Omega |
Half periods |
pari.fun |
String giving the name of the function passed to
PARI. Values of |
numerical |
Boolean with default |
This function calls PARI via an R system()
call.
Returns an object with the same attributes as z
.
Function translates input into, for example,
“ellwp([1+1*I,2+3*I],1.111+5.1132*I)
” and pipes this string
directly into gp
.
The PARI system clearly has more powerful syntax than the basic version that I'm using here, but I can't (for example) figure out how to vectorize any of the calls.
Robin K. S. Hankin
## 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"))
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.