I am trying to write a program using Maple software,
the sys of delay DEs are simplified in the form
x'=constant * x(t-tau)+nonlinear function1 of x(t) and, y(t)
y'=constant* y(t-tau)+ nonlinear function2 of x(t)and y(t)
t is the time from 0 to t_N, x'=dx/dt, tau is the delay, ICs are x(0)=y(0)=0
Mathematica did it and worked fine after the assumption
x(-tau->0)=y(-tau->0)=t^2 worked. Maple's company said they
have no subroutine to do it.
Any help will be highly appreciated