/* 
   Simulation of the discrete logistic map.

   x[i+1] = h*x[i](1.0 - x[i])
   
   Last Modified: August 23, 1999 by Bill Goodwine 
*/

#include<stdio.h>
#include<stdlib.h>

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

  int i,steps;
  double x,h;
  char str[255];
  FILE *fp;

  if(argc != 4) { 
    printf("Usage: %s h xinit steps\n",*argv); 
    exit(1); 
  } 

  h = atof(argv[1]);
  x = atof(argv[2]);    
  steps = atoi(argv[3]);

  fp = fopen(".store.d","w");

  for(i=0;i<=steps;i++) {
    fprintf(fp,"%f\n",x);
    x = x*h*(1-x);
  }     

  fclose(fp);

  strcpy(&str[0],"gnuplot <<EOF
set terminal postscript eps
set output 'logistic.eps'
set nokey
set title 'h =");
  strcat(&str[0],argv[1]); 
  strcat(&str[0],"'\nplot '.store.d'\nEOF\n"); 

  system(str); 
  system("ghostview logistic.eps &");
  system("\\rm .store.d");
}