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.9

#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 );