This preview shows page 1. Sign up to view the full content.
Unformatted text preview: 100 XCM(I) = XC(I) 110 K = K + 1 CALL FUNCGEN (N, XA, XC) IF (IPR.EQ.1) CALL IPRNT (K, XA, XC, N) D0 200 I = 1, N IF (ABS ((XA(I) - XC(I))/XC(I)).GE.EPS) GOTO 300 200 CONTINUE C Convergence RETURN 300 IF(K.EQ.KMAX) GOTO 500 DO 400 I = 1, N W = (XC(I) XCM(I))/(XA(I) XAM(I)) Q = W/(W 1.) IF (Q.GT.0.5) Q = 0.5 IF (Q.LT.5) Q = 5. XCM(I) = XC(I) XAM(I) = XA(I) 400 XA(I) = Q = XAM(I) + (1. Q)*XCM(I) GOTO 110 500 WRITE (6, 900) 900 FORMAT (' CONVG did not converge') STOP END...
View Full Document
This note was uploaded on 02/27/2012 for the course CHEMICAL E 312 taught by Professor Cheung during the Fall '11 term at The University of Akron.
- Fall '11