XProfan

Titre de la fenêtre "Ersatzformel pour Gauß´sches Fehlerintegral, sur 8 se mettre genau"
CLS:font 2:set("decimals",9)
declare MU!,sigma!,x!,x_schon_normiert&,a$,resultat![3]

Repeat

    Imprimer "\n Mittelwert qui Grundgesamtheit MU = ";:input MU!
    Imprimer "          Standardabweichung sigma = ";:input sigma!
    Imprimer "   Testwert (Merkmalsausprägung) x = ";:input x!
    a$=»:x_schon_normiert&=0:Imprimer " Ist x déjà normiert? ";
    input a$:a$=lower$(left$($ Trim(a$),1))
    cas (a$="j") or (a$="y") or (a$="1") or (a$="+"):x_schon_normiert&=1
    Gausskurvenintegral(MU!,sigma!,x!,x_schon_normiert&)
    'Ergebnisse eh bien dans globaler Variable Resultat![0..3]
    imprimer "\n  x-Wert      Dichte      Wahrscheinl.pour: Wert<x        Wert dans +\-x \n"+mkstr$("-",75)
    imprimer tab(2);resultat![0];" ";tab(22);resultat![1];" ";tab(42);resultat![2];" ";tab(62);resultat![3]
    imprimer mkstr$("-",75):imprimer

until 0

proc Gausskurvenintegral :parameters MU!,sigma!,x!,x_schon_normiert&

    Resultat![0]=x!'Rückgabearray extern deklariert
    casenot x_schon_normiert&: x!=si(sigma!<>0,(x!-MU!)/sigma!,(x!-MU!))
    declare F!,S2!,Einseitig!,Zweiseitig!
    F!=exp(-0.5*sqr(x!))/2.5066283:S2!=x!:x!=1/(1+0.2316419*abs(x!))
    Resultat![1]=F!
    Einseitig!=1-F!*x!*(0.31938153+x!*(-0.356563782+x!*(1.781477937+x!*(-1.821255978+x!*1.330274429))))
    Resultat![2]=Einseitig!
    cas s2!<0:Einseitig!=1-Einseitig!
    Zweiseitig!=abs(1-2*Einseitig!)
    Resultat![3]=Zweiseitig!

endproc