procedure GAUSSJ(VAR A:matrx2; N:integer; VAR B:array of real);
var
IPIV,INDXR,INDXC:array[1..50] of integer;
J,I,K,L,LL:integer;
BIG,PIVINV,DUM:real; IROW,ICOL:integer;
begin
Fo
procedure PENDAG(A,B,C,D,E,R:array of real;var U:array of real; N:integer);
var
W,BETA,ALPHA,CG,H:array[0..100] of real;
K:integer;
begin
W[1]:=C[1];
BETA[1]:=0;
BETA[2]:=D[
procedure SSOR(A:matrx2; N:integer; B:array of real;
var X:array of real; EPS,OM:real;var II:integer);
var
IMAX,J,I,L:integer; RX,R:real;
begin
IMAX:=200;
For I
procedure TRIDAG(A, B, C, R:array of real;var U:array of real; N:integer);
var
GAM:array[0..100] of real;
J:integer; BET:real;
begin
If B[1] = 0 Then Exit;
BET:=B[1];
U[1]:
procedure QRDCMP(var A:matrx2; M, N:integer;var Q:matrx2);
var
I,J,K:integer; S,T,Sgn,H,F:real;
begin
For I:=1 To M do
begin
For J:=1 To M do
Q[I, J]:=0;
procedure LUBKSB(var A:matrx2; N:integer; var INDX,B:array of real);
var
I,LL,J,II:integer;
SUM:real;
BEGIN
II:=0;
For I:=1 To N DO
begin
LL:=Trunc(INDX[I]);
procedure CHODCM(A:matrx2; N:integer;var D:array of real; T:array of real);
var
I,J,K:integer; SUM:real;
begin
For I:=1 To N do
begin
Sum:=A[I, I];
For J:=1 To I -
procedure CHOBSB(A:MATRX2; N:integer; D:array of real;var B:array of real);
var
I,J:integer; SUM:real;
begin
For I:=1 To N do
begin
Sum:=B[I];
For J:=1 To I - 1 do
unit Unit2;
interface
uses
unit1;
Function GASDEV:real;
implementation
Function GASDEV:real;
var
V1,V2,FAC,R:real;
begin
If ISET^= 0 Then
begin
repeat
V1:=2
