diff --git a/samples/tutorial_1.6/MakeDef.py b/samples/tutorial_1.6/MakeDef.py
index 467a9e087..a864c777d 100644
--- a/samples/tutorial_1.6/MakeDef.py
+++ b/samples/tutorial_1.6/MakeDef.py
@@ -15,7 +15,7 @@
 #[s] set input.
 list_input = ['Ns','exct']
 dict_input = read.func_input(list_input)
-max_site = 2
+max_site = int(dict_input["Ns"])
 exct     = 2
 All_N    = max_site
 #[e] set input.
@@ -61,7 +61,7 @@
         InterAll[I][1][I][0][J][1][J][0] += complex(0,0.25*para[cnt])
     if intT1[cnt] == 'x' and intT2[cnt] == 'z':
         I = siteI[cnt]
-        J = siteI[cnt]
+        J = siteJ[cnt]
         # xz
         InterAll[I][0][I][1][J][0][J][0] +=  0.25*para[cnt]
         InterAll[I][0][I][1][J][1][J][1] += -0.25*para[cnt]
@@ -74,7 +74,7 @@
         InterAll[I][1][I][1][J][1][J][0] += -0.25*para[cnt]
     if intT1[cnt] == 'y' and intT2[cnt] == 'z':
         I = siteI[cnt]
-        J = siteI[cnt]
+        J = siteJ[cnt]
         # yz
         InterAll[I][0][I][1][J][0][J][0] += complex(0,-0.25*para[cnt])
         InterAll[I][0][I][1][J][1][J][1] += complex(0,0.25*para[cnt])