mackey-glass_DDE_lyap.h

/*-----------------------------------------------------------------

マッキー・グラス遅延微分方程式

Mackey-Glass delay differential equation

dx/dt = A * x(t-tau)/( 1 + x(t-tau)^C ) - B * x(t)

mackey-glass_DDE_lyap.h

-----------------------------------------------------------------*/

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#define ITERATION 50000

#define SKIP      20000

#define DIM 100 /* 想定した変数の数 (assumed number of the variables) */

#define X_TAU   (DIM-1)

#define DELTA_T (TAU/(double)DIM) /* step size of the fourth-order Runge-Kutta method */

/*------------------------------------------------------------------------*/

#define TAU  17.0

#define A     0.2

#define B     0.1

#define C    10.0

/*------------------------------------------------------------------------*/

extern void mackey_glass_DDE( double x, double delay, double (*dxdt) );

extern void solve_mackey_glass_DDE_using_RK4( double x[DIM], double (*dxdt), double u[DIM][DIM] );

extern void quick_sort( int n, double x[DIM] );

extern void calc_lyap_dim( double lambda[DIM], double (*lyap_dim), int (*iii) );

extern void gram_schmidt_orth( double u[DIM][DIM], double e[DIM][DIM] );

extern double log2( double x );