#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include "lib.h"

/*********************** main ********************/
int main( int argc, char** argv ){
  float result;
  int val,i;
  float *lista_modelo, *lista_elem1, *lista_elem2;
  char dado[30],dado1[30];
  FILE *fp,*fp1;
  float XB=0.0, YB=0.0, S11=0.0, S22=0.0, B=0.0, A=0.0, *Y, S3=0.0, SE=0.0;
  char t[30];
  int aux;
    
   //--------------------------

   fp1 = fopen(argv[1],"r");
   if (fp1==NULL) {
	 printf("Nao consigo abrir o arquivo %s",argv[1]);
	 exit(-1);
   }

   fscanf(fp1,"%s\n",dado1);
   lista_elem1 = (float *) calloc(atoi(dado1),sizeof(float));

   XB=0.0;
   for(i=0; i < atoi(dado1); i++)
   {
	 fscanf(fp1,"%s ",t);
	 lista_elem1[i] =  atof(t);
       XB = XB + lista_elem1[i];
   }
   XB = XB / atoi(dado1);

   fp = fopen(argv[2],"r");
   if (fp==NULL) {
	 printf("Nao consigo abrir o arquivo %s",argv[2]);
	 exit(-1);
   }

   aux=0;
   while(!feof(fp))
   {
      fscanf(fp,"%s\n",dado);
      if(atoi(dado)!=0)
	{
      lista_elem2 = (float *) calloc(atoi(dado),sizeof(float));
      YB=0.0;
      for(i=0; i < atoi(dado); i++)
      {
	 fscanf(fp,"%s ",t);
	 lista_elem2[i] =  atof(t);
       YB = YB + lista_elem2[i];
      }
      YB = YB / atoi(dado);

      //--------------------------

      S22=0.0;
	for(i=0;i<atoi(dado);i++)
   	{
         S22 = S22 + (lista_elem1[i]-XB)*(lista_elem2[i]-YB);
   	}
	
      S11=0.0;
	for(i=0;i<atoi(dado);i++)
   	{
        	S11 = S11 + (lista_elem1[i]-XB)*(lista_elem1[i]-XB);
   	}

	B = S22 / S11;
	A = YB - B*XB;

	Y = (float *) calloc(atoi(dado),sizeof(float));
	for(i=0;i<atoi(dado);i++)
   	{
		Y[i] = A + B*lista_elem2[i];
	}
	S3=0.0;
	for(i=0;i<atoi(dado);i++)
   	{
		S3 = S3 + (lista_elem2[i]-Y[i])*(lista_elem2[i]-Y[i]);
	}	
	S3 = S3/(atoi(dado)-2);
	S3 = sqrt(S3);
	SE=0.0;
	for(i=0;i<atoi(dado);i++)
   	{
		SE = SE + (lista_elem1[i] - XB)*(lista_elem1[i] - XB);
	}
	SE = sqrt(SE);
	SE = S3 / SE;
        
      printf("regres=%.5f ",SE);

      //--------------------------
	}
   }

   printf("\n");
}