From 392a611d547287920c46edbdaa8c2025596c8edb Mon Sep 17 00:00:00 2001 From: bnewbold Date: Wed, 21 May 2014 00:59:10 -0400 Subject: hw02 first attempt --- assignments/hw02/euler.jl | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100644 assignments/hw02/euler.jl (limited to 'assignments/hw02/euler.jl') diff --git a/assignments/hw02/euler.jl b/assignments/hw02/euler.jl new file mode 100644 index 0000000..a9513bf --- /dev/null +++ b/assignments/hw02/euler.jl @@ -0,0 +1,29 @@ + +# Use Euler's method to integrate simple one variable ODE +# Use program as template for more interesting models + +# Note: produces slightly different output compared with +# example shown in class. Different values of b, dt, tlast + +using Winston + +a = 20 +b = 0.5 +c = 5 + +dt = 0.2 +tlast = 20 + +iterations = integer(round(tlast/dt)) +xall = zeros(Float64, (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)' +figure +plot(time,xall) -- cgit v1.2.3