diff options
Diffstat (limited to 'notes/lec07_intro2_eulers')
-rw-r--r-- | notes/lec07_intro2_eulers | 46 |
1 files changed, 0 insertions, 46 deletions
diff --git a/notes/lec07_intro2_eulers b/notes/lec07_intro2_eulers deleted file mode 100644 index 17630d2..0000000 --- a/notes/lec07_intro2_eulers +++ /dev/null @@ -1,46 +0,0 @@ - -euler's method: - -dx/dt =~ ( x(t + Dt) - x(t) ) / Dt, for very small Dt - -so, x(t + Dt) = x(t) + f(x) * Dt, which is how to integrate the system - -"if Dt is too large, becomes highly unstable" (duh) - -scary MATLAB advice: -- watch out of totally crazy values (very high) -- non-negative values go negative - -MATLAB built-in ODE solvers: ode23, ode15s - -runge-kutta (use dxdt from between t_n and t_(n+1) ) -variable time-step methods - -summary: Euler's method sucks, news at 11. - -### misc julia notes - -the only place to find history (?!?!!) is in ~/.julia_history - -### codes - -using Winston - -a=20 -b=2 -c=5 -dt = 0.05 -tlast = 2 - -iterations = int( round(tlast/dt) ) -xall = zeros(iterations, 1) -x = c - -for i = 1:iterations - xall[i] = x - dxdt = a - b*x - x = x + dxdt*dt -end - -time = dt * [0:iterations-1]' -plot(time,xall) |