Void main()

{ double S1= simpson(f,0.0,pi,20), // обчислення интегралу

S2= cub(pi)/4.0; // точне значення

printf("Calculated Integral = %f \n",S1);

printf("Exact value = %f \n",S2);

}

double erfc(double x);

Функція повертає значення додаткового інтеграла помилок:

,

де

.