Il Linguaggio Fortran 90/95

36 Tipi ed espressioni
risultati imprevedibili. A titolo di esempio, l’espressione SQRT(2) è formalmente scorretta in
quanto la funzione radice quadrata, SQRT, non può essere usata con argomenti di tipo intero.
Le seguenti espressioni costituiscono esempi di impiego di funzioni intrinseche:
Espressione
SQRT(678.36)
Significato
√
678.36
EXP(y) ey a+SIN(0.5*teta) a + sin 1
2ϑ 1 a
0.5*(b+a/SINH(x))+LOG10(3) (b + ) + log 3
2 sinh x
b/(ABS(x-y)*eps) b/|x − y|ε
EXP(ABS(h/2.)) e |h/2|
SQRT(b*c+SQRT(a))
√
bc + a
Come si vede da questi esempi, una qualunque espressione può essere usata come argomento di
una funzione intrinseca: in particolare essa può ridursi ad una semplice costante, ad una variabile,
oppure può coinvolgere operazioni aritmetiche e/o riferimenti ad altre funzioni intrinseche.
E’, inoltre, permesso utilizzare come argomento di una funzione intrinseca la funzione stessa.
Se in un’espressione aritmetica compare un riferimento ad una funzione intrinseca, questa
viene valutata prima di effettuare qualunque altra operazione. Da ciò segue che le espressioni
che compaiono come argomenti di una funzione saranno le prime ad essere valutate e, se esse
utilizzano a loro volta una funzione, questa avrà la precedenza.
Il seguente programma mostra un semplice esempio di utilizzo di alcune delle più comuni
funzioni intrinseche matematiche: ATAN2 (per la valutazione della funzione arcotangente) eSQRT
(per il calcolo della radice quadrata). Il suo scopo è quello di convertire coordinate cartesiane
in coordinate sferiche:
PROGRAM globo
! Fornisce le coordinate sferiche di un punto note le sue coordinate
! cartesiane e verifica il risultato
IMPLICIT NONE
REAL :: r, theta, phi, x, y, z
! start
WRITE(*,*) " Inserisci x, y e z: "
READ(*,*) x,y,z
r = SQRT(x**2+y**2+z**2)
theta = ATAN2(y,x)
phi = ATAN2(SQRT(x**2+y**2),z)
WRITE(*,*) " Coordinate sferiche: "
WRITE(*,*) r, theta, phi
x = r*COS(theta)*SIN(phi)
y = r*SIN(theta)*SIN(phi)
z = r*COS(phi)
WRITE(*,*) " Coordinate cartesiane corrispondenti: "
WRITE(*,*) x,y,z