/////////////////////////////////////////////////////////////////////////////
//
// logistic.c  Logistic Map f(x)=mu x (1-x)
//
// Author: Evelyn Sander
//         Department of Mathematical Sciences
//         George Mason University
//         June 2009
//
/////////////////////////////////////////////////////////////////////////////
//#define DEBUG

#include "auto_f2c.h"

/////////////////////////////////////////////////////////////////////////////
//
// Required by AUTO. Calculates RHS of MAP.
//
/////////////////////////////////////////////////////////////////////////////

int func (integer ndim, const doublereal *u, const integer *icp,
          const doublereal *par, integer ijac,
          doublereal *f, doublereal *dfdu, doublereal *dfdp)

{
  int i;
  double temp;
  temp=u[0];
  for(i=0;i<24; i++){ 
    temp=par[0]*temp*(1-temp);
  }
  f[0]=temp-u[0];
  return(0);
}

/////////////////////////////////////////////////////////////////////////////
//
// Required by AUTO.
// Input:
//    ndim: Dimension of the ODE system
// Output:
//    u: A starting solution vector
//  par: The corresponding equation-parameter values
//
/////////////////////////////////////////////////////////////////////////////

int stpnt (integer ndim, doublereal t, doublereal *u, doublereal *par) 
{  
   par[0] = 3.63;  // mu
   u[0]=0.505947826;
   return 0;
}

/* ---------------------------------------------------------------------- */
/* ---------------------------------------------------------------------- */
int pvls (integer ndim, const doublereal *u,
          doublereal *par)
{

  return 0;
}
/* ---------------------------------------------------------------------- */
/* ---------------------------------------------------------------------- */
int bcnd (integer ndim, const doublereal *par, const integer *icp,
          integer nbc, const doublereal *u0, const doublereal *u1, integer ijac,
          doublereal *fb, doublereal *dbc)
{
  return 0;
}
/* ---------------------------------------------------------------------- */
/* ---------------------------------------------------------------------- */
int icnd (integer ndim, const doublereal *par, const integer *icp,
          integer nint, const doublereal *u, const doublereal *uold,
          const doublereal *udot, const doublereal *upold, integer ijac,
          doublereal *fi, doublereal *dint)
{
  return 0;
}
/* ---------------------------------------------------------------------- */
/* ---------------------------------------------------------------------- */
int fopt (integer ndim, const doublereal *u, const integer *icp,
          const doublereal *par, integer ijac,
          doublereal *fs, doublereal *dfdu, doublereal *dfdp)
{
  return 0;
}
/* ---------------------------------------------------------------------- */
/* ---------------------------------------------------------------------- */