/************************************************************
*  Henon.h
*  說明：各種Henon映射
*  2004.11.12 zya
***********************************************************/

//Henon吸引子函數
void henon(double *X, double *Y, double x0, double y0, int len)
{
	int i;
	double a=1.4;
	double b=0.3;
	X[0]=x0;
	Y[0]=y0;

	for (i=1; i<len; i++)
	{
		X[i] = 1+b*Y[i-1]-a*X[i-1]*X[i-1];
		Y[i] = X[i-1];
	}
}

void hyperhenon(double *X, double *Y, double *Z, double x0, double y0, double z0,int len)
{
	int i;
	double a=1.76;
	double b=0.1;
	X[0]=x0;
	Y[0]=y0;
	Z[0]=z0;

	for (i=1; i<len; i++)
	{
		X[i] = a-Y[i-1]*Y[i-1]-b*Z[i-1];
		Y[i] = X[i-1];
		Z[i] = Y[i-1];
	}
}