Il Linguaggio Fortran 90/95

102 Istruzioni di controllo
Tabellina del 5
5 per 1 = 5
5 per 2 = 10
5 per 3 = 15
5 per 4 = 20
5 per 5 = 25
5 per 6 = 30
5 per 7 = 35
5 per 8 = 40
5 per 9 = 45
5 per 10 = 50
Un altro esempio di applicazione del ciclo a conteggio è rappresentato dal seguente programma
che stima il valore di π con il Metodo Monte Carlo. Si rammenta brevemente il metodo in
esame. Generate tn coppie (x, y) di numeri random compresi tra 0 e 1 e detto in il numero di
punti che cadono nel quarto di cerchio inscritto nel quadrato di lato unitario, si approssima il
valore d π mediante il rapporto 4in/tn.
PROGRAM pi_greco
! Stima il valore di pi greco con il metodo Monte Carlo
IMPLICIT NONE
REAL :: x, y ! coordinate del colpo
INTEGER :: in=0, tn=0 ! colpi ’a segno’ e colpi totali
INTEGER :: loop ! variabile di controllo del ciclo
! start
WRITE(*,*) " Quanti ’colpi’? "
READ(*,*) tn
DO loop = 1,tn
CALL RANDOM_NUMBER(HARVEST=x)
CALL RANDOM_NUMBER(HARVEST=y)
IF(x**2 + y**2 < 1.0) in = in+1.0
END DO
WRITE(*,*) "Il valore stimato di pi_grego e’: ", 4.0*in/tn
STOP
END PROGRAM pi_greco
La procedura intrinseca RANDOM_NUMBER utilizzata nel programma serve a generare numeri
pseudo-random nell’intervallo 0÷1 e sarà descritta nel capitolo 5. Vale la pena controllare
la bontà del risultato prodotto:
C:\MYPROG>pi_greco
Quanti ’colpi’?
7000
Il valore stimato di pi_grego e’: 3.145714