diff --git a/lu_tridiag.m b/lu_tridiag.m
index 85367f9..081e926 100644
--- a/lu_tridiag.m
+++ b/lu_tridiag.m
@@ -1,4 +1,3 @@
-
 function [ud,uo,lo]=lu_tridiag(e,f,g);
   % lu_tridiag calculates the components for LU-decomposition of a tridiagonal matrix
   % given its off-diagonal vectors, e (a) and g (b)
@@ -12,15 +11,15 @@
 ud=zeros(m,n);
 uo=zeros(m,n);
 lo=zeros(m,n);
-%ud(1)=f(1);
-%lo(1)=0;
-%uo=e;
+ud(1)=f(1);
+lo(1)=0;
+uo=e;
 
 
-%for i=2:n
- %   lo(i)=g(i)/(ud(i-1));
-  %  ud(i)=f(i)-(lo(i)*e(i-1));
-%end 
+for i=2:n
+    lo(i)=g(i)/(ud(i-1));
+    ud(i)=f(i)-(lo(i)*e(i-1));
+end