Question: Delay Differential Equations and Dsolve

Hello,
I'm working on coupled differential equation.
The first system is : 
y₁''+y₁'+y₁=q₁
q₁''+e(q₁²-1)q₁'+q₁=y₁"

And the second one : 
y₂''+y₂'+y₂=q₂
q₂''(t)+e(q₂²-1)q₂'(t)+q₂(t) = y₂"+f(P) q₁(t-tau)

This is a parametric system, f(P) and tau(P) are given function of the parameters P.
e is a constant

I have solved the first system with dsolve (using numeric option).
But when I try to solve the second one (with dsolve, numeric, setting P as a parameter), maple returns an error : 

"Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {t, t-1}"

I think Maple doesn't like " q1(t-tau)".
I have tried to create a new function q where :
q(t)=q1(t-tau)
But Maple returns the same error.

How can I fix it ?

Thanks for reading

EDIT : I have read there is no function in maple that solve delay differential equation.
But this is not a true DDE because q₂ has no effect on q₁. 
So I hope there is a way to "fool" maple and still use dsolve.

EDIT 2 :
I have found how to make it works.
I was using dsolve with the option compile (which increase (a lot) the efficiency of computation).
I delete this option and that's working.
Nevertheless, without the option compile, the computation is very very slow.
MapleHelp recommands to combine the 2 systems for more efficiency. But, when I combine, maple return the previous error.
How can I make it quicker ?

Here is the code : 

test_2_cylindre_sans_compile.mw