I would like to give a theoretical estimation of local truncation error (and then for the global one) for a solution to a numerical initial value problem by NDSolve.
For this estimate I need to know what method has been used and the size of the steps (or an average size for the global error, since I think it uses adaptive size step). How can I get this information?