!GUSSE JORDEN METHOD

USE MSIMSL

REAL :: A(15,15) ,AINV(15,15), TEST(15,15) ,ACOPY(15,15) ,B(15,1),BCOPY(15,1)

A(1,:)=(/1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,1 ,1 ,1 ,1 ,4,5/)

A(2,:)=(/2 ,3 ,2 ,1 ,3 ,2 ,4 ,1 ,0 ,6 ,3 ,1 ,4 ,5 ,6 /)

A(3,:)=(/3 ,13,11,5 ,4 ,7 ,9 ,2 ,3 ,1 ,9 ,2 ,0 ,7 ,7 /)

A(4,:)=(/4 ,7 , 2,3 ,5 ,0 ,13, 2,1 ,22, 3,7 ,6 ,9 ,6 /)

A(5,:)=(/5 ,45,37,1 ,1 ,9 ,3 ,2 ,5 , 4,8 , 6,9 ,22,8 /)

A(6,:)=(/6 ,2 ,4 ,7 ,8 ,4,5 ,3 ,7 ,2 ,1 ,4 ,5 ,2 ,9 /)

A(7,:)=(/7 ,4 ,3 ,4 ,1 ,2 ,7 ,5 ,8 ,8 ,4 ,5 ,1 ,2 ,4 /)

A(8,:)=(/8 ,6 ,5,2 ,2 ,5 ,3 ,9 ,4 ,3 ,5 ,7 ,7 ,4 ,3 /)

A(9,:)=(/9 ,5 ,15,17,0 ,14,2 ,5 ,8 ,7 ,7 ,4,1 ,8 ,7 /)

A(10,:)=(/10,13,9 ,5 ,2 ,9 ,1 ,2 ,0 ,7 ,5,4 ,7 , 9,14/)

A(11,:)=(/11,7 ,12,2,15,2 ,12,7 ,16,2 ,6 ,3 ,2 ,5 ,3 /)

A(12,:)=(/12,6 ,6 ,7,2 ,2 ,5 ,3 ,9 ,4 ,3 ,5 ,7 ,3 ,9 /)

A(13,:)=(/13,5 ,5 ,1,17,0 ,14,2 ,5 ,8 ,7 ,7 ,14,1 ,8 /)

A(14,:)=(/14,2 ,1,13,1 ,9 ,8 ,7 ,1 ,2 ,2 ,4 ,7 ,7 ,8 /)

A(15,:)=(/15,8 ,3 ,0 ,7 ,9 ,1 ,9 ,4 ,2 ,1 ,2 ,1 ,2 ,6 /)

B(15,1 )=(/4 , 7,6,8,3,1 ,5,3 ,2 ,1 ,1,9,7,7,1 /)

n=15

AINV(1,:)= (/ 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0/)

AINV(2,:)= (/ 0,1,0,0,0,0,0,0,0,0,0,0,0,0,0/)

AINV(3,:)= (/ 0,0,1,0,0,0,0,0,0,0,0,0,0,0,0/)

AINV(4,:)= (/ 0,0,0,1,0,0,0,0,0,0,0,0,0,0,0/)

AINV(5,:)= (/ 0,0,0,0,1,0,0,0,0,0,0,0,0,0,0/)

AINV(6,:)= (/ 0,0,0,0,0,1,0,0,0,0,0,0,0,0,0/)

AINV(7,:)= (/ 0,0,0,0,0,0,1,0,0,0,0,0,0,0,0/)

AINV(8,:)= (/ 0,0,0,0,0,0,0,1,0,0,0,0,0,0,0/)

AINV(9,:)= (/ 0,0,0,0,0,0,0,0,1,0,0,0,0,0,0/)

AINV(10,:)=(/ 0,0,0,0,0,0,0,0,0,1,0,0,0,0,0/)

AINV(11,:)=(/ 0,0,0,0,0,0,0,0,0,0,1,0,0,0,0/)

AINV(12,:)=(/ 0,0,0,0,0,0,0,0,0,0,0,1,0,0,0/)

AINV(13,:)=(/ 0,0,0,0,0,0,0,0,0,0,0,0,1,0,0/)

AINV(14,:)=(/ 0,0,0,0,0,0,0,0,0,0,0,0,0,1,0/)

AINV(15,:)=(/ 0,0,0,0,0,0,0,0,0,0,0,0,0,0,1/)

BCOPY=B

ACOPY=A

n=15

DO K=1,n

zarib=1./A(k,k)

A(K,:)=zarib*A(K,:)

AINV(K,:)=zarib*AINV(K,:)

DO m=1,n

IF (m==k) CYCLE

zarib=-a(m,k)

A(m,:)=A(m,:)+zarib*A(k,:)

AINV(m,:)=AINV(m,:)+zarib*AINV(K,:)

end do

END DO

TEST= MATMUL(ACOPY,AINV)

print*,'A'

PRINT*,TRANSPOSE(A)

PRINT*,'AINV'

PRINT*,TRANSPOSE(AINV)

print*,'TEST'

PRINT*,TRANSPOSE (TEST)

END Matrix size (15 * 15) (to be solved by three methods of Gauss Gordon (g. J. M), method of removing the method of Gaussian eliminatiin method (g. E. M) and method of the upper triangle and the lower triangle (l u - factorization).. Hi please help me I really need it

n 15 Real 8 Gauss Jordan Method Gaussian Elimination LU Factorization X X2 X9 X + XS X6 X2 Xe X9 X 10 X 11 X 12 X 13 *14 n 15 Real 8 Gauss Jordan Method Gaussian Elimination LU Factorization X X2 X9 X + XS X6 X2 Xe X9 X 10 X 11 X 12 X 13 *14