procedure QTRAP(A, B:real; var S:real);
const
EPS = 0.000001; JMAX = 20;
var
OLDS:real; J:integer;
begin
OLDS:= -1E+30;
For J:= 1 To JMAX do
begin
TRAPZD(A, B,
unit Unit2;
interface
procedure TRIDAG(A, B, C, R:array of real;var U:array of real; N:integer);
implementation
procedure TRIDAG(A, B, C, R:array of real;var U:array of real; N:integer);
var
unit Unit2;
interface
procedure VANDER(X:array of real;var W:array of real;
Q:array of real; N:integer);
implementation
procedure VANDER(X:array of real;var W:a
Procedure FGAUSS(X:real; A:array of real;
var Y:real;var DYDA:array of real; NA:integer);
var
I,II:integer; ARG,EX,FAC:real;
begin
Y:=0;
For II:=1 To (NA DIV 3) do
