Nikita Beloshenko Create Your Badge		Shop (under construction).

Beloshenko + 2007
------------------------------------------------------------------
Gaus.for
------------------------------------------------------------------
Program Gaus
c Gaus metodc go down 
integer i,j,n
real M(20,20),r
open(unit=1,file='Gaus.txt',status='Old')
read(1,*)n
read(1,*)((M(i,j),j=1,N+1),i=1,N)

6 r=M(1,1)

do 1 j=1,N+1 
M(1,j)=M(1,j)/r 
1 continue
do 2 i=2,N
r=M(i,1)
if (r.EQ.0) goto 9
do 3 j=1,N+1 
M(i,j)=M(i,j)/r
M(i,j)=M(i,j)-M(1,j)
3 continue
9 continue
2 continue
do 4 i=2,N
do 5 j=1,N+1
M(i-1,j)=M(i,j)
5 continue
4 continue
do 7 i=1,N
do 8 j=2,N+1
M(i,j-1)=M(i,j)
8 continue
7 continue
 N=N-1 
if (N.GE.2) goto 6

r=M(1,2)/M(1,1)
c r-last root of the system
close(1)
end

 


------------------------------------------------------------------ 
 Gaus.txt
 ------------------------------------------------------------------ 

3 

1 1 1 6 

1 -1 1 2
 
1 1 -1 0