dimanche 26 mai 2013

Algorithme de la Méthode du Point Fixe en Langage C

Algorithme de la Méthode du Point Fixe en Langage C :

Dans l’exemple ci-dessous Nous cherchons une solution à l’équation x = e-x avec une
précision de 10^(-15).

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 1000
#define EPS 1.0e-15 // Précision voulue
double g(double x)
{
return exp(-x); // fonction test
}
voidmain()
{
double x0,*x;
x = (double *)calloc(N,sizeof(double)); // allocation dynamique
printf("Entrez une valeur initiale X0:");
scanf("%lf",&x0);
x[0] = x0; // Initialisation de la valeur de départ
int k=0;
do
{
x[k+1] = g(x[k]); // calcul des termes de la série x(k+1)
if(fabs(x[k+1] - x[k]) <= EPS ) // critère de convergence
{
printf("\n\nNombre d'Iterations K = %d \t La solution X =
%18.16g\n", k, x[k+1]);
exit(0);
}
k++;
}
while(fabs(x[k+1] - x[k]) > EPS );// continuer les calcul tant que le critère n'est
pas atteind
delete [] x;
}

1 commentaire: