7300 V»IEAS-VMEAS*OX«Z( I .d)»(0.5*( Y( I . J)+y( I . J+ 1 ) )-VDEL( I ) ) 



7300 GO TO 200 



740O 204 IF(J.EO. JMAX) GO TO 202 



7500 VMEAS«VMEAS+DX*0.5*(Y(I . J+ 1 ) -Y ( I . J- 1 ) ) •Z( I . J ) 



7 GOO GO TO 200 



7700 202 VMEAS=VMEAS*OX*Z( I , J)*(Y(I . J)-0.5*(Y( I ,d)+Y(I . J-1 )) ) 



7800 200 CONTINUE 



7900 C'PRIOR TO EOS, COMPUTE AND STORE SEVERAL VALUES WE NEED OVER AND OVER 



8000 C'BECAUSE COMPUTER CAN'T RAISE A NEG VALUE TO AN EXPONENT 



8100 C'MUST PRESERVE THE SIGN. 



8200 DO 400 II=1.IMAX 



8300 00 401 JJ='JBEGIN(II ).iJMAX 



8400 ARG1-V(II . Jd)-YDEL(Ili 



8500 DYSIGN=SIGN( 1 . ,ARG1) 



8600 DY=ABS(Y(II,dJ)-YDEL(II)) 



8700 DYTWOdI . JJ)-DY**EXPON 



8800 OVONEdl .dJ)=OYSIGN*OY»* THIRD 



8900 OVMTWOdi . JJ)=DY**(-EXPON) 



9000 OVMONE(II . dJ) *DYSIGN*DY**( -THIRD ) 



9100 DYMF0R(II,dd)=DY**(-2.*EXPGN) 



9200 DYFORdI ,dd)=DY**(2. *EXPON) 



9300 401 CONTINUE 



9400 ARG2'1400. -YDEL( II ) 



9500 DSIGN=>SIGN( 1 . . ARG2 ) 



9600 OYE'^ABSCARGZ) 



9700 YETWOdI j^DYE'^EXPON 



9800 YEONE( II )>=DSIGN*DYE»*THIRD 



9900 YEMONEdl ) =DSIGN*DYE" ( -THIRD ) 



lOOOO YEMTWO( II)=DYE*»(-EXPON) 



10100 YEMFORdI )=DYE**(-2. 'EXPON) 



10200 YEFIVEdl )=DSIGN'DYE»*(5. 'THIRD) 



10300 400 CONTINUE 



I0400 C'LET'S INPUT THE FIRST ROW OF THE MATRIX. A 



10500 SUM1B=0.0 



I0600 DO 300 II=1,IMAX 



10700 DO 300 vJd = dBEGIN( II ) .JMAX 



10800 300 SUM1B'SUMlB+2. »DYF0R(II ,dJ) 



10900 AMATRX( 1 , 1 ) = SUM1B 



11000 SUMLAM=0.0 



1 1 100 DO 305 K="1 , IMAX 



11200 SUM1K=0.0 



11300 DO 306 JJ=JBEGIN(K) .JMAX 



11400 306 SUMtK=SUM1K+2. •EXP0N*DYM0NE(K,dJ)»(Z(K,Ud)-2.*B* 



11500 • DYTWO(K,JJ)) 



1 1600 SUMLAM'SUMLAM-0.6»DX«YEFIVE(K) 



11700 305 AMATRX( 1 ,K+1 )-SUM1K 



11800 AMATRX( 1 , IMAX+2)»SUMLAM 



11900 C*NOW THE MIDDLE ROWS OF THE AMATRX. 



12000 DO 410 LR0W»2. IMAX+1 



12100 SUM2B-0.0 



12200 II^LROW-1 



12300 00 415 Jd=JBEGIN( II ) . JMAX 



12400 415 SUM2B=SUM2B+2. •EXPON«Z( II ,dd)*DYMQNE( II .dJ)-4 . *EXPON* 



12500 • B*DYONE( I I .dJ) 



12600 AMATRX(LROW. 1 ) =SUM2B + XLAMB*DX • YETWO( I I ) 



12700 DO 430 11 = 1. IMAX 



12800 SUM2Y»0.0 



12900 00 425 Od = JBEGINd I ) . JMAX 



13000 425 SUM2V«SUM2Y + 2. •EXP0N»THIRD*B*Z( 1 1 , dd ) ♦DYMFOR( I I . dd ) + THI RD*EXP0N 



13100 • •2. •B*B*DYMTWO( II .dd) 



13200 IF(dI*1).E0.LR0W) GO TO 431 



13300 *«ATHX(LROW. II+l )=0.0 



13400 GO TO 430 



13SOO 431 AMATRX(LROW, II + 1) = SUM2Y-XLAMB»EXP0N*DX*B*YEM0NEdI ) 



13600 430 CONTINUE 



13700 410 AMATRX(LROW. IMAX+2)=DX*B*YETW0(LR0W-1) 



13800 C*NOW THE LAST ROW OF THE MATRIX A 



13900 SUMFB^'O.O 



14000 DO 450 1 1' I. IMAX 



14100 450 SUMFB»SUMFB*0.6*DX*YEFIVEdI ) 



14200 AMATRX<IMAX>2. 1)=SUMFB 



109 



