diff --git a/NregNsyn.hoc b/NregNsyn.hoc
index f9c93d8..c1bea91 100644
--- a/NregNsyn.hoc
+++ b/NregNsyn.hoc
@@ -9,6 +9,7 @@ printf("irregular & asynchronous input. frequency = %d \n", freq)
 
 objref r
 r = new Random()
+r.Random123()
 
 
 for i=0, totSyn-1 {  //totSyn is the total # of synapses
diff --git a/createNewSyn.hoc b/createNewSyn.hoc
index 20a964b..85cd18f 100644
--- a/createNewSyn.hoc
+++ b/createNewSyn.hoc
@@ -60,6 +60,7 @@ for (i=0; i <= totSyn-1; i +=1) {
     dend[comp].sec {
 	dist = distance(tmp)
         Ens[i] = new NetStim(0.5)
+        Ens[i].noiseFromRandom123(i, 1, 0)
         Ens[i].interval=200
         Ens[i].number = 1000000000000000000000000000000
         Ens[i].start=0
diff --git a/createNewSyn2.hoc b/createNewSyn2.hoc
index 3540a15..1423629 100644
--- a/createNewSyn2.hoc
+++ b/createNewSyn2.hoc
@@ -27,6 +27,7 @@ totSyn = 1000
 
 objref Ens
 Ens = new NetStim(0.5)
+Ens.noiseFromRandom123(0, 2, 0)
 Ens.interval=200
 Ens.number = 1000000000000000000000000000000
 Ens.start=0
diff --git a/readme.txt b/readme.txt
index c69cf42..5260b82 100644
--- a/readme.txt
+++ b/readme.txt
@@ -32,3 +32,8 @@ will open the simulation window.
 
 Questions on how to use this model
 should be directed to xsli2@yahoo.com
+
+Updates
+
+20240226 Use noiseFromRandom123 for result identity between 8.2 and 9.0
+2024-11: Updated hoc files to use Random123
diff --git a/regNsyn.hoc b/regNsyn.hoc
index a19f42a..091b46b 100644
--- a/regNsyn.hoc
+++ b/regNsyn.hoc
@@ -9,6 +9,7 @@ printf("regular & asynchronous input. frequency = %d \n", freq)
 
 objref r
 r = new Random()
+r.Random123()
 
 for i=0, totSyn-1 {  //totSyn is the total # of synapses
 	Ens[i].start=0 //Ens: the reference to the NetStim