/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E   L A   F O N C T I O N   S I N U S  :                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrC/Fsinus.01$vv$I' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20130516094553).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

double    f3=1.0/(1.0*2.0*3.0);
double    f5=1.0/(1.0*2.0*3.0*4.0*5.0);
double    f7=1.0/(1.0*2.0*3.0*4.0*5.0*6.0*7.0);
double    f9=1.0/(1.0*2.0*3.0*4.0*5.0*6.0*7.0*8.0*9.0);

double    FonctionPrivee_Fsinus(double x)
          {
          double    x1=x,x2;
          double    signe=+1;

          x1=FonctionPrivee_Fmodulo(x1,6.283185307179586);

          if        (x1 > 3.14159265358979312)
                    {
                    signe=-1;
                    x1=SOUSf(x1,3.14159265358979312);
                    }
          else
                    {
                    }

          x2=MUL2f(x1,x1);

          return(MUL2f(signe,MUL2f(ADD2f(MUL2f(ADD2f(MUL2f(ADD2f(MUL2f(ADD2f(MUL2f(f9,x2),-f7),x2),f5),x2),-f3),x2),1),x1)));
          }



Copyright (c) Jean-François Colonna, 2013-2015.
Copyright (c) CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2013-2015.