mg - plot the phase space of the Mackey-Glass system
mg [-width integer] [-height integer] [-skip integer]
[-points integer] [-delta integer] [-tau integer]
[-A double] [-B double] [-dt double] [-x0 double]
[-factor double] [-data] [-inv] [-mag integer]
The phase space of the Mackey-Glass system, which is
described by the delay differential equation
dx(t)/dt = A * x(t-Tau) / (1 + x(t-Tau)^10) - B * x(t),
is plotted according to the specified parameters. The x-
coordinate of the plot is determined by x(t) while the y-
coordinate is determined by x(t-delta).
Width of the plot in pixels.
Height of the plot in pixels.
Number of initial points to skip.
Number of points to plot.
Time steps to delay for.
Value of the Tau parameter.
Value of the A parameter.
Value of the B parameter.
Time step size.
Initial X value.
Auto-scale expansion factor.
-data Don't plot, but print points.
-inv Invert all colors?
How to plot points.
The plot region is determined by the points that are ini-
tially skipped. If this number is too small (i.e., it is
not very representative of the range of the plotted val-
ues), then you may need to increase the number specified
by the -skip option. Alternatively, you can adjust the
value given to -factor, which simply fractionally
increases the border of the plot.
The program uses a second-order Euler's method to perform
the numerical integration, which is sufficient for simple
tasks such as this.
No sanity checks are performed to make sure that any of
the options make sense.
Copyright (c) 1997, Gary William Flake.
Permission granted for any use according to the standard
GNU ``copyleft'' agreement provided that the author's com-
ments are neither modified nor removed. No warranty is
given or implied.