diff --git a/contents/gaussian_elimination/code/julia/gaussian_elimination.jl b/contents/gaussian_elimination/code/julia/gaussian_elimination.jl
index 22a87cdb5..24f47d91f 100644
--- a/contents/gaussian_elimination/code/julia/gaussian_elimination.jl
+++ b/contents/gaussian_elimination/code/julia/gaussian_elimination.jl
@@ -52,10 +52,7 @@ function back_substitution(A::Array{Float64,2})
     # Creating the solution Vector
     soln = Vector{Float64}(rows)
 
-    # initialize the final element
-    soln[rows] = A[rows, cols] / A[rows, cols-1]
-
-    for i = (rows - 1):-1:1
+    for i = rows:-1:1
         sum = 0.0
         for j = rows:-1:i
             sum += soln[j]*A[i,j]
diff --git a/contents/gaussian_elimination/gaussian_elimination.md b/contents/gaussian_elimination/gaussian_elimination.md
index f7e31bfdb..34c4d3ca1 100644
--- a/contents/gaussian_elimination/gaussian_elimination.md
+++ b/contents/gaussian_elimination/gaussian_elimination.md
@@ -385,13 +385,13 @@ Here it is in code:
 
 {% method %}
 {% sample lang="jl" %}
-[import:70-96, lang:"julia"](code/julia/gaussian_elimination.jl)
+[import:67-93, lang:"julia"](code/julia/gaussian_elimination.jl)
 {% sample lang="c" %}
 This code does not exist yet in C, so here's Julia code (sorry for the inconvenience)
-[import:70-96, lang:"julia"](code/julia/gaussian_elimination.jl)
+[import:67-93, lang:"julia"](code/julia/gaussian_elimination.jl)
 {% sample lang="rs" %}
 This code does not exist yet in rust, so here's Julia code (sorry for the inconvenience)
-[import:70-96, lang:"julia"](code/julia/gaussian_elimination.jl)
+[import:67-93, lang:"julia"](code/julia/gaussian_elimination.jl)
 {% sample lang="hs" %}
 [import:38-46, lang:"haskell"](code/haskell/gaussianElimination.hs)
 {% endmethod %}