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).
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; } |
Et si on ne veut pas enregistrer les termes dans un tableau dynamique !!
RépondreSupprimer