/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T E S T   D ' U N E   F O N C T I O N   E N   D E N T S   D E   S C I E                                                    */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/DentScie.01$c' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  "INCLUDES.01.I"
                                        /* Introduit le 20051116092839...                                                            */

extern    double    sin();

#define   PAS       0.01
                                        /* Pas de l'incrementation.                                                                  */

#define   NTERMES                                                                                                                       \
                    (2*3)                                                                                                               \
                                        /* Nombre d'iterations qui doit etre de preference pair.                                     */
main()
          {
          double    angle;
          double    minimum=+1e200,maximum=-1e200;

          for       (angle=0 ; angle<DEUX_PI ; angle=angle+PAS)
                    {
                    double    fonction=0;
                    int       n;

                    for       (n=0 ; n<(NTERMES/2) ; n++)
                              {
                              double    p1=(4*n)+1;
                              double    p3=(4*n)+3;

                              fonction=fonction+(sin(p1*angle)/(p1*p1));
                              fonction=fonction-(sin(p3*angle)/(p3*p3));
                              }

                    minimum=MIN2(minimum,fonction);
                    maximum=MAX2(maximum,fonction);
                    }

          printf("min=%.16f   max=%.16f\n",minimum,maximum);
          }