# Setting parameters for \`ReachabilityAnalysis\`'s \`@ivp\`

**URL:** <https://discourse.julialang.org/t/setting-parameters-for-reachabilityanalysis-s-ivp/67776>\
**Category:** Numerics\
**Tags:** question, reachabilityanalysis, reachability\
**Created:** [September 6, 2021, 9:52pm UTC](https://discourse.julialang.org/t/setting-parameters-for-reachabilityanalysis-s-ivp/67776 "2021-09-06T21:52:53Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![cadojo](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/cadojo/32/25328_2.png) [@cadojo](https://discourse.julialang.org/u/cadojo)\
**Post date:** [September 6, 2021, 9:52pm UTC](https://discourse.julialang.org/t/setting-parameters-for-reachabilityanalysis-s-ivp/67776/1 "2021-09-06T21:52:53Z")

</div>

I see in documentation for `ReachabilityAnalysis` that parameters are typically set in a function, as is the case in the [duffing oscillator example](https://juliareach.github.io/ReachabilityAnalysis.jl/dev/models/DuffingOscillator/#Duffing-Oscillator). I’m working with a system where the `ODEFunction` was automatically generated with `ModelingToolkit`, so I didn’t “write” my `ODEFunction`. Can I somehow set the parameters in the `@ivp` with something like `@ivp(..., p=[1,2,3])`?

---

<div class="post-metadata">

**Author:** ![mforets](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mforets/32/298_2.png) [@mforets](https://discourse.julialang.org/u/mforets)\
**Post date:** [September 7, 2021, 12:53am UTC](https://discourse.julialang.org/t/setting-parameters-for-reachabilityanalysis-s-ivp/67776/2 "2021-09-07T00:53:54Z")

</div>

That functionality that you mention, about hooking MTK generated functions with JuliaReach, is not implemented.

The `@ivp` macro cannot handle the params arg as you tried (seems like a good idea). But it is possible to pass params vectors through the solve function like this:

```julia
@taylorize function duffing2!(du, u, p, t)
    α, β, γ, δ = p
    x, v = u
    f = γ * cos(ω * t)
    du[1] = u[2]
    du[2] = - α*x - δ*v - β*x^3 + f
end
# ...
sol = solve(prob, T=10.0, alg=TMJets21a(), params=(-1.0, 1.0, 0.3, 0.37))

```