diff --git a/Fig13A.hoc b/Fig13A.hoc
index a50bc23..d9f6068 100644
--- a/Fig13A.hoc
+++ b/Fig13A.hoc
@@ -175,6 +175,7 @@ assign_synapses(0.01)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig13B.hoc b/Fig13B.hoc
index 9c1fa3f..036f303 100644
--- a/Fig13B.hoc
+++ b/Fig13B.hoc
@@ -176,6 +176,7 @@ assign_synapses(0.005)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig5A.hoc b/Fig5A.hoc
index dfb8b83..9fc723c 100644
--- a/Fig5A.hoc
+++ b/Fig5A.hoc
@@ -175,6 +175,7 @@ assign_synapses(1)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig5B.hoc b/Fig5B.hoc
index 171de8d..aae5781 100644
--- a/Fig5B.hoc
+++ b/Fig5B.hoc
@@ -177,6 +177,7 @@ assign_synapses(0.5)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig6A.hoc b/Fig6A.hoc
index e760860..062186a 100644
--- a/Fig6A.hoc
+++ b/Fig6A.hoc
@@ -209,6 +209,7 @@ assign_synapses(0.5)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig6B.hoc b/Fig6B.hoc
index 5fb32c8..4ebbe4c 100644
--- a/Fig6B.hoc
+++ b/Fig6B.hoc
@@ -209,6 +209,7 @@ assign_synapses(0.333)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig6C.hoc b/Fig6C.hoc
index 1dfc0ba..d39474d 100644
--- a/Fig6C.hoc
+++ b/Fig6C.hoc
@@ -214,6 +214,7 @@ assign_synapses(0.25)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/Fig6D.hoc b/Fig6D.hoc
index 784da03..62fc4b7 100644
--- a/Fig6D.hoc
+++ b/Fig6D.hoc
@@ -209,6 +209,7 @@ assign_synapses(0.2)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"
diff --git a/README b/README
index 93fc683..618ce35 100644
--- a/README
+++ b/README
@@ -122,3 +122,4 @@ ModelDB-Note: This model updated for NEURON version 4.3.1 and higher.
 10-2007: AD, better synaptic mechanisms using counters
 20120106 changed mod files for better solve methods (avoid euler)
 as per web page http://www.neuron.yale.edu/phpbb/viewtopic.php?f=28&t=592
+2024-11: Updated hoc files to use Random123