From a0711c3670c79d40295d8f02d1388cca9d5b4dc8 Mon Sep 17 00:00:00 2001
From: Ju Eun Kang <133086206+jueun429@users.noreply.github.com>
Date: Mon, 14 Jul 2025 15:08:31 +0900
Subject: [PATCH] Make simulations deterministic

---
 GynCycle_Original code/CreateFollicles.m           | 4 +++-
 GynCycle_Original code/HumanSimulationFollGrowth.m | 4 +++-
 GynCycle_Original code/poissonproc.m               | 4 +++-
 3 files changed, 9 insertions(+), 3 deletions(-)

diff --git a/GynCycle_Original code/CreateFollicles.m b/GynCycle_Original code/CreateFollicles.m
index c925d06..f81afa6 100644
--- a/GynCycle_Original code/CreateFollicles.m	
+++ b/GynCycle_Original code/CreateFollicles.m	
@@ -1,4 +1,6 @@
-function [FSHVec, StartVec] = CreateFollicles(para,paraPoi,tb,te)
+function [FSHVec, StartVec] = CreateFollicles(para,paraPoi,tb,te)
+% use a fixed random seed for reproducibility
+rng(0,"twister");
 
 %create normal distributed fsh sensitivities for each foll
 fileID2 = fopen('FSH.txt','w+');
diff --git a/GynCycle_Original code/HumanSimulationFollGrowth.m b/GynCycle_Original code/HumanSimulationFollGrowth.m
index 02ab5b9..92bab32 100644
--- a/GynCycle_Original code/HumanSimulationFollGrowth.m	
+++ b/GynCycle_Original code/HumanSimulationFollGrowth.m	
@@ -1,7 +1,9 @@
 %set specifications for simulations a call calculation and solving of the ODE
 %needed:'Model28_Parameter.dat', 'InitialValues.txt'
 
-function HumanSimulationFollGrowth
+function HumanSimulationFollGrowth
+% initialise RNG for deterministic behaviour
+rng(0,"twister");
 %
 %-----------------------------------------------------------------------
 %
diff --git a/GynCycle_Original code/poissonproc.m b/GynCycle_Original code/poissonproc.m
index 91003cf..f170e0a 100644
--- a/GynCycle_Original code/poissonproc.m	
+++ b/GynCycle_Original code/poissonproc.m	
@@ -1,4 +1,6 @@
-function S =poissonproc(lambda,tspan)
+function S =poissonproc(lambda,tspan)
+% deterministic Poisson process based on a fixed seed
+rng(0,"twister");
 
 tb=tspan(1);    %start time
 te=tspan(2);  %end time