g01eu returns the value of the Vavilov distribution function

Φ_{V} (λ; κ, β^{2})

It is intended to be used after a call to g01zu.

Syntax

C#
public static double g01eu( double x, double[] rcomm, out int ifail )

Visual Basic
Public Shared Function g01eu ( _ x As Double, _ rcomm As Double(), _ <OutAttribute> ByRef ifail As Integer _ ) As Double

Visual C++
public: static double g01eu( double x, array<double>^ rcomm, [OutAttribute] int% ifail )

F#
static member g01eu : x : float * rcomm : float[] * ifail : int byref -> float

Parameters

x: Type: System..::..Double
On entry: the argument $λ$ of the function.

rcomm: Type: array<System..::..Double>[]()[][]
An array of size [ $322$ ]
On entry: this must be the same parameter rcomm as returned by a previous call to g01zu.

ifail: Type: System..::..Int32%
On exit: $ifail = 0$ unless the method detects an error or a warning has been flagged (see [Error Indicators and Warnings]).

Return Value

g01eu returns the value of the Vavilov distribution function

Φ_{V} (λ; κ, β^{2})

It is intended to be used after a call to g01zu.

Description

g01eu evaluates an approximation to the Vavilov distribution function

Φ_{V} (λ; κ, β^{2})

given by

Φ_{V} (λ; κ, β^{2}) = \int_{- \infty}^{λ} ϕ_{V} (λ; κ, β^{2}) d λ,

where

ϕ (λ)

is described in g01mu. The method used is based on Fourier expansions. Further details can be found in Schorr (1974).

References

Schorr B (1974) Programs for the Landau and the Vavilov distributions and the corresponding random numbers Comp. Phys. Comm. 7 215–224

Error Indicators and Warnings

Errors or warnings detected by the method:

$ifail = 1$: Either the initialization method has not been called prior to the first call of this method or a communication array has become corrupted.

$ifail = -9000$: An error occured, see message report.
$ifail = -8000$: Negative dimension for array $〈value〉$
$ifail = -6000$: Invalid Parameters $〈value〉$

Accuracy

At least five significant digits are usually correct.

Parallelism and Performance

None.

Further Comments

g01eu can be called repeatedly with different values of

λ

provided that the values of

κ

and

β^{2}

remain unchanged between calls. Otherwise, g01zu must be called again. This is illustrated in [Example].

Example

This example evaluates

Φ_{V} (λ; κ, β^{2})

λ = 0.1

κ = 2.5

and

β^{2} = 0.7

, and prints the results.

Example program (C#): g01eue.cs

Example program data: g01eue.d

Example program results: g01eue.r