1) First of all, can you describe in a  bit more detailed way the usage of AppCtx class of matfree.py module to solve ODE systems <br>(determined as in rober.py), without jacobian initialisation, in other words how can change rober.py to solve this issue?<br>
<br>2) Does THETA integration implement time step adaptation?<br><br>3) Suppose i have a large ODE system, how can i implement multiprocessor (parallel) integration in a way similar with those (function definition and plotting) in rober.py?<br>