Hello,

Where in the MoL is the error estimates for RK45 calculated?  I can’t seem to find that in the source code.  Somehow they are written to the MOL::ERRORESTIMATE grid function but I can’t see where that is done.

-DG

-- 
David Garrison, Ph.D. 
Associate Dean for the College of Science and Engineering, 
Professor and Founding Chair of Physics
University of Houston-Clear Lake 
Bayou 3611
Houston, TX 77058 

Phone: 281-283-3796 
https://www.uhcl.edu/science-engineering/faculty/garrison-david 
http://www.uhcl.edu/physics

"If we knew what it was we were doing, it would not be called research, would it?" — Albert Einstein.