NAG Library Routine Document

g01etf (prob_landau)

1
Purpose

g01etf returns the value of the Landau distribution function Φλ.

2
Specification

Fortran Interface
Function g01etf ( x)
Real (Kind=nag_wp):: g01etf
Real (Kind=nag_wp), Intent (In):: x
C Header Interface
#include <nagmk26.h>
double  g01etf_ (const double *x)

3
Description

g01etf evaluates an approximation to the Landau distribution function Φλ given by
Φλ=-λϕλdλ,  
where ϕλ is described in g01mtf, using piecewise approximation by rational functions. Further details can be found in Kölbig and Schorr (1984).

4
References

Kölbig K S and Schorr B (1984) A program package for the Landau distribution Comp. Phys. Comm. 31 97–111

5
Arguments

1:     x – Real (Kind=nag_wp)Input
On entry: the argument λ of the function.

6
Error Indicators and Warnings

None.

7
Accuracy

At least 7 significant digits are usually correct, but occasionally only 6. Such accuracy is normally considered to be adequate for applications in experimental physics.
Because of the asymptotic behaviour of Φλ, which is of the order of exp-exp-λ, underflow may occur on some machines when λ is moderately large and negative.

8
Parallelism and Performance

g01etf is not threaded in any implementation.

9
Further Comments

None.

10
Example

This example evaluates Φλ at λ=0.5, and prints the results.

10.1
Program Text

Program Text (g01etfe.f90)

10.2
Program Data

Program Data (g01etfe.d)

10.3
Program Results

Program Results (g01etfe.r)