# Fixedpoints DynamicalSystems.jl

**URL:** <https://discourse.julialang.org/t/fixedpoints-dynamicalsystems-jl/96004>\
**Category:** Modelling & Simulations\
**Created:** [March 13, 2023, 1:20pm UTC](https://discourse.julialang.org/t/fixedpoints-dynamicalsystems-jl/96004 "2023-03-13T13:20:46Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![Sergey\_Novak](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/sergey_novak/32/37716_2.png) [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Post date:** [March 13, 2023, 1:20pm UTC](https://discourse.julialang.org/t/fixedpoints-dynamicalsystems-jl/96004/1 "2023-03-13T13:20:46Z")

</div>

Hello, how can i automaticaly calculate Jacobian?  
I have next system and manual calculate Jacobian is hard  
 ![image](https://global.discourse-cdn.com/julialang/original/3X/3/d/3dd803404023a63816d8c6b26147484b5c5e72be.png)  
Code:

```julia
function TM(u, p, t)
    U(y) = p[8] + p[9] / ( 1.0 + exp( -50.0 * (y - p[7]) ) )
    σ(x) = 1.0 / ( 1.0 + exp( -20.0 * (x-p[6]) ) )
    du1 = (-u[1] + p[1] * log( 1.0 + exp( (p[5] * U(u[3]) * u[2] * u[1] + p[11] ) / (p[1]) ) ) ) / p[2]
    du2 = (1.0 - u[2])/p[3] - U(u[3])*u[2]*u[1]
    du3 = (-u[3])/p[4] + p[10] * σ(u[2])
    return SA[du1, du2, du3]
end

t = 5000.0; tstep = 0.001; ttr = 100000.0
integ_set = (alg = RK4(), adaptive = false, dt = tstep)

const τ = 0.013; const τD = 0.080; const τy = 3.3; const J = 3.07; const β = 0.300
const xthr = 0.75; const ythr = 0.4
const α = 1.58; const ΔU0 = 0.305;
I0_local = -1.6102257124
p = SA[α, τ, τD, τy, J, xthr, ythr, U0_, ΔU0, β, I0_local]

ds_A1 = CoupledODEs(TM, SA[3.614150444316901, 0.7995140316515059, 0.4256000782173303], p, diffeq = integ_set)
tr_A1, _ = trajectory(ds_A1, t, Δt = tstep; Ttr = ttr)

fp, eigs, stable = fixedpoints(ds_A1, box)

```

---

<div class="post-metadata">

**Author:** ![Datseris](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/datseris/32/13406_2.png) [@Datseris](https://discourse.julialang.org/u/Datseris)\
**Post date:** [March 13, 2023, 1:56pm UTC](https://discourse.julialang.org/t/fixedpoints-dynamicalsystems-jl/96004/2 "2023-03-13T13:56:29Z")

</div>

[https://symbolics.juliasymbolics.org/dev/](https://symbolics.juliasymbolics.org/dev/)

---

<div class="post-metadata">

**Author:** ![Sergey\_Novak](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/sergey_novak/32/37716_2.png) [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Post date:** [March 13, 2023, 2:05pm UTC](https://discourse.julialang.org/t/fixedpoints-dynamicalsystems-jl/96004/3 "2023-03-13T14:05:47Z")

</div>

Thank you. Why is general information not displayed when I write `println(ds_A1) `? Am I doing something wrong? For example, I want jupyter notebook to display information about `ds_A1` and `ds_A2` in the one cell