@lion137: @_13th_Dragon @enedil
Dziękuje za każdą pomoc, było jeszcze kilka błędów, które udało mi się rozwiązać i program działa zgodnie z oczekiwaniami:
#include <stdio.h>
#include <math.h>
double fun (double t, int c) {
double wynik;
double s = exp((gamma((c+1)/2.0)) - (gamma(c/2.0))) * (1 / sqrt(M_PI*c));
double o = (1+((t*t)/c));
double l = (-(c+1)/2);
wynik = s*(pow(o, l));
return wynik;
}
int main(){
double a, b, dx, ai, bi, suma_pole, i, n;
int c = 1089498;
/* Wczytanie granic przedziału oraz precyzji całkowania. */
a = 0;
b = 5633.107;
n = 5633107;
/* Obliczenie podprzedziału */
dx=(b-a)/n;
suma_pole = 0.0;
for(i=0.0; i<n; i++) {
ai=a+i*dx;
bi=ai+dx;
suma_pole+=dx*((fun(ai, c)+fun(bi, c))/2);
}
/* Wydruk wyniku na ekran */
printf("Wartość całki funkcji f(t) dla przedziału (a, b) wynosi %.6lf.\n", suma_pole);
return 0;
}