Mercurial > repos > public > sbplib_julia

comparison TimeStepper.jl @ 12:c4bc1165fdf7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
remove white space
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 17 Dec 2018 14:47:24 +0100
parents bed51234616b
children 1e845cd91cd3
comparison
equal deleted inserted replaced
11:7f075bacbd68 12:c4bc1165fdf7
1 abstract type TimeStepper end 1 abstract type TimeStepper end
2 2
3 # Returns v and t 3 # Returns v and t
4 function getState(ts::TimeStepper) 4 function getState(ts::TimeStepper)
5 error("not implemented") 5 error("not implemented")
6 end 6 end
49 k1 = ts.F(ts.v,ts.t) 49 k1 = ts.F(ts.v,ts.t)
50 k2 = ts.F(ts.v+0.5*ts.k*k1,ts.t+0.5*ts.k) 50 k2 = ts.F(ts.v+0.5*ts.k*k1,ts.t+0.5*ts.k)
51 k3 = ts.F(ts.v+0.5*ts.k*k2,ts.t+0.5*ts.k) 51 k3 = ts.F(ts.v+0.5*ts.k*k2,ts.t+0.5*ts.k)
52 k4 = ts.F(ts.v+ ts.k*k3,ts.t+ ts.k) 52 k4 = ts.F(ts.v+ ts.k*k3,ts.t+ ts.k)
53 ts.v = ts.v + (1/6)*(k1+2*(k2+k3)+k4)*ts.k 53 ts.v = ts.v + (1/6)*(k1+2*(k2+k3)+k4)*ts.k
54 54
55 ts.n = ts.n + 1 55 ts.n = ts.n + 1
56 ts.t = ts.t + ts.k 56 ts.t = ts.t + ts.k
57 57
58 return nothing 58 return nothing
59 end 59 end