// Régression linéaire avec Scilab
// Dominique Lefebvre - TangenteX.com
// Septembre 2010

// Définition des deux vecteurs de données
X = [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15];
Y = [0,9,19,29,41,50,62,69,83,95,101,109,119,130,142,149];

// calcul de la régression
coeff = regress(X,Y);

// extraction des coefficients
// Attention à l'inversion: Scilab pose Y = a + b*X alors
// que je pose Y = a*X + b, conventionnellement
a = coeff(2,:);
b = coeff(1,:);
printf("a: %f  b: %f \n",a,b);

// tracé du nuage de points et de la droite de régression
clf;
xtitle("Méthode des moindres carrés","Tension (V)","Intensité (mA)");
plot(X,Y,"+r");
Yr = a*X + b;
plot(X,Yr);