/*
Henon Map
X_{k+1}=Y_k+1-a X^2_k
Y_{k+1}=b X_k
*/
#include
#include
#include
#define a 1.5
#define b 0.25
main()
{
int j; /*j;step number*/
double x=1.0,nx=0.0;
double y=0.0,ny=0.0;
for(j=0;j<=100000;j++)
{
nx=y+1.0-a*x*x;
ny=b*x;
printf("%f %f\n",nx, x);
/*printf("%d %f\n",j, nx);*/
x=nx;
y=ny;
}
return 0;
}
出力結果
ZOOM
ZOOM ZOOM
Back to C Language
