Marsich and Merega 



B(l>iSAMP*CGt3,l)*(AA(1,L)*(CG(2, l)-F3^1*CG(%,l))+FC*AA(2,L)*CG(5 

 1,1)) 

 iG(L+6, l)-B(l) 



25 CG(L+8! I )-- CAMP*V*FC*CG(3, » )*CG(it. I )*AA(2,L) 



IIiDEF-5 

 GO TC 301 



26 ZZ(1.L)-DEF 

 DC 27 I^I.K 

 PI»6.2C*CG(l.l)/CLUf: 



AAdjVsiilFfPI) • ' 



AA(1,2)"-CCSF(PI) 



27 C( I )"n( I )*CG(1 , I )-CAMP*CG(3, I )*V*F3*CG (U, I )*AA(1 ,L) 



inDEr-6 



GO TC 301 



28 ZZ(2.L)=DEF 



29 cc;!t1::ue 

 c 



C ,"CTIC;; DIFFEr^.EilTIAL EOl'ATIC;; SCLUTICi! 



c 



DC 30 1=1,2 



Kl-=5-2*l 



Pli=2-I 



PJ«=2*I-3 



XA(1,I)=F1**2 *(CC(K1 U)-vCC(K1 2))-CC(!:i,l)+PI*V*(CC(2,3)+V*CC(1.2 



D) 



XA(2,I)=F1*CC(K1 3) 



X3(1.I)»F1**2 *CC(2,2)-CC(2,1)+PI*V*CC(1,3) 



30 X3(2, l)=F1*(CC(2,3)+PJ*V*CCn,2)) 

 DC 31 J"1,2 



KK«= 3-J 

 DC 31 1=1,2 

 K1=3-l 

 PI-= 2*1-3 



31 C3(JJ)-=XA(1,J)*ZZ(J,I)-PI*XA(2,J)*ZZ(J,K1)-XF.(1,J)*ZZ{K!' I)+PI*XB 

 1(2,jWz(KKh) 



PI=XA(2,1)*XA(2,2)-XA(1,1)*XA(1,2)+X3(1,1)*XB(1,2)-X3(2,1)*XB(2,2) 



PJ"XA(1,1)*XA(2,2)+XA(2,l)*XA(1,2)-XE(l,1)*XB(2,2)-X3(2,l)*Xn(1,2) 



DLT»PI*M+PJ*PJ 



DC 32 J-1,2 



DC 32 1-1,2 



K1-3-I 



D 1-3-2* I 



32 ZZ(J J )-(PI*BC(J, I )+01*PJ*BB(J,K1))/DLT 

 DC 5b \t^},k 



KrV(IM+1)/2 



Kr^(IM-1)/2 



KPh.|M-2*KN 

 56 VCM(IM, IC)-ZZ(KM,KP) 

 C 



C V/RITE (PvEAD) CN (FROM) AUXILIARY TAPE 2 (PARTIAL RESULTS) 



C 



WRITE TAPE 2,((CG(J3,J4),J>1,5).J4-1,K),((CG(J5,J6),J5-7,10),J6-1 



i.K).zzjc.Ft;(vcM(j7;ic);j>i;«f)' 



33 read" TAPE 2,((CG(J3,J'*),J>1.5),J^1,K),((CG(J5,J6),J5-7,10),J6-1 



READ TAPE 2, ((CG(J3,J^),J>1.5),. 

 |,K),ZZ,FC,Fl,(VGM(J7,IC),J7^1,if) 



866 



