constante entier N <- 15

{ Type agrege de stockage de deux tableaux    }
{ de 15 reel                                  }
  
structure deuxSeries
  reel [N] x
  reel [N] y
fin structure

{ Fonction de calcul et retour de la moyenne  }
{ des valeurs contenues dans un tableau       }
{ de reel                                     }
{ t : Le tableau de calcul                    }

réel fonction moyenne(-> réel [] t)
  reel moyenne
  entier i
  moyenne <- 0.0
  pour i de 0 à longueur(t)-1 faire
    moyenne <- moyenne+t[i]
  fait
  moyenne <- moyenne/longueur(t)
  retourner moyenne
fin fonction

{ Fonction de calcul et retour de l'ecart     }
{ quadratique des valeurs contenues           }
{ dans un tableau de reel                     }
{ t : Le tableau de calcul                    }
{ m : La valeur par rapport à laquelle        }
{     le calcul est réalisé                   }

réel fonction ecartQuadratique(-> réel [] t,-> réel m)
  réel v
  entier i
  v <- 0.0
  pour i de 0 à longueur(t)-1 faire
    v <- v + (t[i]-m)*(t[i]-m)
  fait
  v <- sqrt(v)
  retourner v
fin fonction

{ Fonction de calcul et retour du coefficient }
{ de correlation lineaire existant            }
{ entre deux series de valeurs reelles        }
{ ds : Le paramètre de type deuxSeries        }
{      dans lequel sont stockées              }
{      les deux séries                        }    

réel fonction coefficientCorrelation(-> deuxSeries ds)
  réel cc
  entier i
  réel mx
  réel my
  réel eqmx
  réel eqmy
  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