Re: Looking for Mittag-Leffler function [message #37983] |
Wed, 11 February 2004 13:21  |
Craig Markwardt
Messages: 1869 Registered: November 1996
|
Senior Member |
|
|
mzkiss@wisc.edu (Miklos Kiss) writes:
> Hi everyone! I am trying to write a code which will generate the
> Mittag-Leffler function,
>
> E_alpha(x) = sum(k = 0 to infinity) (x^k)/gamma(alpha*k + 1).
>
> There are more general cases, but in this particular application,
> x is real, and alpha is between 0 and 1. Oh, and as a special case,
> if alpha = 1, then this reduces to exp(x). My problem is that coding
> it up as is works up to a point before reaching machine limits (x^k
> for large x and large k, as well as large values of the gamma
> function), but I need solutions for large x (x >= 20).
Since your question really isn't specific to IDL, perhaps you could
investigate how other authors might have done it. Possibilities to
check into include special function libraries like CEPHES or the Gnu
Scientific Library.
Good luck,
Craig
--
------------------------------------------------------------ --------------
Craig B. Markwardt, Ph.D. EMAIL: craigmnet@REMOVEcow.physics.wisc.edu
Astrophysics, IDL, Finance, Derivatives | Remove "net" for better response
------------------------------------------------------------ --------------
|
|
|