Elementary Principles 400 - 9.35 (cont'd) C *DECREMENT BY...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
9- 48 C* *DECREMENT BY 50 DEG. AND LOOK FOR A SIGN IN PSI DO 10I = 1, 20 CALL PSICAL (T, PHI, PSI) IF ((PSIL*PSI).LT.0.0) GO TO 40 TLAST = T PSIL = PSI T = T – 50. 10 CONTINUE 40 IF (T.GE.0.0) GO TO 45 WRITE (3, 2) 2 FORMAT (1X, 'T LESS THAN ZERO -- ERROR') STOP *APPLY REGULA-FALSI 45 DO 50 I = 1, 20 IF (I.NE.1) T2L = T2 T2 = (T*PSIL-TLAST*PSI)/(PSIL-PSI) IF (ABS(T2-T2L).LT.0.01) GO TO 99 CALL PSICAL (T2, PHIT, PSIT) IF (PSIT.EQ.0) GO TO 99 IF ((PBIT*PBIL).GT.0.0) PSIL = PSIT IF ((PSIT*PSIL).GT.0.0) TLAST = T2 IF ((PSIT*PSI).GT.0.0) PSI = PSIT IF ((PSIT*PSI).GT.0.0) T = T2 50 CONTINUE IF (I.EQ.20) WRITE (3, 3) 3 FORMAT ('0', 'REGULA-FALSI DID NOT CONVERGE IN 20 ITERATIONS') 93 STOP END 1 * SUBROUTINE PSICAL (T, PHI, PSI)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the 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.

Ask a homework question - tutors are online