#include <iostream.h>
#include <math.h>
void regression(float x[], float y[], int nb, float &m, float &b)
{
	float xy=0,sx=0,sy=0,x2=0;
	int i;
	for(i=0;i<nb;i++)
	{
		xy=xy+x[i]*y[i];
		sx=sx+x[i];
		sy=sy+y[i];
		x2=x2+x[i]*x[i];
	}
	b=(nb*xy-sx*sy)/(nb*x2-sx*sx);
	m=(sy-b*sx)/nb;
}
void main()
{
	float x[10]={2,4,6,8,10,12,14,16,18,20};
	float y[10]={-2.84247,-2.0781,0.416201,0.928943
				,2.568643,2.070369,4.520165,3.96398
				,5.672622,6.246995};
	float m,b;
	regression(x,y,10, m,b);
	cout<<"m="<<m<<"\tb="<<b<<endl;
}