unit2.pas

来自「1、说明： 本书中所有的常用数值算法子过程按书中的章数分别放在以C开头」· PAS 代码 · 共 32 行

unit Unit2;

interface
Function PROBKS(ALAM:real):real;

implementation
Function PROBKS(ALAM:real):real;
label 99;
const
    EPS1 = 0.001;    EPS2 = 0.00000001;
var 
    A2,FAC,TERMBF,TERM,AAA,PRO:real;  J:integer;
begin
    A2:=-2 * Sqr(ALAM);
    FAC:=2; 
    PRO:=0; 
    TERMBF:=0; 
    For J:=1 To 100 do
    begin
        TERM:=FAC * Exp(A2 * J*J);
        PRO:=PRO + TERM;
        AAA:=Abs(TERM);
        If (AAA < EPS1 * TERMBF) Or (AAA < EPS2 * PRO) Then
            goto 99;
        FAC:=-FAC;
        TERMBF:=Abs(TERM);
    end;
    PRO:=1;
99: PROBKS:=PRO;
end;
end.