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

Soyez le premier à donner votre avis sur cette source.

Snippet vu 26 800 fois - Téléchargée 29 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

cs_orius
Messages postés
16
Date d'inscription
samedi 8 novembre 2003
Statut
Membre
Dernière intervention
3 juillet 2004
-
Jolie source compile sans problème sous Visual C++.

ça aurait été mieux je pense de mettre des long double dans les déclarations plutot que des doubles. Car si l'on change les bornes a et b le double et inssufisant je trouve :p
cs_kall333
Messages postés
7
Date d'inscription
samedi 15 mars 2003
Statut
Membre
Dernière intervention
16 juillet 2005
-
Tout petit algo simple et limpide.
La simplicité est le reflet de l'intelligence.

Superbe!
cs_donmilo1
Messages postés
1
Date d'inscription
mercredi 24 novembre 2004
Statut
Membre
Dernière intervention
24 novembre 2004
-
je voudrai bien avoir la methode de samson 1/3 et samson 3/8 pour la resolution des integrales en c++
sanamanar
Messages postés
1
Date d'inscription
samedi 3 janvier 2009
Statut
Membre
Dernière intervention
17 janvier 2011
-
merci
christelsahatsop
Messages postés
2
Date d'inscription
dimanche 23 juin 2013
Statut
Membre
Dernière intervention
25 juin 2013
-
je voudrai bien avoir la methode de dichotomie pour la resolution des integrales en c .merci

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.