Add a seed in the header of output file to ensure reproductibility

hybrid.pl

@@ -109,8 +109,10 @@ hybrid_static_simulation(ODEFilename, StochFilename, OutFilename, Volume, Time)
          nl(Stream)
          ),
          _),
+  format(Stream, "% Hybrid version of ~w and ~w ~n", ODEFilename, StochFilename),
   write(Stream, 'parameter(tau = 0).'), nl(Stream), 
   Seed is 0,
+  format(Stream, "seed(~d).~n", [Seed]),
   write(Stream, 'parameter(ran = '),
   write(Stream, Seed),
   write(Stream, ').'), nl(Stream),
@@ -1170,8 +1172,10 @@ hybrid_dynamic_simulation(InputFile,OutFilename,Volume,Time,PropTresh,PartTresh)
   sort(Reactants1,SortedReactants), % Used for later hybrid-reactants finding and constraints lists 
   sort(Species1,SortedSpecies),
 
+  format(Stream, "% Hybrid version of ~w~n", InputFile),
   write(Stream, 'parameter(tau = 0).'), nl(Stream), 
   Seed is 0,
+  format(Stream, "seed(~d).~n", [Seed]),
   write(Stream, 'parameter(ran = '),
   write(Stream, Seed),
   write(Stream, ').'), nl(Stream),