Funciones

//Funciones auxiliares de solucion del laboratorio 3

function[a1,a0]=linregre(x,y)
// Obtiene pendiente e intercepto de una regresion lineal
//
// Calling Sequence
//
//[m,c]=linregre(z,w)
//
// Parameters
//
//x:arreglo con valores de la variable independiente
//y: arreglo con valores de la variable dependiente
//
// Return
//
//a1:coeficiente del temrino lineal de la regresion.
//a0:termino independiente de la regresion
//
//Description
//
//Calcula la pendiente e intercepto de la recta de mejor
//ajuste cuadratico de un conjunto de pares (x,y).
//
n=length(x);
Sxy= (n.*(sum(x.*y)));
Sx=sum (x);
Sy=sum (y);
Sx2=sum(x.*x);
xmedia=(sum(x) )/n;
ymedia=(sum (y))/n ;

a1=(n.*Sxy - Sx.*Sy)/(n.*Sx2 - Sx.*Sx);
a0=ymedia-a1*xmedia;

endfunction