/*-----------------------------------------------------------------
ティンカーベル写像のリアプノフスペクトラムの計算
calculating the Lyapunov spectrum of the Tinkerbell map
x[t+1] = x[t] * x[t] - y[t] * y[t] + A * x[t] + B * y[t]
y[t+1] = 2.0 * x[t] * y[t] + C * x[t] + D * y[t]
tinkerbell_lyap.h
-----------------------------------------------------------------*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define ITERATION 100000
#define SKIP 10000
// parameters
#define A 0.9
#define B -0.6
#define C 2.0
#define D 0.5
// 初期値
// initial values
#define X0 0.164
#define Y0 0.124
// 変数の数
// number of variables
#define DIM 2
/*------------------------------------------------------------------------*/
extern void tinkerbell_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 );