Mathc initiation/c34ca

Un livre de Wikilivres.
Aller à la navigation Aller à la recherche


Sommaire


Installer et compiler ces fichiers dans votre répertoire de travail.


Crystal Clear mimetype source c.png c17a.c
'
/* ---------------------------------- */
/* save as c17a.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "fa.h"
/* --------------------------------- */
int main(void)
{
double a,b,n,ny,M;
 
 clrscrn();
 printf(" f : (x,y)-> %s\n\n", feq);
 printf(" u :   (x)-> %s  \n", ueq);
 printf(" v :   (x)-> %s\n\n\n", veq);

 a = 0;     b = 1;
 n = 2*10; ny = 2*10;

 M = simpson_dydx(f, a,b,n, u,v,ny);

 printf(" With the simpson's rule.\n\n");
 printf("    (%.3f  (%s\n", b, veq);
 printf(" int(   int(       %s  dy dx = %.6f\n", feq, M);
 printf("    (%.3f  (%s\n\n\n", a, ueq);

 stop();

 return 0;
}
/* --------------------------------- */
Même résultat que la version avec l'algorithme des trapèzes à quelques décimales près.

Exemple de sortie écran :

 f : (x,y)-> x**2+y**2+1

 u :   (x)-> 0  
 v :   (x)-> 2-2*x


 With the simpson's rule.

    (1.000  (2-2*x
 int(   int(       x**2+y**2+1  dy dx = 1.833333
    (0.000  (0


 Press return to continue.