Schlogl’s model¶
- class pints.toy.stochastic.SchloglModel(initial_molecule_count=20)[source]¶
Schlogl’s system of chemical reactions has a single type of molecules and starts with an initial count \(A(0)\). The evolution of the molecule count is defined through the rates \(k_1\), \(k_2\), \(k_3\) and \(k_4\) and the following equations:
- ..math::
2A xrightarrow{k_1} 3A 3A xrightarrow{k_2} 2A 0 xrightarrow{k_3} A A xrightarrow{k_4} 0
Extends
pints.ForwardModel
,pints.toy.ToyModel
.- Parameters:
initial_molecule_count – The initial molecule count \(A(0)\).
- interpolate_mol_counts(time, mol_count, output_times)¶
Takes raw times and inputs and mol counts and outputs interpolated values at output_times
- n_outputs()¶
Returns the number of outputs this model has. The default is 1.
- n_parameters()¶
- simulate(parameters, times)¶
- simulate_raw(rates, max_time)¶
Returns raw times, mol counts when reactions occur.