/*Life is always simple, if you get going complex, you definitely are not living*/
/*Just see how patterns emerge from total chaos*/
/*enjoy, if you like maths or even drawing*/
#include <stdio.h>
main()
{
float x0, e, r, r1=3.987;
int i, T, n, M=100;
printf("Please enter the value of r(2,4] : ");
scanf("%f",&r);
printf("Please enter the number of time steps : ");
scanf("%d",&T);
float C[T][M];
printf("Please enter the value of x(0) (0,1] : ");
scanf("%f",&x0);
printf("Please enter the value of coupling strength, e : ");
scanf("%f",&e);
C[0][0]=x0;
for (i=1; i<=M-1; ++i)
C[0][i]=r1*C[0][i-1]*(1-C[0][i-1]);
for (n=1; n<=T-1; ++n)
{
C[n][0]= (1-e)*r*C[n-1][0]*(1-C[n-1][0]) + e*(r*C[n-1][M-1]*(1-C[n-1][M-1]) + r*C[n-1][1]*(1-C[n-1][1]))*0.5;
C[n][M-1]= (1-e)*r*C[n-1][M-1]*(1-C[n-1][M-1]) + e*(r*C[n-1][M-2]*(1-C[n-1][M-2]) + r*C[n-1][0]*(1-C[n-1][0]))*0.5;
for(i=1; i<=(M-2); i++)
C[n][i]= (1-e)*r*C[n-1][i]*(1-C[n-1][i]) + e*(r*C[n-1][i-1]*(1-C[n-1][i-1]) + r*C[n-1][i+1]*(1-C[n-1][i+1]))*0.5;
}
FILE *fp;
fp = fopen ("chaos.dat", "w");
int j, k;
for (j=0; j<=T-1; ++j)
{
fprintf (fp, "\n");
for (k=0; k<=(M-1); ++k)
fprintf (fp, "%f\n", C[j][k]);
}
fclose(fp);
printf("Now please open Chaos.dat and plot the data in gnuplot using Splot command (Help : google for the command syntax)\n Enjoy the patterns\n");
getch();
}
