comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Looking for Mittag-Leffler function
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Looking for Mittag-Leffler function [message #37983] Wed, 11 February 2004 13:21 Go to previous message
Craig Markwardt is currently offline  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
------------------------------------------------------------ --------------
[Message index]
 
Read Message
Read Message
Previous Topic: Re: Savgol - Differentiation problems
Next Topic: Embedding IDL in Java

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sat Oct 11 02:32:37 PDT 2025

Total time taken to generate the page: 0.15391 seconds