XProfan

Título de la ventana "Mittels Newton-Raphson en el Genauigkeit verbesserte Lambert-W-Funktion"
Cls:Conjunto("decimals",17)
imprimir $profver:waitinput 7000

proc Wu :parámetros y!

    declarar Wu!

    if y!>1

        Wu!=ln(y!)-ln(ln(y!)) + 0.5*ln(ln(y!))/ln(y!)

    más

        Wu!=1

    endif

    volver Wu!

ENDPROC

proc Wo :parámetros x!

    declarar Wo!
    Wo!=ln(x!)-ln(ln(x!))+ exp(1)/(exp(1)-1)*ln(ln(x!))/ln(x!)
    volver Wo!

ENDPROC

Declarar xn!,x!,y!,epsilon!,i&,max&
y!=1.295' Stabil sólo al Haupt-Ast

Repeat

    if  y! > if($profver="11.2a-NT",10^305 , 10^53)

        sound 500,200:sound 300,200
        waitinput 20000:end

    endif

    y!=y!*1.1
    max&=4000'Verbesserungs-Läufe
    epsilon!=5/10^16:caso y!>227.93:epsilon!=5.687/10^14
    xn!=Wu(y!)'Startwert (a klein) para Newton-Raphson
    i&=0

    Repeat

        inc i&
        x!=xn!

        if ((x!+1)*exp(x!)) = 0

            xn!=(Wu(y!)+Wo(y!))/2'entonces eben de hier de Näherungsformel... de x~18.56
            romper

        más

            xn! = x!-(x!*exp(x!)-y!)/((x!+1)*exp(x!))

        endif

    until (abs(xn!-x!)<=epsilon!) or (i&>max&)

    if i&>max&

        imprimir " x-Búsqueda para y=";format$("%g",y!);" nach",i&,"Durchläufen abgebrochen!"
        imprimir tab(2);i&,
        imprimir tab(10);"Y=";format$("%g",y!),
        imprimir tab(23);"x=";format$("%g",xn!),
        imprimir tab(50);"err=";format$("%g",abs(x!-xn!))
        '   imprimir tab(50);"Y2=";format$("%g",xn!*exp(xn!))
        sound 2000,200:waitinput

    más

        imprimir i&,
        imprimir tab(5);"Y=";format$("%g",y!),
        imprimir tab(29);"x=";format$("%g",xn!),
        imprimir tab(52);"err=";format$("%g",abs(x!-xn!))
        '   imprimir tab(50);"Y2=";format$("%g",xn!*exp(xn!))

    endif

    if %csrlin>22:waitinput :locate 1,1:endif' cls:endif

    Until %key=27

    FIN