Problems with Initial Solutions

[drsaroj]

I have a query as to why the jitcdde solution is different than obtained by Mathematica NDSolve when it applies to the same problem with the same set of parameters and the same initial history function. 1a.jpeg and 1b.jpeg appear to be the same except the initial solution abound time=0.


python.txt

[Wrzlprmft]

The reason is that step_on_initial continuities integrates the dynamics up to t=10. You then start integrating again from 0 to 10 (instead of starting at 0). Therefore the solutions you get here are extrapolated from the last step (i.e., from 10−ε to 10). (And yes, this should raise a warning.) If you really need the beginning, you have to use the output of several steps of integrate_blindly.