Sambine and Thompson Solver for DDE

This solver is used for Delayed Differential Equations. The code has been copied from the github repository without any change.

This method is also addressed in:
Thompson, S., and Shampine, L.F., Solver DDEs in Matlab, Applied Numerical Mathematics, 37, 441-458 (2001).
Control parameters for Sambine and Thompson DDE Solver
Name Descriptions Constraints and Notes
R tolerance Relative Error Tolerance Decreasing this value will increase accuracy, but add to the computation time
A Tolerance Absolute Error Tolerance Decreasing this value will increase accuracy, but add to the computation time
Max Steps Maximum number of steps Must be greater than 1
Init Step Size Initial step size Must be greater than 0
Max Step Size Maximum step Size Must be greater than 0 and less than whole time range