Demo Dose Response Command
This scripts demonstrates how to use the doseResponse command. This command is used to find the steady state concenetrations for a model given a particular input u. Mathematically this means solving for x that is consitent with the initial conditions such that:
We solve this equation using Newton's method which allows the method to find both stable and unstable fixed points.
Contents
Load the model
First we load the model using the loadModel command. In this case it is a model of a a distributive MAPK phosphorylation/dephosphophorylation loop.
model = loadModel('OneStep');
Graph the model structure
You can check how the process has proceeded by using the model graphing feature. This shows all species, and their connectivity in the reaction network.
graphReactions(model);
