-
DIMENSION PHINEW(NTOT), PHIOLD(NTOT)
-
NTOT1 = NTOT-1
1 |
BEGIN TEST=0
-
DO 1 I=2, NTOT1
-
TEMP = 0.5 * (PHIOLD(I-1) + PHIOLD(I+1))
-
TEST = AMAX1 (TEST, ABS (TEMP-PHIOLD(I)))
|
2 |
1 PHINEW (I) = TEMP
|
3 |
2 PHIOLD (I) = PHINEW(I)
-
IF (TEST > CONVG) GO TO BEGIN
-
ELSE STOP
|