matrixA: array[1..3,1..3] of real =
                    ( (-19/20,1/5, 3/5),
                      (-1 ,0.1, 0.5),
                      (-0.01 ,0 ,1/200)
                    );
     One: array [1..3,1..3] of real = (
                      (1,0,0),
                      (0,1,0),
                      (0,0,1)
                    );
     U:array[1..3] of real = (1,1,1.1);

double matrixA[3][3] =
	{{-19/20,1/5, 3/5},
	 {-1 ,0.1, 0.5},
	 {-0.01 ,0 ,1/200}
        };
double One[3][3] =
        {{1,0,0},
         {0,1,0},
         {0,0,1}
        };
double U[3] ={1,1,1.1};

var
i,j,k,q:byte;
A,At,A1,A2,Ar,One1:array[1..3,1..3] of real;
delta,Det,S,alpha:real;
B,Z,U1:array[1..3] of real;
f:text;

BYTE i,j,k,q;
DOUBLE A[3][3],At[3][3],A1[3][3],A2[3][3],Ar[3][3],One1[3][3]
DOUBLE delta,Det,S,alpha
DOUBLE B[3],Z[3],U1[3]
FILE *f