diff --git a/matrices/laplacian/julia-laplaciano.jl b/matrices/laplacian/julia-laplaciano.jl
index 6ca4fcb..e833cf5 100644
--- a/matrices/laplacian/julia-laplaciano.jl
+++ b/matrices/laplacian/julia-laplaciano.jl
@@ -8,17 +8,17 @@ nx = N
 ny = N
 nz = N
 
-ex = fill(1, nx)
-ey = fill(1, ny)
-ez = fill(1, nz)
+ex = fill(1, nx - 1)
+ey = fill(1, ny - 1)
+ez = fill(1, nz - 1)
 
-Dxx = spdiagm(-1 => ex, 0 => -2 * ex, +1 => ex)
-Dyy = spdiagm(-1 => ey, 0 => -2 * ey, +1 => ey)
-Dzz = spdiagm(-1 => ez, 0 => -2 * ez, +1 => ez)
+Dxx = spdiagm(-1 => ex, 0 => -2 * [ex; 1], +1 => ex)
+Dyy = spdiagm(-1 => ey, 0 => -2 * [ey; 1], +1 => ey)
+Dzz = spdiagm(-1 => ez, 0 => -2 * [ez; 1], +1 => ez)
 
-Ix = spdiagm(0 => [ex; 1])
-Iy = spdiagm(0 => [ey; 1])
-Iz = spdiagm(0 => [ez; 1])
+Ix = spdiagm(0 => ex)
+Iy = spdiagm(0 => ey)
+Iz = spdiagm(0 => ez)
 
 L = kron(Dxx, Iy, Iz) + kron(Ix, Dyy, Iz) + kron(Ix, Iy, Dzz)
 
diff --git a/matrices/laplacian/laplacian-discretization-3d.mat b/matrices/laplacian/laplacian-discretization-3d.mat
index d33a882..eac550e 100644
Binary files a/matrices/laplacian/laplacian-discretization-3d.mat and b/matrices/laplacian/laplacian-discretization-3d.mat differ