Max-plus algebra delay
[z] = maxplusdelay(x,d,T,A)
This function returns z, the matrix that contains propagation of delays until the delay is gone.
The model of the system is given by
x(k+1) = A⊗x(k)⊕d(k+1), x(0)=x,
where d(k) is the schedulling vector. The timetable d(k) is realistic if for all k≥0
A⊗d(k) ≤ d(k+1).
For detail see: J.G. Braker, "Algorithms and Applications Timed Discrete Event Systems”, Ph.D thesis, Department of Technical Mathematics an Informatics Delft University of Technology, (1993).