Parabolic differential equations with time delay

Let $d_1=1,d_2=2,a_{11}=frac{5}{13},a_{12}=frac{22}3,a_{21}=-2,a_{22}=frac{6}7,tau=frac{5}7$, $psi(t,x)=cos^42x,phi(t,x)=frac{3}{13}x^4sin^2 3x$, $Omega=[0,200]$

How to solve:

$$left{
begin{array}{lc}
dfrac{partial u(t,x)}{partial t}=d_1triangle u(t,x)+u(t,x)left(r_1-a_{11}u(t-tau,x)-a_{12}v(t,x)right),& t>0,xinOmega
dfrac{partial v(t,x)}{partial t}=d_2triangle v(t,x)+v(t,x)left(-r_2+a_{21}u(t,x)-a_{22}v(t,x)right),& t>0,xinOmega
dfrac{partial u}{partial n}=dfrac{partial v}{partial n}=0,quad tge0,xinpartialOmega quad(text{Neumann conditions})
u(t,x)=phi(t,x)ge 0,qquad v(t,x)=psi(t,x)ge 0, &(t,x)in[-tau,0]timesOmega
end{array}
right.$$