#second order linear constant coefficient
# to demonstrate canonical forms in phasespace 

# Define the differential equations
dx1/dt=a*x1+b*x2
dx2/dt=c*x1+d*x2

# Declare the parameters
par a=4,b=0,c=0,d=2

# Defining certain numbers
#numbers q=1,r=2

# Defining initial conditions (not necessary)
init x1=.1
init x2=.1

# Setting options  (not necessary)

@ maxstor=100000,xlo=-1,ylo=-1,xhi=1,yhi=1
@ xp=x1,yp=x2,zp=t,total=50,dt=.01
# xlo=-1,ylo=1,xhi=1,yhi=1
#zmin=4,zmax=36
#zhi=48,zlo=0
#xmin=-1,xmax=1,ymin=-1,ymax=1,
#@H0crv=0.01,Hmxcrv=0.05,Tint=6.28318530717959

#Tell xpp we are finished
done