Forum de mathématiques - Bibm@th.net
Vous n'êtes pas identifié(e).
- Contributions : Récentes | Sans réponse
#1 01-06-2010 18:10:01
- yoshi
- Modo Ferox
- Inscription : 20-11-2005
- Messages : 16 946
[C] Approximation de Pi (Méthode de Simpson)
Bonjour,
Approximation de [tex]\pi=\int_0^1 \frac{4}{1+x^2}\,dx[/tex] par la formule de Simpson.
Pour des explications se reporter au programme Python.
Je débute en C, donc le traitement des erreurs (et des accents) m'est inconnu : pour changer le pas, changer la valeur de p dans le code.
# include <stdlib.h>
# include <math.h>
int main()
{
double s=0;
double p = 0.001 ;
double pi =M_PI,nb=1/p;
int a = 0,b = 1,k=0;
printf(" Valeur Approchee de Pi\n\n");
printf("Par approximation de l'integrale sur [0 ; 1] de 4/(1+x^2)\n");
printf(" Formule de Simpson composee\n\n\n");
for(k=1;k<nb;k++)
{
s+=(1+k%2)*4/(1+pow(a+k*p,2));
}
s=p/3*(4/(1+pow(a,2))+4/(1+pow(b,2))+2*s);
printf("Valeur connue de Pi = %.16f\n",pi);
printf("Valeur obtenue de Pi = %.16f\n\n",s);
getchar();
return 0;
}
@+
Arx Tarpeia Capitoli proxima...
Hors ligne