diff options
Diffstat (limited to 'hw3/Lorenz.jl')
| -rw-r--r-- | hw3/Lorenz.jl | 13 |
1 files changed, 4 insertions, 9 deletions
diff --git a/hw3/Lorenz.jl b/hw3/Lorenz.jl index 4a7fa0e..b1dd8f0 100644 --- a/hw3/Lorenz.jl +++ b/hw3/Lorenz.jl @@ -1,3 +1,6 @@ +# FOR PROBLEM 3.25 +# author: sotech117 + using DifferentialEquations using Plots @@ -25,8 +28,7 @@ r_maxes = [] z_maxes = [] r_mins = [] z_mins = [] -let - r_diverge_1 = -1 +let for i in 1:r_steps r = r_values[i] p = [10.0, r, 8/3] # parameters of the Lorentz 63 system @@ -49,12 +51,6 @@ let push!(r_mins, r_values[i]) push!(z_mins, z_values[j]) end - - # calcualte dz - # if r_diverge_1 < 0 && length(z_mins) > 1 && length(z_maxes) > 1 && abs(z_mins[end] - z_maxes[end]) > 20.0 - # r_diverge_1 = r - # println("Divergence starts at r = $r") - # end end if i % 500 == 0 @@ -70,6 +66,5 @@ let vline!([122], label="Critical r value ≈ 122", lc=:orange, lw=1.5, ls=:dash) vline!([146], label="Critical r value ≈ 146", lc=:purple, lw=1.5, ls=:dash) - savefig("hw3/test11.png") end
\ No newline at end of file |