In particular, I cannot understand what is the "matrix" that I need to provide as HISTORY in the input argument of any of the octave/dde-solvers -- i.e. any of the odeXXd solvers in the odepkg. There is no explanation on the specifics of this matrix and why it seems to be required to be of size n_var * n_delay. This is very different from MATLAB implementations. Any insight would be helpful.
Thanks in advance.
Maurizio