//
// Programmer:    Craig Stuart Sapp <craig@ccrma.stanford.edu>
// Creation Date: Fri May 19 18:54:58 PDT 2000
// Last Modified: Fri May 19 18:55:01 PDT 2000
// Filename:      ...sig/doc/examples/all/Restrepotune/Restrepo.cpp
// Syntax:        C++ 
// $Smake:        g++ -O -o %b %f && strip %b
//
// Description:   This program generates a sequence of values
//                based on the 2-D invertible iterated map with
//                chaotic solutions proposed by Alfredo Restrepo
//                <arestrep@uniandes.edu.co>.
//
//                The equations for generating the sequence are:
//
//                x[n+1] = x[n](x[n]-1)/x[n-1]
//                and
//                y[n] = x[n]
// 
//                Where y[n] is the current output and x[n] is the
//                current input.
//

#include <iostream.h>
#include <stdlib.h>

void generateRestrepoValues(int count, double x = 1.99, double oldx = 1.99);

int main(int argc, char* argv[]) {

   if (argc < 2) {
      cout << "Usage: "<< argv[0] << " repetitions x[0] x[-1]" << endl;
      cout << "   Suggested values: \n";
      cout << "      x[0]  = 1.99\n";
      cout << "      x[-1] = 1.99\n";
      cout << endl;
      exit(1);
   }
   
   int count = atoi(argv[1]);

   double x = 1.99;
   double oldx = 1.99;

   if (argc >= 3) {
      x = strtod(argv[2], NULL);
   }
   if (argc >= 4) {
      oldx = strtod(argv[3], NULL);
   }

   generateRestrepoValues(count, x, oldx);

   return 0;
}

///////////////////////////////////////////////////////////////////////////

void generateRestrepoValues(int count, double x, double oldx) {
   double newx;
   cout << oldx << '\n';
   for (int i=0; i<count; i++) {
      cout << x << '\n';
      newx = x * (x - 1) / oldx;
      oldx = x;
      x = newx;
   }
}


/*

Date: Thu, 18 May 2000 19:12:46 -0500 (GMT+5)
From: Alfredo Restrepo <arestrep@uniandes.edu.co>
To: Juan Reyes <juanig@ccrma.stanford.edu>
Subject: Re: sugerencias

Hola Juan,

de pronto pueden usar este mapa:

x(n+1) = x(n) [x(n) - 1] / x(n-1)

La idea es darle los valores iniciales x(0) y x(1) e iterar de ahi en adelante.
Este mapa tambien se puede ver en dimension 2, como el de Henon.

No todos los valores sirve, por ejemplo, no usen ni cero ni uno. Algunos
que sirven son:
1.99, 1.99
0.035, 1.1

Esta ecuacion sale del estudio de las senales raiz del filtro Teager y son
resultados preliminares.
De pronto si usan la ecuacion en un paper, nos podrian referenciar.

Si tienes mas preguntas, a la orden.

Saludos,
Alfredo                                     

*/