diff --git a/final_project/README.md b/final_project/README.md
index cee3ee3..7a7b0ac 100644
--- a/final_project/README.md
+++ b/final_project/README.md
@@ -56,7 +56,7 @@ b. Use a Monte Carlo model to determine the mean and standard deviation  for the
 maximum deflection $\delta x$ if b and h are normally distributed random variables
 with 0.1 % standard deviations at q=50 N/m. 
 
-3. Now use the central difference approximation to set up a system of equations for the
+2. Now use the central difference approximation to set up a system of equations for the
 beam for q(x)=cst, P=0, and $\omega=0$. Use the boundary conditions with a numerical
 differentiation to determine the valuea of the end points
 
@@ -71,7 +71,7 @@ differentiation to determine the valuea of the end points
     e. Comment on the results from the analytical and numerical approaches (if you used
     functions then provide help files, if you used scripts, then describe the steps used)
 
-4. Now set up the system of equations using a central difference method if P>0 and
+3. Now set up the system of equations using a central difference method if P>0 and
 $\omega=0$
 
     a. set up the system of equations for 6 segments as a function of q and P
@@ -83,7 +83,7 @@ $\omega=0$
     d. solve a-c for q=1,10,20,30,50 and plot the numerical results of q vs $\delta x$ for
     P=0, 100, 200, 300 (4 lines, labeled as P=0,P=100,...)
 
-5. Now set up an eigenvalue problem to solve for the natural frequencies of the simply
+4. Now set up an eigenvalue problem to solve for the natural frequencies of the simply
 supported beam if P=0 and q=0. 
 
     a. set up the system of equations for 6 segments
@@ -96,7 +96,7 @@ supported beam if P=0 and q=0.
 
     e. Plot the shape of the beam for the first 3 natural frequencies
 
-6. (Bonus 5pt) Create a function to return the system of equations for the eigenvalue
+5. (Bonus 5pt) Create a function to return the system of equations for the eigenvalue
 problem as a function of P, if P>0. Then, plot the lowest natural frequency vs the applied
 load P. 
 
diff --git a/final_project/final_project.pdf b/final_project/final_project.pdf
index 3df5f1f..7c6fa9e 100644
Binary files a/final_project/final_project.pdf and b/final_project/final_project.pdf differ
diff --git a/lecture_24/lecture_24.ipynb b/lecture_24/lecture_24.ipynb
index 9ecde82..66b8be6 100644
--- a/lecture_24/lecture_24.ipynb
+++ b/lecture_24/lecture_24.ipynb
@@ -23,7 +23,9 @@
     "\n",
     "- On the final project, to get the GitHub bonus, do you have to solve the issue? Or do the points go to the one who opens the issue?\n",
     "\n",
-    "- can we go over the final project"
+    "- can we go over the final project\n",
+    "\n",
+    "Tues - more background and help"
    ]
   },
   {
diff --git a/lecture_24/octave-workspace b/lecture_24/octave-workspace
new file mode 100644
index 0000000..8c437bb
Binary files /dev/null and b/lecture_24/octave-workspace differ
diff --git a/lecture_25/final_notes.pdf b/lecture_25/final_notes.pdf
new file mode 100755
index 0000000..ab5c0df
Binary files /dev/null and b/lecture_25/final_notes.pdf differ
diff --git a/lecture_25/lecture_25.pdf b/lecture_25/lecture_25.pdf
new file mode 100755
index 0000000..8141f37
Binary files /dev/null and b/lecture_25/lecture_25.pdf differ
diff --git a/myode.m b/myode.m
new file mode 100644
index 0000000..cedefb6
--- /dev/null
+++ b/myode.m
@@ -0,0 +1,7 @@
+function dudt = myode(t,u)
+    x=u(1);
+    y=u(2);
+    dudt=zeros(2,1);
+    dudt(1)=0.1*x+0.3*x*y;
+    dudt(2)=0.1*x*y-2*y;
+end
\ No newline at end of file