Il Linguaggio Fortran 90/95

6.6 Procedure intrinseche per gli array 355
1. Riceve il numero n delle coppie di dati (x,y) ed i relativi valori dall’unità chiamante.
2. Calcola le sommatorie di x, y, xy x 2 usando allo scopo le funzioni intrinseche SUM e
DOT_PRODUCT.
3. Calcola il valore dei coefficienti m e c.
4. Restituisce in uscita i valori di m e di c all’unità chiamante.
SUBROUTINE least_sq(x,y,n,m,c)
! Scopo: calcolare i coefficienti della retta interpolante
! ai minimi quadrati della dispersione di dati x-y
IMPLICIT NONE
! Variabili dummy
INTEGER, INTENT(IN) :: n
REAL, DIMENSION(n),INTENT(IN) :: x, y
REAL, INTENT(OUT) :: m, c
! Variabili locali
REAL :: sum_x, sum_xsq, sum_xy, sum_y
! Calcolo delle sommatorie
sum_x=SUM(x)
sum_y=SUM(y)
sum_xy=DOT_PRODUCT(x,y)
sum_xsq=DOT_PRODUCT(x,x)
! Calcolo dei coefficienti
m=(n*sum_xy-sum_x*sum_y)/(n*sum_xsq-sum_x*sum_x)
c=(sum_y-m*sum_x)/n
END SUBROUTINE least_sq
MATMUL(MATRIX_A,MATRIX_B)
Esegue il prodotto righe×colonne di due matrici numeriche o logiche.
Gli argomenti MATRIX_A che MATRIX_B devono essere array dello stesso tipo e di rango uno o
due (almeno uno dei due, però, deve avere rango due) ed inoltre è necessario che l’ampiezza della
prima (o dell’unica) dimensione di MATRIX_B sia uguale all’ampiezza dell’ultima (o dell’unica)
dimensione di MATRIX_A.
Il rango e la forma, invece, dipendono dal rango e dalla forma degli argomenti, alla maniera
seguente:
• Se MATRIX_A ha forma (/n,m/) e MATRIX_B ha forma (/m,k/), allora il risultato sarà un
array di rango due di forma (/n,k/).
• Se MATRIX_A ha forma (/m/) e MATRIX_B ha forma (/m,k/), allora il risultato sarà un
array di rango uno di forma (/k/).