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

#include "pngwriter.h"

#define nx 15000
#define ny 9000

int main()
{

time_t seconds1;             //Time counter
time_t seconds2;
seconds1 = time (NULL);
srand(seconds1);

//*************** FILES *****************
FILE * f1;
FILE * f2;
FILE * susce;

f1 = fopen("f1.txt","w");
f2 = fopen("f2.txt","w");
susce =  fopen("susce.txt","w");

pngwriter png(nx,ny,0,"bifu2.png");

//***************************************

int jx,jy;
double x,r;
double r1 = 3.4;
double r2 = 4.0;
double x1 = 0.0;
double x2 = 1.0;

for(jx = 0; jx < nx; jx++)
	{
	for(jy=0;jy<ny;jy++) png.plot(jx,jy,1.0,1.0,1.0);

	r =  r1 + 1.0*(r2-r1)*jx/nx;
	x = 0.25;
	
	int n = 0;
	while(n < 1000)
	{
	    x = r*x*(1.0-x);
		n++;
	}
	while(n < 2500)
	{
	    x = r*x*(1.0-x);
		jy = floor(ny*x);
		png.plot(jx,jy, 0.0,0.0,0.0);
		n++;
	}
    }
//printf("%f %f\n",a,b);

seconds2 = time (NULL);

printf("\n\nTime spent: %d seconds = %f minutes\n", (seconds2-seconds1),(seconds2-seconds1)/60.0);

fclose(f1); 
fclose(f2);
fclose(susce);

png.close();


return(0);

}