/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        E T U D E   D E   L A   P E R T E   D E S   P R O P R I E T E S   D ' A S S O C I A T I V I T E                            */
/*        D E   L A   M U L T I P L I C A T I O N   E T   D E   L A   D I S T R I B U T I V I T E   D E                              */
/*        L A   M U L T I P L I C A T I O N   P A R   R A P P O R T   A   L ' A D D I T I O N   L O R S                              */
/*        D E   L ' E T U D E   D E   L A   D Y N A M I Q U E   D E   V E R H U L S T  :                                             */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*             *     * * * * * *   * * * * * *   * * * * * *   *         *   * * * * * *   *   * * * * * *   *         *             */
/*                        *             *        *             **        *        *        *   *         *   **        *             */
/*            * *         *             *        *             * *       *        *        *   *         *   * *       *             */
/*                        *             *        *             *  *      *        *        *   *         *   *  *      *             */
/*           *   *        *             *        *             *   *     *        *        *   *         *   *   *     *             */
/*                        *             *        * * *         *    *    *        *        *   *         *   *    *    *             */
/*          * * * *       *             *        *             *     *   *        *        *   *         *   *     *   *             */
/*                        *             *        *             *      *  *        *        *   *         *   *      *  *             */
/*         *       *      *             *        *             *       * *        *        *   *         *   *       * *             */
/*                        *             *        *             *        **        *        *   *         *   *        **             */
/*        *         *     *             *        * * * * * *   *         *        *        *   * * * * * *   *         *             */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        ATTENTION :                                                                                                                */
/*                                                                                                                                   */
/*                    Ce programme ('$xtc/verhulst.94$c') est                                                                        */
/*                  reference dans '$xiMd/verhulst.94$c.$m4' a                                                                       */
/*                  a des fins de demonstration 'WWW'.                                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Nota :                                                                                                                     */
/*                                                                                                                                   */
/*                    Il existe une version en "precision                                                                            */
/*                  infinie" de ce programme ('v $xtb/verhulst.94').                                                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/verhulst.94$c' :                                                                                           */
/*                                                                                                                                   */
/*                    John F. Colonna (LACTAMME, 1992MMJJhhmmss).                                                                    */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  < stdio.h >

#define   PRECIS    double                                                                                                              \
                                        /* Definition de la precision ('float' ou 'double') des calculs. Le phenomene de la          */ \
                                        /* sensibilite a l'ordre des operations ne se manifeste qu'en 'double'...                    */
#define   pas       10                                                                                                                  \
                                        /* Pas d'edition des valeurs (dx1,dx2,dx3,dx4,dx5).                                          */
#define   N         10                                                                                                                  \
                                        /* Nombre d'instants a editer (le nombre d'iterations est donc 'N*pas').                     */
#define   R         3.0                                                                                                                 \
                                        /* Taux de croissance.                                                                       */
#define   X0        0.5                                                                                                                 \
                                        /* Conditions initiales.                                                                     */

main()
     {
     PRECIS    r=R;
     PRECIS    dx1=X0;
     PRECIS    dx2=X0;
     PRECIS    dx3=X0;
     PRECIS    dx4=X0;
     PRECIS    dx5=X0;
     int       n;

     for  (n=0 ; n<(N*pas) ; n++)
          {
          if   ((n%pas) == 0)
               {
               printf("\n X(%04d) = %9.6f       %9.6f       %9.6f       %9.6f       %9.6f",n,dx1,dx2,dx3,dx4,dx5);
               }
          else
               {
               }

          dx1 = ((r+1)*dx1) - (r*(dx1*dx1));
          dx2 = ((r+1)*dx2) - ((r*dx2)*dx2);
          dx3 = (((r+1) - (r*dx3))*dx3);
          dx4 = (r*dx4) + ((1 - (r*dx4))*dx4);
          dx5 = dx5 + (r*((dx5) - (dx5*dx5)));
                                        /* Le resultat depend de la syntaxe...                                                       */
          }
     }



Copyright (c) Jean-François Colonna, 1995-2012.
Copyright (c) France Telecom R&D and CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 1995-2012.