Il Linguaggio Fortran 90/95

138 Array
! *** Sezione esecutiva ***
! Valori iniziali di tentativo
T = 0
! Condizioni al contorno
T(1:10,1) = 1.0
T(1,1:10) = (/(0.1*j,j=10,1,-1)/)
! Inizio iterazioni
Told = T
niter = 0
DO
T(2:9,2:9) = (Told(1:8,2:9)+Told(3:10,2:9) &
+Told(2:9,1:8)+Told(2:9,3:10))/4.0
diff = MAXVAL(ABS(T(2:9,2:9)-Told(2:9,2:9)))
niter = niter+1
! Test di convergenza
PRINT *, "Iter n. ",niter, diff
IF (diff < 1.0E-4) THEN
EXIT
END IF
Told(2:9,2:9) = T(2:9,2:9)
END DO
! Stampa dei risultati
PRINT *, "Campo finale: "
DO i = 1,10
PRINT "(10F7.3)", Told(1:10,i)
END DO
END PROGRAM laplace
Si riporta, per completezza, uno stralcio dell’output:
Iter n. 1 0.4750000
Iter n. 2 0.1812500
Iter n. 3 0.1125000
...
Iter n. 94 1.116693E-04
Iter n. 95 1.049638E-04
Iter n. 96 9.861588E-05
Campo finale:
1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.900 0.890 0.879 0.866 0.848 0.820 0.774 0.686 0.494 0.000
0.800 0.782 0.762 0.737 0.705 0.660 0.590 0.477 0.291 0.000
0.700 0.675 0.648 0.617 0.577 0.524 0.449 0.342 0.192 0.000
0.600 0.570 0.539 0.504 0.462 0.409 0.340 0.249 0.134 0.000
0.500 0.467 0.434 0.399 0.359 0.312 0.254 0.181 0.095 0.000