Integration numérique par la méthode des trapèzes (c)

Soyez le premier à donner votre avis sur cette source.

Snippet vu 34 930 fois - Téléchargée 32 fois

Contenu du snippet

C’est une méthode de Newton-Cotes : On subdivise l’intervalle d’intégration [a ;b] en n sous intervalles de longueur égale h.
h=(b-a)/n
Sur chaque sous intervalle [xi;xi+1] on remplace f par son polynôme d’interpolation de degré 1 (sa corde) passant par les points (xi,f(xi)) et
(xi+1, f(xi+1)).

Algorithme : Intégrale f(x)dx entre [a,b] = h[ (f(a)+f(b))/2 + SOMME de i=1 à n-1 de f(a+i*h)]

Source / Exemple :


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

double f(double x)
{
return sqrt(x); 
}

int main(int argc, char *argv[])
{
  int i,n;
  double s,t,x;
  double a,b,h;

  a=1.0; /* Borne inferieure de l'intervalle */
  b=2.0; /* Borne superieure de l'intervalle */
  n=10;  /* Nombre d'itération */

  h=(b-a)/n;
  s=f(a)+f(b);

  for (i=0;i<=n-1;i++)
  {
   x=a+i*h;
   s=s+2*f(x);
   t=s*h/2;

  printf("\n %.16e",t);
  }

 return 0;

}

A voir également

Ajouter un commentaire Commentaires
merci
christelsahatsop
Messages postés
2
Date d'inscription
dimanche 23 juin 2013
Statut
Membre
Dernière intervention
25 juin 2013

24 juin 2013 à 07:03
je voudrai bien avoir la methode de dichotomie pour la resolution des integrales en c .merci
sanamanar
Messages postés
1
Date d'inscription
samedi 3 janvier 2009
Statut
Membre
Dernière intervention
17 janvier 2011

17 janv. 2011 à 20:19
merci
cs_donmilo1
Messages postés
1
Date d'inscription
mercredi 24 novembre 2004
Statut
Membre
Dernière intervention
24 novembre 2004

24 nov. 2004 à 14:29
je voudrai bien avoir la methode de samson 1/3 et samson 3/8 pour la resolution des integrales en c++
cs_kall333
Messages postés
7
Date d'inscription
samedi 15 mars 2003
Statut
Membre
Dernière intervention
16 juillet 2005

3 juil. 2004 à 08:23
Tout petit algo simple et limpide.
La simplicité est le reflet de l'intelligence.

Superbe!
Afficher les 6 commentaires

Vous n'êtes pas encore membre ?

inscrivez-vous, c'est gratuit et ça prend moins d'une minute !

Les membres obtiennent plus de réponses que les utilisateurs anonymes.

Le fait d'être membre vous permet d'avoir un suivi détaillé de vos demandes et codes sources.

Le fait d'être membre vous permet d'avoir des options supplémentaires.