/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T E S T   D E   L A   F O N C T I O N   ' fFCgamma_sommation_Euler_MacLaurin(...) '  :                                     */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtKg/gamma.21$K' :                                                                                             */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 19990621133356).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        I N T E R F A C E   ' listG '  :                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        :Debut_listG:                                                                                                              */
/*        :Fin_listG:                                                                                                                */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D I R E C T I V E S   S P E C I F I Q U E S   D E   C O M P I L A T I O N  :                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
@define   PRAGMA_CPP_____MODULE_NON_DERIVABLE

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        F I C H I E R S   D ' I N C L U D E S  :                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  INCLUDES_BASE

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M A C R O S   U T I L E S  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   PRINT(x,y)                                                                                                                    \
                    Bblock                                                                                                              \
                    DEFV(complexe,z);                                                                                                   \
                    DEFV(complexe,gamma);                                                                                               \
                    Cinitialisation(z,x,y);                                                                                             \
                                        /* Calcul de z=x+iy.                                                                         */ \
                    CALS(fFCgamma_sommation_Euler_MacLaurin(z));                                                                        \
                    Cegal(gamma,fFCgamma_sommation_Euler_MacLaurin_____valeur);                                                         \
                                        /* Calcul de l'approximation 'approximation' de 'gamma(z). On procede en deux temps afin     */ \
                                        /* de ne pas evaluer deux fois la fonction 'fFCgamma_sommation_Euler_MacLaurin(...)' a cause */ \
                                        /* de la programmation de 'Cegal(...)' ('v $ximcd/operator$FON FCegal').                     */ \
                    CAL2(Prin4("gamma(%+g,%+g)=(%+g,%+g)\n"                                                                             \
                              ,Reelle(z),Imaginaire(z)                                                                                  \
                              ,Reelle(gamma),Imaginaire(gamma)                                                                          \
                               )                                                                                                        \
                         );                                                                                                             \
                    Eblock

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T E S T   D E   L A   F O N C T I O N   ' fFCgamma_sommation_Euler_MacLaurin(...) '  :                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     /*..............................................................................................................................*/
     GET_ARGUMENTS_(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_F("precision_Bernoulli=""pb=",calcul_des_nombres_de_Bernoulli_____precision);
                         GET_ARGUMENT_F("pz=""precision_famma=",fFCfamma_sommation_Euler_MacLaurin_____precision_du_calcul);
                         GET_ARGUMENT_L("editer=",fFCfamma_sommation_Euler_MacLaurin_____editer_les_parametres_de_calcul);
                         GET_ARGUMENT_I("Vmin="
                                       ,fFCfamma_sommation_Euler_MacLaurin_____valeur_minimale_de_nombre_V_lors_du_calcul
                                        );
                         GET_ARGUMENT_I("Vmax="
                                       ,fFCfamma_sommation_Euler_MacLaurin_____valeur_maximale_de_nombre_V_lors_du_calcul
                                        );
                         )
                    );

     CAL2(Prin0("\n"));

     PRINT(+1.0,0.0);
     PRINT(+2.0,0.0);
     PRINT(+3.0,0.0);
     PRINT(+4.0,0.0);
     PRINT(+5.0,0.0);
     PRINT(+6.0,0.0);
     PRINT(+7.0,0.0);
     PRINT(+8.0,0.0);
     PRINT(+9.0,0.0);
                                        /* Valeur attendue gamma(n)=(n-1)!.                                                          */

     CAL2(Prin0("\n"));

     PRINT(+3.0,+2.0);
     PRINT(+3.0,-2.0);

                                        /* Valeur attendue : les deux valeurs obtenus doivent etre conjuguees.                       */
     CAL2(Prin0("\n"));

     RETU_Commande;
     Eblock
ECommande