%% Demo Dose Response Command
%
% <<../../../images/LogoSmall.jpg>>
%
% 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:
%
% $$\dot{x} = A_1 x + A_2 x \otimes  x + B_1 u + B_2 x \otimes u + K = 0$$
%
% We solve this equation using Newton's method which allows the method to
% find both stable and unstable fixed points.

%% 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);