I've been trying this in university:

ε = 0.005

p = 2

deq1 = x'[t] == (x[t] + y[t] - q*(x[t]^2) - x[t]*y[t])

deq2 = y'[t] == -y[t] + z[t] - x[t]*y[t]

deq3 = z'[t] ==   (x[t] - z[t])/p

tmin = 0; tmax = 20;

sol1 = NDSolve[{deq1, deq2, deq3, x[0] ==  8, y[0] ==   6, 
   z[0] ==   23}, {x, y, z}, {t, tmin, tmax}]

and it keeps giving me an error of

NDSolve::deqn: Equation or list of equations expected instead of True in the first
argument {(x'[t]==x[t]-q x[t]^2+y[t]-x[t] y[t],(y'[t]==-y[t]-x[t] y[t]+z[t],
(z'[t]==1/2 (x[t]-z[t]),True,True,True}.

I mean, I know what it says, but I have no idea how to fix it...Isn't this the way to write NDSolve? It's the first time I get something like that Please help