#---------------------------------------------------------
#
#     Sample code
#                ----------- The toda A2 equation
#
# Note: Need about 10 minutes to solve this equation
#
#---------------------------------------------------------

function list {
	ROOT3_2 = 0.8660254037844386
	TOdy1dt(x,y,z,w)= z
	TOdy2dt(x,y,z,w)= w
	TOdy3dt(x,y,z,w)=
	 	-(exp(x)-0.5*(exp(-0.5 *x+ROOT3_2*y)+exp(-0.5*x-ROOT3_2*y)))
	TOdy4dt(x,y,z,w)=
		-ROOT3_2*(exp(-0.5*x+ROOT3_2*y)-exp(-0.5*x-ROOT3_2*y)) 
	}

toda2 = eqn{
         [TOdy1dt(x,y,z,w),TOdy2dt(x,y,z,w),TOdy3dt(x,y,z,w),TOdy4dt(x,y,z,w)]
	 [initials = 0.9:-1.01:52.02:-1.03]
	 [time = 0.0:100.0]  [step = 0.005] [skip 5]
	 [method=rkqc]
	 [small = 0.00000001]
	 [scale = 1.0,1.0,20.0,20.0] 
         [axes: automatic]
	 };
#---------------------------------------------------------
center;
plot toda2;
#---------------------------------------------------------

