#include <iostream.h>

float SSxy(float *x, float xbar,float *y,float ybar)
{
  float r = 0;
  for(int i = 0; i < 6; i++)
    r= r+(x[i]-xbar)*(y[i]-ybar);
  return r;
}

main()
{
  float w[6]={14.5, 17.9, 18.9, 19.1, 21, 23.2};
  float c[6]={20.9, 23.6, 27.3, 25.3, 28.1, 32};
  float wbar, cbar;
  float beta0,beta1;

  cout << "c  " << endl;
  for(int i = 0; i<6; i++)
    {wbar = wbar + w[i];
    cbar =cbar+c[i];
    };
  
  cbar = cbar/6;
  wbar = wbar/6;
  
  cout << "wbar: " << wbar << endl;
  cout << "cbar: " <<cbar << endl;

  cout << "d  " << endl;

  beta1=SSxy(w,wbar,c,cbar)/SSxy(w,wbar,w,wbar);
  beta0=cbar-beta1*wbar;

  cout << "beta0 " << beta0 << endl;
  cout << "beta1 " << beta1 << endl;

  float sse,s;

  sse = SSxy(c,cbar, c, cbar)-beta1*SSxy(w,wbar,c,cbar);
  s = sqrt(sse/4);
  cout<< "e"<<endl;
  cout << "sse "<< sse<< endl;
  cout << "s "<< s<< endl;
  cout << "s_beta1 "<<s/sqrt(SSxy(w,wbar,w,wbar))<<endl;


}