Mercurial > repos > public > sbplib
comparison +time/Rk4SecondOrderNonlin.m @ 30:cd2e28c5ecd2
Corrected parameter treatment in nonlin timesteppers.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 28 Sep 2015 13:23:52 +0200 |
| parents | 16acb2775aca |
| children | b5e5b195da1e c6fcee3fcf1b |
comparison
equal
deleted
inserted
replaced
| 29:32b39dc44474 | 30:cd2e28c5ecd2 |
|---|---|
| 17 methods | 17 methods |
| 18 function obj = Rk4SecondOrderNonlin(D, E, S, k, t0, v0, v0t) | 18 function obj = Rk4SecondOrderNonlin(D, E, S, k, t0, v0, v0t) |
| 19 default_arg('S',0); | 19 default_arg('S',0); |
| 20 default_arg('E',0); | 20 default_arg('E',0); |
| 21 | 21 |
| 22 if isnumeric(S) && S == 0 | 22 if isnumeric(S) |
| 23 S = @(v,t)0; | 23 S = @(v,t)S; |
| 24 end | 24 end |
| 25 | 25 |
| 26 if isnumeric(E) && E == 0 | 26 if isnumeric(E) |
| 27 E = @(v,t)0; | 27 E = @(v)E; |
| 28 end | 28 end |
| 29 | 29 |
| 30 obj.k = k; | 30 obj.k = k; |
| 31 obj.t = t0; | 31 obj.t = t0; |
| 32 obj.w = [v0; v0t]; | 32 obj.w = [v0; v0t]; |
