XProfan

Título de la ventana "   Asym-Erf(), el Gauss-Fehlerfunktion en asymptotischer Näherung "
' Quelle: https://jean-pierre.moreau.pagesperso-orange.fr/Basic/asymerf_bas.txt
' (D) Demo 2016/11 transponiert después de XProfan 11.2a by P. Pájaro carpintero, Vienna / Austria
' Keine sin embargo geartete Gewähr! No warranty whatsoever!
'****************************************************
'*         Program to demonstrate ASYMERF           *
'* ------------------------------------------------ *
'* Reference: BASIC Scientific Subroutines, Vol. II *
'* by F.R. Ruckdeschel, BYTE/McGRAWW-HILL, 1981 [1].*
'* ------------------------------------------------ *
'* SAMPLE RUN:                                      *
'* Find the value of ERF(X)=2*Exp(-X*X)/SQRT(PI)    *
'*                                                  *
'* Entrada X ? 3                                      *
'* ERF(X)= .9999779 with error estimate= -.00000000 *
'* Number of terms evaluated qué 10                 *
'*                                                  *
'* Entrada X ? 4                                      *
'* ERF(X)= 1.0000000 with error estimate= 0.0000000 *
'* Number of terms evaluated qué 17                 *
'****************************************************
'  DEFINT I-N : DEFDBL A-H, O-Z
Ventana de Estilo 24:CLS:font 2:Declarar x!,y!,e!,n&:Conjunto("decimals",17)
Tests:
Imprimir "\n\n  El Gauss-Fehlerfunktion Erf(x) para valor x = ";:Entrada x!
'Gosub "S1000"
S1000
Imprimir
Imprimir "\n  Erf(x) para valor x= ";x!;"  beträgt: ";format$("0.######",y!)
PRINT "  con uno Fehleabschätzung a ";format$("##0.##########",e!)
PRINT
PRINT "          ";n&;" Terme wurde berücksichtigt.";
WaitInput
Goto "Tests"
'***********************************************************
'* Asymptotic series expansion of the integral of          *
'* 2 EXP(-X*X)/(X*SQRT(PI)), the normalized error function *
'* (ASYMERF). This program determines the values of the    *
'* above integrand using a asymptotic series which is     *
'* evaluated to the level of maximum accuracy.             *
'* The integral is from 0 to X. The input parameter, X     *
'* must be > 0. The results are returned en Y and Y1,      *
'* with the error measure en E. The number of terms used   *
'* is returned en N. The error is roughly equal to first   *
'* term neglected en the series summation.                 *
'* ------------------------------------------------------- *
'* Reference: A short table of integrals by B.O. Peirce,   *
'* Ginn and Company, 1957.                                 *
'***********************************************************

Proc S1000

    Declarar c1!,c2!,y1!
    n& = 1 : y! = 1 : c2! = 1 / (2 * Sqr(x!))
    'G1100:

    Repeat

        y! = y! - c2!
        n& = n& + 2: c1! = c2!
        c2! = -c1! * n& / (2 * x! * x!)
        'Test for divergence - The romper point is roughly N=X*X
        Case ABS(c2!) > Abs(c1!) : BREAK':= Goto "G1200"
        ' Continue summation

    Until 0':=  Goto "G1100"

    ' G1200:
    n& = (n& + 1) / 2
    e! = EXP(-x! * x!) / (x! * 1.772453850905516)
    y1! = y! * e!
    y! = 1 - y1!
    e! = e! * c2!
    RETORNO
    'End of file Asymerf.prf