diff --git a/README.md b/README.md
index f91027b..3ce5658 100644
--- 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
new file mode 100644
index 0000000..5993ace
Binary files /dev/null and b/assets/README-66741291.png differ
diff --git a/prob2_ode.m b/prob2_ode.m
new file mode 100644
index 0000000..79e28f0
--- /dev/null
+++ b/prob2_ode.m
@@ -0,0 +1,3 @@
+function dy= prob2_ode(t,y)
+dy=-y;
+end
\ No newline at end of file
diff --git a/prob2ptB.m b/prob2ptB.m
new file mode 100644
index 0000000..4f336d6
--- /dev/null
+++ 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');
\ No newline at end of file