ikeda_lyap.h
/*-----------------------------------------------------------------
池田写像のリアプノフスペクトラムの計算
calculating the Lyapunov spectrum of the Ikeda map
x = 1.0 + MYU * ( x * cos(theta) - y * sin(theta) )
y = MYU * ( x * sin(theta) + y * cos(theta) )
theta = A - B/(double)( 1.0 + x * x + y * y )
ikeda_lyap.h
-----------------------------------------------------------------*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define ITERATION 1000000
#define SKIP 10000
// parameters
#define MYU 0.83
#define A 0.4
#define B 6.0
// 初期値
// initial values
#define X0 0.111
#define Y0 0.222
// 変数の数
// number of variables
#define DIM 2
/*------------------------------------------------------------------------*/
extern void ikeda_eqs( double (*x), double (*y), double u[DIM][DIM] );
extern void gram_schmidt_orth( double u[DIM][DIM], double e[DIM][DIM] );
extern double calc_lyap_dim( double lambda[DIM] );
extern double log2( double x );