elliptic integrals [message #74909] |
Wed, 09 February 2011 09:48  |
Jeremy Bailin
Messages: 618 Registered: April 2008
|
Senior Member |
|
|
Are there really no in-built functions to compute incomplete elliptic integrals? Or does someone have a pre-built routine already (save myself the effort of typing out and translating some NR code...)?
-Jeremy.
|
|
|
|
Re: elliptic integrals [message #75003 is a reply to message #74909] |
Wed, 09 February 2011 13:05  |
Foldy Lajos
Messages: 268 Registered: October 2001
|
Senior Member |
|
|
On Wed, 9 Feb 2011, Jeremy Bailin wrote:
> Are there really no in-built functions to compute incomplete elliptic integrals? Or does someone have a pre-built routine already (save myself the effort of typing out and translating some NR code...)?
>
> -Jeremy.
>
Have you tried GSL?
phi=...
k=...
val=call_external('/usr/lib64/libgsl.so', 'gsl_sf_ellint_f', $
double(phi), double(k), long(1), $
/all_value, /d_value, /auto_glue)
regards,
Lajos
|
|
|
|