Algorithme de la Méthode des Rectangles pour le calcul d’intégral en Langage C

Algorithme de la Méthode des Rectangles pour le calcul d’intégral en Langage C :

On souhaite calculer l’intégrale de la fonction suivante :

f(x) = exp(-x)

Dans l’intervalle : [0, 1].

#include <stdio.h> #include <stdlib.h> #include<math.h> #define a 0.0 // Borne inférieure d'intégration #define b 1.0 // Borne supérieure d'intégration #define N 100 // Nombre de points double f(double x) { return exp(-x); // fonction à intégrer } void main() { double *x; x = (double *)calloc(N, sizeof(double)); // Allocation dynamique des abscisses xi double Somme, I; double h = (b-a)/N; // Pas d'intégration Somme = 0.0; for(int i=1;i<=N;i++) // l'indice s'arrête à N { x[i] = a +i*h; Somme += f(x[i] + h/2); // calcul de la somme } I = Somme*h; // valeur finale de l'intégrale printf("\n L'integrale I = %18.16g \n\n",I); // Affichage du résultat delete [] x; // Libération de la mémoire }