/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        E T U D E   D E   L ' A P P L I C A T I O N   L O G I S T I Q U E  :                                                       */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/logistik.02$c' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 1997MMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#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         4.0                                                                                                                 \
                                        /* Taux de croissance.                                                                       */
#define   X0        0.4                                                                                                                 \
                                        /* Conditions initiales.                                                                     */

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

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

          dx1 = r*(dx1*(1-dx1));
          dx2 = (r*dx2)-(r*(dx2*dx2));
          dx3 = (r*dx3)-((r*dx3)*dx3);
          dx4 = (r*dx4)-(dx4*(dx4*r));
                                        /* Application logistique...                                                                 */

          }
     }