> 
restart;
# Call Maple DEtools Package:
>  with(DEtools):
------------------------------
# Define System of 2ODEs:
> ODE2:=[diff(x(t),t)= 0.25*x*(1-x/1.4)-1.5*x*y/(1+0.16*x),\
diff(y(t),t)=0.9*1.5*x*y/(1+0.16*x)-0.8*y]:
------------------------------
# DEplot Format:
------------------------------
# DEplot([rhs of eqns],[vars],{[i.c]},stepsize=h,
#   method=`Euler` or default is RK4)
# Remark:  Press Enter to Generate the DEplot Plots below:
------------------------------
#  DEplot resulting in 3D plot in txy-space by default RK4 method:
> DEplot(ODE2,[t,x,y],t=0..50,{[0,1,2]},stepsize=0.1,\
title=`DEplot with default RK4`);
------------------------------
------------------------------
# DEplot by Euler method:
> DEplot(ODE2,[t,x,y],t=0..50,{[0,1,2]},
stepsize=0.1,\
method=`Euler`,title=`DEplot with Euler`);
------------------------------
# DEplot by Modified or Improved Euler method:
> DEplot(ODE2,[t,x,y],t=0..50,{[0,1,2]},stepsize=0.1,\
method=`impeuler`,title=`DEplot with Modified Euler`);
------------------------------