<?xml version="1.0" encoding="UTF-8"?>
<!-- Created by Gepasi 3.30 on September 03, 2002, 23:14 -->
<sbml xmlns="http://www.sbml.org/sbml/level1" level="1" version="1">
	<model name="The_Brusselator">
		<notes>
			<body xmlns="http://www.w3.org/1999/xhtml">
				<p>The famous Brusselator model in which a chemical species X catalyses</p>
				<p>its own formation. In a certain region of parameter space the Brusselator</p>
				<p>can function in a limit cycle regime around an unstable steady state.</p>
				<p></p>
				<p>Run this simulation and then check the report file to see how Gepasi can find an unstable steady state, and calculate the eigenvalues of the Jacobian matrix.</p>
			</body>
		</notes>
		<listOfCompartments>
			<compartment name="compartment" volume="1"/>
		</listOfCompartments>
		<listOfSpecies>
			<specie name="A" compartment="compartment" initialAmount="0.5" boundaryCondition="true"/>
			<specie name="X" compartment="compartment" initialAmount="3" boundaryCondition="false"/>
			<specie name="Y" compartment="compartment" initialAmount="3" boundaryCondition="false"/>
			<specie name="B" compartment="compartment" initialAmount="3" boundaryCondition="true"/>
			<specie name="D" compartment="compartment" initialAmount="0" boundaryCondition="true"/>
			<specie name="E" compartment="compartment" initialAmount="0" boundaryCondition="true"/>
		</listOfSpecies>
		<listOfReactions>
			<reaction name="R1" reversible="false">
				<listOfReactants>
					<specieReference specie="A" stoichiometry="1"/>
				</listOfReactants>
				<listOfProducts>
					<specieReference specie="X" stoichiometry="1"/>
				</listOfProducts>
				<kineticLaw formula="k_1*A">
					<listOfParameters>
						<parameter name="k_1" value="1"/>
					</listOfParameters>
				</kineticLaw>
			</reaction>
			<reaction name="R2" reversible="false">
				<listOfReactants>
					<specieReference specie="X" stoichiometry="2"/>
					<specieReference specie="Y" stoichiometry="1"/>
				</listOfReactants>
				<listOfProducts>
					<specieReference specie="X" stoichiometry="3"/>
				</listOfProducts>
				<kineticLaw formula="k_2*X*X*Y">
					<listOfParameters>
						<parameter name="k_2" value="1"/>
					</listOfParameters>
				</kineticLaw>
			</reaction>
			<reaction name="R3" reversible="false">
				<listOfReactants>
					<specieReference specie="X" stoichiometry="1"/>
					<specieReference specie="B" stoichiometry="1"/>
				</listOfReactants>
				<listOfProducts>
					<specieReference specie="Y" stoichiometry="1"/>
					<specieReference specie="D" stoichiometry="1"/>
				</listOfProducts>
				<kineticLaw formula="k_3*X*B">
					<listOfParameters>
						<parameter name="k_3" value="1"/>
					</listOfParameters>
				</kineticLaw>
			</reaction>
			<reaction name="R4" reversible="false">
				<listOfReactants>
					<specieReference specie="X" stoichiometry="1"/>
				</listOfReactants>
				<listOfProducts>
					<specieReference specie="E" stoichiometry="1"/>
				</listOfProducts>
				<kineticLaw formula="k_4*X">
					<listOfParameters>
						<parameter name="k_4" value="1"/>
					</listOfParameters>
				</kineticLaw>
			</reaction>
		</listOfReactions>
	</model>
</sbml>
