-
Notifications
You must be signed in to change notification settings - Fork 0
/
userimg.jl
57 lines (52 loc) · 1.94 KB
/
userimg.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
using Plots
using DifferentialEquations
# cache things for the first plot
gr()
plot(1:4,1:4)
# cache things for ODE solvers
let
function f(du,u,p,t)
du .= 0
end
# u0 is 2-D array
u0 = [0.0 0.0 0.0]# Array{Float64,2}
solve( ODEProblem(f,u0,(0.0,1.0)) )
solve( ODEProblem(f,u0,(0.0,1.0)),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0)),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0)),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0)),AutoVern7(Rodas5()) )
# u0 is 1-D array
u0 = vec([0.0 0.0 0.0])# Array{Float64,1} / 1-D array
solve( ODEProblem(f,u0,(0.0,1.0)) )
solve( ODEProblem(f,u0,(0.0,1.0)),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0)),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0)),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0)),AutoVern7(Rodas5()) )
# add parameters
p = [1.0 2.0]
u0 = [0.0 0.0 0.0]
solve( ODEProblem(f,u0,(0.0,1.0),p) )
solve( ODEProblem(f,u0,(0.0,1.0),p),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0),p),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern7(Rodas5()) )
u0 = vec([0.0 0.0 0.0])
solve( ODEProblem(f,u0,(0.0,1.0),p) )
solve( ODEProblem(f,u0,(0.0,1.0),p),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0),p),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern7(Rodas5()) )
p = vec([1.0 2.0])
u0 = [0.0 0.0 0.0]
solve( ODEProblem(f,u0,(0.0,1.0),p) )
solve( ODEProblem(f,u0,(0.0,1.0),p),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0),p),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern7(Rodas5()) )
u0 = vec([0.0 0.0 0.0])
solve( ODEProblem(f,u0,(0.0,1.0),p) )
solve( ODEProblem(f,u0,(0.0,1.0),p),Rodas4() )
solve( ODEProblem(f,u0,(0.0,1.0),p),QNDF() )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern9(Rodas4P()) )
solve( ODEProblem(f,u0,(0.0,1.0),p),AutoVern7(Rodas5()) )
end