#include<stdio.h>
#include<math.h>

int main(void) {

  double x[100], y[100], ey[100],
    S1=0, Sx=0, Sx2=0, Sy=0, Sxy=0,
    D, a, sa, b, sb;
  int i,Ndata;

  printf("Δώσε το πλήθος των μετρήσεων: ");
  scanf("%d",&Ndata);

  if (Ndata>100) {
    printf("Το πλήθος των μετρήσεων δεν μπορεί να ξεπερνά τις 100.");
    return 1;
  }

  // εισαγωγή των δεδομένων
  for (int i=0; i<Ndata; i++) {
    printf("---> %dη μέτρηση:\n",i+1);
    printf("Δώσε το x: ");
    scanf("%lf",&x[i]);
    printf("Δώσε το y: ");
    scanf("%lf",&y[i]);
    printf("Δώσε το σφάλμα στο y: ");
    scanf("%lf",&ey[i]);
  }

  // υπολογισμός των αθροισμάτων S1, Sx, Sx2, Sy, Sxy
  for (i=0; i<Ndata; i++) {
    double sigmaFactor;
    sigmaFactor=1/ey[i]/ey[i];
    S1+=sigmaFactor;
    Sx+=x[i]*sigmaFactor;
    Sx2+=x[i]*x[i]*sigmaFactor;
    Sy+=y[i]*sigmaFactor;
    Sxy+=x[i]*y[i]*sigmaFactor;
  }

  // παράγοντας D
  D=S1*Sx2-Sx*Sx;

  // εκτύπωση παραγόντων
  printf("Πλήθος μετρήσεων = %d\n",Ndata);
  printf("   Sum(1/s2)     = %f\n",S1);
  printf("   Sum(x/s2)     = %f\n",Sx);
  printf("   Sum(x2/s2)    = %f\n",Sx2);
  printf("   Sum(y/s2)     = %f\n",Sy);
  printf("   Sum(xy/s2)    = %f\n",Sxy);
  printf("   D             = %f\n",D);

  // κλίση b
  b=(S1*Sxy-Sx*Sy)/D;
  sb=sqrt(S1/D);

  // τεταγμένη επί την αρχή a
  a=(Sx2*Sy-Sx*Sxy)/D;
  sa=sqrt(Sx2/D);

  // εκτύπωση αποτελεσμάτων
  printf("Ευθεία ελαχίστων τετραγώνων: y = a + b x\n");
  printf(" τεταγμένη επί την αρχή a = %f ± %f\n",a,sa);
  printf("                  κλίση b = %f ± %f\n",b,sb);

  return 0;
}