# Time to reach steady-state (ODEProblem)

**URL:** https://discourse.julialang.org/t/time-to-reach-steady-state-odeproblem/81501
**Category:** New to Julia
**Created:** [May 23, 2022, 10:23am UTC](https://discourse.julialang.org/t/time-to-reach-steady-state-odeproblem/81501 "2022-05-23T10:23:59Z")
**Posts on this page:** 2
**Page:** 1

<div class="post-metadata">

### Author: ![SiddhantDeshmukh](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/siddhantdeshmukh/32/35324_2.png) [@SiddhantDeshmukh](https://discourse.julialang.org/u/SiddhantDeshmukh)
#### Post date: [May 23, 2022, 10:23am UTC](https://discourse.julialang.org/t/time-to-reach-steady-state-odeproblem/81501/1 "2022-05-23T10:23:59Z")

</div>

Hello,  
I would like to find the point at which all variables in a given ODESystem reach a steady-state, i.e. their relative change is below some defined threshold. I came across SteadyStateProblem to calculate a steady-state, but I was not able to get this working with modelingtoolkitize() to automatically compute a Jacobian matrix. In my problem setup, I solve a chemical kinetics problem with different initial conditions many times, and I would further like to find the point at which various species reach a steady-state. Is it possible to / does it make sense to calculate the time at which species reach a steady-state?

---

<div class="post-metadata">

### Author: ![Wen-Wei\_Tseng](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/wen-wei_tseng/32/28839_2.png) [@Wen-Wei\_Tseng](https://discourse.julialang.org/u/Wen-Wei_Tseng)
#### Post date: [May 24, 2022, 8:26am UTC](https://discourse.julialang.org/t/time-to-reach-steady-state-odeproblem/81501/2 "2022-05-24T08:26:52Z")

</div>

[`TerminateSteadyState`](https://diffeq.sciml.ai/stable/features/callback_library/#TerminateSteadyState) callback might work for your case.