Смотреть видео: 2.5.3. Решение оду с помощью функции odesolve(?