Mathc initiation/c33ca
Installer et compiler ces fichiers dans votre répertoire de travail.
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\n", feq);
printf(" u : (x)-> %s \n", ueq);
printf(" v : (x)-> %s\n\n\n", veq);
printf("\n\n\n");
a = 0;
b = 1;
n = 2*150;
ny = 2*150;
M = trapezoid_dydx(f, a,b,n, u,v,ny);
printf(" With the trapezoidal'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 de Sympson à 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 trapezoidal's rule.
(1.000 (2-2*x
int( int( x**2+y**2+1 dy dx = 1.833343
(0.000 (0
Press return to continue.