XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX SOURCE PROGRAM XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



*16 6 



C 



C 



C 



C 



10 



1000 



3000 

 2000 



5000 



6000 

 4000 



500 

 100 

 200 



1620, FORTRAN II. PRODUCE MATRICES FOR DERIVING EQUATIONS OF PRMAT 01 



LINEAR SURFACES. PRMAT 02 



THEORY FROM W.C.KRUM8EIN . PROGRAMED BY BETTY BENSON, PRMAT 03 



MODIFIED BY W.S.WILSON, JULY 17. 1964. PRMAT 04 



DIMENSION X(4), Y(4). EM(4.3), S(3»3) PRMAT 05 



READ 10. (X{I),I=1»4) PRMAT 06 



READ 10, (Yd), 1 = 1, 4) PRMAT 07 



FORMAT (4F5.0) PRMAT 08 



DO 1000 L=1.4 PRMAT 09 



EMtL.l) = 1. PRMAT 10 



EM(L,2) = X(L) PRMAT 11 



EM(L»3) = Y(L) PRMAT 12 



DO 2000 1=1,3 PRMAT 13 



DO 2000 J=l,3 PRMAT 14 



S( I, J) = 0.0 PRMAT 15 



DO 3000 L=l,4 PRMAT 16 



S(I.J) = S(I.J) +EM(L»I)»EM(L.J) PRMAT 17 



CONTINUE PRMAT 18 



DO 4000 K=l,3 PRMAT 19 



DIV = S(K»IC) PRMAT 20 



S(K.K) = 1.0 PRMAT 21 



DO 5000 J=l,3 PRMAT 22 



S(K.J) = S(K.J)/DIV PRMAT 23 



DO 4000 1=1,3 PRMAT 24 



IF (I-O 1,4000,1 PRMAT 25 



DIV = S( I,K) PRMAT 26 



S(I,K) = 0.0 PRMAT 27 



DO 6000 J«l,3 PRMAT 28 



S(I,J) = S(I,J)-DIV*S(K,J) PRMAT 29 



CONTINUE PRMAT 30 



PUNCH 500 PRMAT 31 

 FORMAT (50HMATRICES FOR DERIVING EQUATIONS OF LINEAR SURFACES) PRMAT 32 



PUNCH 100, ((S(I,J),J=1»3) ,1=1,3) PRMAT 33 



FORMAT (6F12.8) PRMAT 34 



PUNCH 200, ( (EM(L,I),L=1,4) ,1=1,3) PRMAT 35 



FORMAT (12F6.2) PRMAT 36 



END PRMAT 37 



XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX INPUT XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



0.0 

 0.0 



1.0 

 0.0 



1.0 

 1.0 



0.0 

 1.0 



(X(I) ,I»1,4) 

 (Yd) ,1 = 1,4) 



XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX OUTPUT XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



MATRICES FOR DERIVING EQUATIONS OF LINEAR SURFACES 



.75000000 -.50000000 -.50000000 -.50000000 1.00000000 0.00000000 

 -.50000000 0.00000000 1.00000000 

 1.00 1,00 1.00 1.00 0.00 1,00 1.00 0.00 0.00 0.00 1.00 1.00 



39 



