![](/user_photo/_userpic.png)
книги / Статистический анализ данных в геологии. Кн. 2
.pdf1004 FORMAT)/’ Column 1 = X Variable' Г Column 2 = Y Variab'e 1‘ Column 3 = Y Value Based On Regression Equation'
2 |
Column 4 = Column 2 - Column 3') |
|
||
2000 FORMAT (У,28H Number of Observed Pairs = .15) |
||||
2001 |
FORMAT (25H To.al Sums of Squares = |
.F20.5) |
||
2002 |
FORMAT (37H Sums of Squares Due To Regression -=,F20 5 |
|||
2003 |
FORMAT (36H |
Sums of Squares Due To Deviation .= Г20 5) |
||
2004 |
FORMAT (19H |
Goodness of Fit = ,F15.6) |
.FI5.6) |
|
О |
5 FORMAT (27H Correlation Coefficient = |
|||
ILINO |
|
|
101 CONTINUE
B(l.l)-1.0
100 CONTINUE
DET-1.0
C
C CALCULATE INVERSE C
DO 102 1-1,N
CDIVIDE ITH ROW OF A AND В BY A(l,I) DIV=A(I,I)
IF(ABS(DIV) ZERO) 99,99,2
2DET=DET*DIV DO 103 J-1,N
A(I,J)=A(I,J)/D!V
B(I,J)=B(I,J)/QIV 103 CONTINUE
CREDUCE THE ITH COLUMN OF A TO ZERO DO 104 J=1 ,N
|
|
IF (l-J) 1,104,1 |
|
|
|
1 |
RATIO=A(J,l) |
|
|
|
DO 105 K-1.N |
|
|
|
|
|
A(J,K)=A(J,K) RATIO*A(l,K) |
|
|
|
|
B(J.K)=B(J,K) RATIO‘B(l,K) |
|
|
|
105 |
CONTINUE |
|
|
|
104 CONTINUE |
|
|
|
|
102 CONTINUE |
|
|
|
|
WRITE(I6,201) DET |
|
||
|
201 |
FORMAT(/' |
Determinant of Matrix A |
» F15 5) |
C |
RETURN |
|
|
|
|
|
|
|
|
C |
MATRIX CANNOT BE INVERTED A DIAGONAL ELEMENT IS APPROXIMATELY |
|||
C |
ZERO |
|
|
|
C |
|
|
|
|
|
99 DET - 0 0 |
|
|
|
|
WRITE(I6,202) |
|
||
|
202 FORMAT)/' |
Matrix Cannoi Be |
Inverted') |
|
|
LZ = 9 |
|
|
|
|
CALL PAUSE(2) |
|
||
|
RETURN |
|
|
|
C |
END |
|
|
|
|
|
|
|
|
CSLE |
SLE |
PROGRAM 4 9 |
|
C
C SUBROUTINE FOR SOLUTION OF N SIMULTANEOUS EQUATIONS
C MATRIX A IS N X N AND В IS A COLUMN VECTOR OF N ELEMENTS
CA IS CONVERTED TO THE IDENTITY MATRIX
СВ CONTAINS SOLUTION
C