Mathc initiation/Fichiers c : c32cd
Apparence
Installer et compiler ces fichiers dans votre répertoire de travail.
c00d.c |
---|
/* ---------------------------------- */
/* save as c00d.c */
/* ---------------------------------- */
#include "x_hfile.h"
#include "fd.h"
/* ---------------------------------- */
int main(void)
{
double ax = 1, ay = -1,
bx = 4, by = 2,
nx = 2*10, ny = 2*10,
Mxy = 0, Myx = 0;
clrscrn();
printf(" f : (x,y)-> %s\n\n\n", feq);
Myx = simpson_dydx(f, ay,by,ny, ax,bx,nx);
Mxy = simpson_dxdy(f, ax,bx,nx, ay,by,ny);
printf(" With the simpsonal's rule.\n\n");
printf(" (%.3f (%.3f\n", bx,by);
printf(" int( int( %s Dy Dx = %.6f\n", feq, Myx);
printf(" (%.3f (%.3f\n\n\n", ax,ay);
printf(" (%.3f (%.3f\n", by,bx);
printf(" int( int( %s Dx Dy = %.6f\n", feq, Mxy);
printf(" (%.3f (%.3f\n\n\n", ay,ax);
stop();
return 0;
}
/* ---------------------------------- */
/* ---------------------------------- */
Vous pouvez vérifier avec Octave version 5.2 :
I = dblquad (f, xa, xb, ya, yb)
>> I = dblquad (@(x, y) (2*x+6*x.*x.*y), 1, 4, -1, 2)
I = 234.00
Exemple de sortie écran :
f : (x,y)-> 2*x+6*x**2*y
With the simpsonal's rule.
(4.000 (2.000
int( int( 2*x+6*x**2*y Dy Dx = 234.000000
(1.000 (-1.000
(2.000 (4.000
int( int( 2*x+6*x**2*y Dx Dy = 234.000000
(-1.000 (1.000
Press return to continue.