constante entier N <- 15

{ Type agrege de stockage de deux tableaux    }
{ de 15 double                                }
  
structure DeuxSeries
  x : Tableau [N] de reel
  y : Tableau [N] de reel
fin structure

{ Methode de calcul de la moyenne des valeurs }
{ contenues dans un tableau de double         }

réel fonction moyenne(t)
  Données
    t : Tableau [N] de réel
  Locales
    moyenne : reel
    i : entier
  moyenne <- 0.0
  pour i de 0 à N-1 faire
    moyenne <- moyenne+t[i]
  fait
  moyenne <- moyenne/N
  retourner moyenne
fin fonction

{ Methode de calcul de l'ecart quadratique    }
{ des valeurs contenues                       }
{ dans un tableau de double                   }

réel fonction ecartQuadratique(t,m)
  Données
    t : Tableau [N] de réel
    m : réel
  Locales
    v : réel
    i : entier
  v <- 0.0
  pour i de 0 à N-1 faire
    v <- v + (t[i]-m)*(t[i]-m)
  fait
  v <- sqrt(v)
  retourner v
fin fonction

{ Methode de calcul du coefficient            }
{ de correlation lineaire existant            }
{ entre deux series de valeurs reelles        }

réel fonction coefficientCorrelation(ds)
  Données
    ds : deuxSeries
  locales
    cc : réel
    i : entier
    mx : réel
    my : réel
    eqmx : réel
    eqmy : réel
  mx <- moyenne(ds.x)
  my <- moyenne(ds.y)
  cc <- 0.0
  pour i de 0 à N-1 faire
    cc <- cc + (ds.x[i]-mx)*(ds.y[i]-my)
  fait
  eqmx <- ecartQuadratique(ds.x,mx)
  eqmy <- ecartQuadratique(ds.y,my)
  cc <- cc/(eqmx*eqmy)
  retourner cc
fin fonction