problem 2 part B

nin13001 · Dec 6, 2017 · fdd7e99 · fdd7e99
1 parent 40a7ed0
commit fdd7e99
Show file tree

Hide file tree

Showing 4 changed files with 32 additions and 1 deletion.
diff --git a/README.md b/README.md
@@ -1,4 +1,21 @@
 # 06_initial_value_ode
 # Problem 2
 ## Part A
-The analytical solution to y' = -y is y = e^-t 
+The analytical solution to y' = -y is y = e^-t
+
+## Part B
+
+```
+dt=[0.5 .1 .001];
+figure();
+hold on
+     x=[0:dt(2):3]';
+     y_n=zeros(size(x));
+    y_n(1)=1;
+for i=2:length(x)
+    dy=prob2_ode(x,y_n(i-1));
+    y_n(i)=y_n(i-1)+ dt(2)*dy;
+end  
+plot(x,exp(-x),'k-',x,y_n,'o');
+```
+![](assets/README-66741291.png)
diff --git a/assets/README-66741291.png b/assets/README-66741291.png
diff --git a/prob2_ode.m b/prob2_ode.m
@@ -0,0 +1,3 @@
+function dy= prob2_ode(t,y)
+dy=-y;
+end
diff --git a/prob2ptB.m b/prob2ptB.m
@@ -0,0 +1,11 @@
+dt=[0.5 .1 .001];
+figure();
+hold on
+     x=[0:dt(2):3]';
+     y_n=zeros(size(x));
+    y_n(1)=1;
+for i=2:length(x)
+    dy=prob2_ode(x,y_n(i-1));
+    y_n(i)=y_n(i-1)+ dt(2)*dy;
+end  
+plot(x,exp(-x),'k-',x,y_n,'o');