how to write fractional system with delay in matlab?

Hello, I want to write this sytem fractional with delay in matlab?

D^{alpha } y_{1}(t) = 10.5 - y_{1} (t) / 1 + 0.0005 y_{4}^{3}(t-4)

D^{alpha} y_{2}(t) = (y_{1}(t) / 1 + 0.0005 y_{4}^{3} (t-4) )-y_{2}(t),

D^{alpha} y_{3} (t) = y_{2}(t) - y_{3}(t)

D^{alpha} y_{4} (t) = y_{3} (t) - 0.5 y_{4}(t)

y(t) = [60, 10,10,20]^{T}, t \leq 0

alpha= 0.95

How to plot the above system?