share/contrib/Eulix/Eulix_T3N.mac

   1 load("Eulix.mac")$
   2
   3
   4 /* Lotka-Volterra
   5     alpha= 0.1  beta= 0.015  gamma= 0.0225  delta= 0.02
   6
   7     y1' =  y1*(alpha-beta*y2)     y1(0)= 1
   8     y2' = -y2*(gamma-delta*y1)    y2(0)= 1
   9 */
  10
  11 alpha:0.1$beta:0.015$gamma:0.0225$delta:0.02$
  12
  13
  14 y1_198:1.01243149501057267529$
  15 y2_198:0.99965472408223507814$
  16
  17 debugmode(true)$
  18 ER: Eulix([y1*(alpha-beta*y2),-y2*(gamma-delta*y1)],[y1,y2],[1,1],['t,0,198],'tabular='none)$
  19 print("solution at final time: ",ER)$
  20