Notes

#include <stdio.h>
#include <math.h>
#define N 50

double integral(double a, double b, int n, double (*f)(double)); //funzione che calcola l'integrale di una funzione
double function(double); //funzione da integrare
int main() {
double a, b; // estremi di integrazione
double I = 0; // valore dell'integrale
double num[N+1]; // numero di intervalli
double integrals[N+1]; // array dei valori di I calcolati per ogni num[i]
double reciproci[N+1]; // ascisse del grafico, ossia 1/num[i]
double A, B; // coefficienti della retta di equazione y=A+Bx
double sigma[N]; // errori nei vari integrali (che puoi sapere solo se sai il cazzo di valore dell'integrale)
double S;.
double Sx;
double Sy;
double Sxx;
double Sxy;
double Stt;
double t;
double tmp;
int i;
printf("inserisci gli estremi di integrazione: ");
scanf("%lf %lf", &a, &b);
for(i=1; i<=N; i++){
num[i] = (double) i;
reciproci[i]=1/num[i];
integrals[i]=integral(a, b, i, function);
printf("%lf %lfn", num[i], integrals[i]);
}
for(i=0; i<N; i++){
sigma[i] = integrals[i+1]-integrals[N+1];
S += 1./(sigma[i]*sigma[i]);
Sx += reciproci[i+1]/(sigma[i]*sigma[i]);
Sy += integrals[i+1]/(sigma[i]*sigma[i]);
Sxx += (reciproci[i+1]*reciproci[i+1])/(sigma[i]*sigma[i]);
Sxy += (reciproci[i+1]*integrals[i+1])/(sigma[i]*sigma[i]);
}
for(i=0; i<N; i++){
t = (1/sigma[i])*(reciproci[i+1]-(Sx/S));
Stt += t*t;
tmp = (t*integrals[i+1])/sigma[i];
}
B = tmp/Stt;
A = (Sy - Sx*B)/S;
printf("%lf %lfn", A, B);
}


double function(double x){
return sin(x);
}

double integral(double a, double b, int n, double (*f)(double)){
double In=0.;
double delta = (b-a)/(double) n;
double x=a;
int j;
for(j=0; j<n; j++){
In += f(x+0.5*delta);
x += delta;
}
In *= delta;
return In;
}