W ramach projektu z modelowania systemów dynamicznych implementuje w C++ algorytm rozwiązywania równań różniczkowych Bogackiego–Shampine'a (ode23). Literatury niestety nie mam zbyt wiele dlatego jeżeli ktoś potrafiłby pomoc to proszę .
Mianowicie, na jakiej zasadzie powinno się dobrać wartość, o która pomniejszany będzie krok metody, w przypadku gdy błąd przekroczy dopuszczalna tolerancje, i ew jak też dobierać samą tolerancję?