1 files changed, 29 insertions, 0 deletions
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)