// Created by libAntimony v2.8.0
model *Rattanakul2003_BoneFormationModel()
// Compartments and Species:
compartment Compartment_;
species $x in Compartment_, $y in Compartment_, $z in Compartment_;
// Rate Rules:
x' = a1/(k1 + y) - b1*x;
y' = epsilon*(((a2 + a3*x)*y*z)/(k2 + x^2) - b2*y);
z' = epsilon*delta*(a4*x - (b3*z + (a5*x*z)/(k3 + x)));
// Species initializations:
x = 2;
y = 1;
z = 0.15;
// Compartment initializations:
Compartment_ = 1;
Compartment_ has litre;
// Variable initializations:
a1 = 0.05;
k1 = 0.1;
b1 = 0.1;
epsilon = 0.1;
a2 = 0.009;
a3 = 0.675;
k2 = 0.5;
b2 = 0.3;
delta = 0.9;
a4 = 0.01;
b3 = 0.01;
a5 = 0.005;
k3 = 0.025;
// Other declarations:
const Compartment_, a1, k1, b1, epsilon, a2, a3, k2, b2, delta, a4, b3;
const a5, k3;
// Unit definitions:
unit minute = 60 second;
unit time_unit = 60 second;
// Display Names:
time_unit is "time";
x is "PTH";
y is "active osteoclasts";
z is "active osteoblasts";
end