/* ***********************************************************************
   *									 *
   *				Bifurcation				 *
   *									 *
   *				Program by				 *
   *			   Christopher D. Watkins			 *
   *									 *
   *			     'C' conversion by				 *
   *                            Larry Sharp				 *
   *									 *
   ***********************************************************************
*/

#include "stdio.h"
#include "dos.h"
#include "conio.h"
#include "math.h"
#include "math.inc"
#include "graph2.inc"

#define min  0.0
#define max  125.0
#define Res  300.0
#define n    300
#define beta 5
#define x0   0.9
#define T    25

float dlambda;
float lambda;
float scale, x;
int   i, Temp;
int   xxx, yyy;
char k;

void main()
{
  InitGraphics(56);
  dlambda=(max-min)/Res;
  scale=Res/max;
  lambda=min;

  do
  {
    x=x0;
    for(i=1; i<=n; i++)
    {
      x=lambda*x/Power((1+x), (beta));
      xxx=Round(30*x);
      yyy=Round(lambda*scale);
      Temp=GetPixel(xxx, yyy);
      Temp+=1;
      if(Temp<T)
	Plot(xxx, yyy, Temp+128);
      else
      {
	Plot(xxx, yyy, 0);
	Temp-=1;
      }
    }
    lambda+=dlambda;
  }
  while(lambda<=max);

  ExitGraphics();
}