# Mathieu's equation ch3-mathieu.ode 
# First, write the 2nd-order van der Pol equation as a system
# of two first-order ODEs as follows:
dx/dt = y
dy/dt = -(kappa + 2*epsilon*cos(t))*x

# Set the default value of the parameter beta:
# Don't put spaces between the parameter name, the equal sign, and the parameter value!
par kappa=0.25,epsilon=0.05

# Set the default initial conditions:
init x=0.1  y=0.0

# Set the axes to plot x versus t as opposed to a phase plane: 
@ xplot=t, yplot=x

# Set the stopping time and the maximum allowable value of the dependent variables
@ total=200,bound=10000

# Set the default time step, the default viewing window, and the maximum number of data points to store.
@ dt=.01,xlo=0.0,xhi=200.0,ylo=-100.0,yhi=100.0,maxstor=100010
done