Mathc initiation/Fichiers h : c25a4

Sommaire ◀ Utilise la commande "Retour en Arrière" de ton navigateur.

Installer ce fichier dans votre répertoire de travail.

	x_fxy.h utilitaire

/* ---------------------------------- */
/* save as x_fxy.h                   */
/* --------------------------------- */
double fxy_x(
double (*P_f)(double x, double y),
double   h,
pt2d  p
)
{
double tplsh;
double tmnsh;

 tplsh = ((*P_f)(p.x+h,p.y));
 tmnsh = ((*P_f)(p.x-h,p.y));

 return(( tplsh-tmnsh)/(2.*h) );
}
/* --------------------------------- */
double fxy_y(
double (*P_f)(double x, double y),
double   h,
pt2d  p
)
{
double tplsh;
double tmnsh;

 tplsh = ((*P_f)(p.x,p.y+h));
 tmnsh = ((*P_f)(p.x,p.y-h));

 return(( tplsh-tmnsh)/(2.*h) );
}
/* --------------------------------- */
double fxy_xx(
double (*P_f)(double x, double y),
double   h,
pt2d  p
)
{
double t;
double tplsh;
double tmnsh;

 t     = ((*P_f)(p.x  , p.y));
 tplsh = ((*P_f)(p.x+h, p.y));
 tmnsh = ((*P_f)(p.x-h, p.y));

 return( (tplsh-2*t+tmnsh)/(h*h) );
}
/* --------------------------------- */
double fxy_yy(
double (*P_f)(double x, double y),
double   h,
pt2d  p
)
{
double t;
double tplsh;
double tmnsh;

 t     = ((*P_f)(p.x, p.y  ));
 tplsh = ((*P_f)(p.x, p.y+h));
 tmnsh = ((*P_f)(p.x, p.y-h));

 return( (tplsh-2*t+tmnsh)/(h*h) );
}
/* --------------------------------- */

Dans ce fichier il y a les fonctions pour calculer les dérivées partielles.

        First partial derivatives.

f_x  =  [f(x+h,y)-f(x-h,y)]/2h

f_y  =  [f(x,y+h)-f(x,y-h)]/2h

         Second partial derivatives.

f_xx =  [f(x+h,y)-2f(x,y)+f(x-h,y)]/h**2

f_yy =  [f(x,y+h)-2f(x,y)+f(x,y-h)]/h**2