This example shows you how to construct a simple model with two species (A and B) and a reaction. The reaction is
A -> B, which follows the mass action kinetics with the forward rate parameter
k. Hence the rate of change is .
Create a SimBiology model named
m1 = sbiomodel('simpleModel');
Add a reaction that involves two species
A is converted to
r1 = addreaction(m1,'A -> B');
SimBiology automatically add species
B to the model.
SimBiology Species Array Index: Compartment: Name: InitialAmount: InitialAmountUnits: 1 unnamed A 0 2 unnamed B 0
Set the initial amount of the first species (
A) to 10.
m1.species(1).InitialAmount = 10;
Define the kinetic law of the reaction to follow the mass action kinetics. You can achieve this by adding a kinetic law object to the reaction
kineticLaw = addkineticlaw(r1,'MassAction');
Add a rate constant parameter to the mass action kinetic law. You must set the
ParameterVariableNames property of the kinetic law object to the name of the parameter
'k1' so that the reaction rate can be determined.
p1 = addparameter(kineticLaw,'k',0.5); kineticLaw.ParameterVariableNames = 'k';
Simulate the model.
sd = sbiosimulate(m1);
Plot the simulation results.