numlib/src/int.pas

54                Var result, abserr, resabs, resasc: ArbFloat);
124    abserr := abs((resk-resg)*hlgth);
125    If (resasc <> 0) And (abserr <> 0)
126     Then abserr := resasc*amin1(1,exp(1.5*ln(200*abserr/resasc)));
128     Then abserr := amax1((50*macheps)*resabs, abserr);
236                Var result, abserr: ArbFloat;
250       abserr := giant;
291            abserr := err2+err3;
292            If abserr<5*macheps*abs(result)
293             Then abserr := 5*macheps*abs(result);
323               If (error<=abserr)
326                  abserr := error;
354          abserr := giant
358          abserr := abs(result-res3la[3]) +
364          If abserr<5*macheps*abs(result)
365           Then abserr := 5*macheps*abs(result)
370                 limit: ArbInt; Var result, abserr: ArbFloat;
395  abserr := 0;
411   qk21(f, a, b, result, abserr, defabs, resabs);
416  elist^[1] := abserr;
418   If (abserr <= 100*macheps*defabs) And (abserr>errbnd) Then ier := 2;
420   If (ier <> 0) Or ((abserr <= errbnd) And (abserr <> resabs)) Or (abserr=0)
432  errmax := abserr;
435   errsum := abserr;
436  abserr := giant;
549                               If (ktmin > 5) And (abserr < 1e-3*errsum)
551                               If abseps < abserr Then
554                                  abserr := abseps;
558                                   If abserr <= ertest Then go_on := false
593   If (abserr=giant) Or smallers Then p0 := true
598       If ierro=3 Then abserr := abserr+correc;
602       If abserr > errsum Then p0 := true
609       If abserr/abs(result) > errsum/abs(area) Then p0 := true
627   If Not p3 Then abserr := errsum;
649                 Var result, abserr, resabs, resasc: ArbFloat);
718       abserr := abs((resk-resg)*hlgth);
720       If (resasc<>0) And (abserr<>0)
723            reskh := 200*abserr/resasc;
725             Then abserr := resasc*reskh*sqrt(reskh)
726            Else abserr := resasc
733            If abserr<reskh Then abserr := reskh
743                 Var result, abserr: ArbFloat;
779       qk15i(f, bound, inf, 0, 1, result, abserr, defabs, resabs);
787       If (abserr<=100*macheps*defabs) And (abserr>errbnd) Then exit;
789       If ((abserr<=errbnd) And (abserr<>resabs)) Or (abserr=0) Then exit;
800       elist^[1] := abserr;
803       errmax    := abserr;
806       errsum    := abserr;
807       abserr    := giant;
927               If (ktmin>5) And (abserr<1e-3*errsum) Then ier := 4;
929               If (abseps<abserr)
933                   abserr := abseps;
938                   If (abserr<=ertest) Then continue := false
966       If (errsum<=errbnd) Or (abserr=giant) Then
969             abserr := errsum;
983       If ierro=3 Then abserr := abserr+correc;
987        If abserr/abs(result)>errsum/abs(area)
991            abserr := errsum;
1004       If abserr>errsum Then
1007             abserr := errsum;