Prev: Please Somebody Help! Digital hologram
Next: Best OS/platform for Matlab
From: Manjeera on 29 Jul 2010 12:59
Hi I'm a Biology student. I am trying to solve a set of 7 ODE simultaneously. These are the signals in a cancerous cell. I think there is something wrong with the code but I am not able to figure it out. Can someone please help?

% parametric values%
k1 = 1; V2p= 0.005; V2dp= 0.25; k3= 0.005;
V25p = 0.017; V25dp= 0.17; Vweep = 0.01; Vweedp= 1;
kA = 0.01; KKA =0.1; kB =1; KKB =1;
kC = 0.13; KKC =0.01;kD = 0.13; KKD = 1;
kE =0.02 ; KKE = 0.1, kF = 0.1; KKF =1;
kG =0.02; KKG =0.01; kH =0.15; KKH = 0.01;
TotCdk = 100;

% initial conditions %
Cyclin(0) =0; MPF (0)=0 ;preMPF(0)= 0 ; Cdc25P(0)= 0; Wee1P(0) = 0;
IEP(0)= 1; APC(0) = 1;

% conservation of Cdk subunits%
Cdk = TotCdk - MPF - preMPF;

%parameters set to 1%
TotCdc25 =1; TotWee1 =1; TotIE =1; TotAPC =1; PPase= 1;

%rate functions embodying the regulatory signals%
k25 = V25p*(TotCdc25-Cdc25P)+(V25dp*Cdc25P);
kwee = Vweep*Wee1P + Vweedp*(TotWee1-Wee1P);
k2 = V2p*(TotAPC-APC)+V2dp*APC;


Cyclin = dsolve('D(Cyclin) = k1 - (k2*Cyclin)*Cdk', 'Cyclin(0)=0', 't');
MPF = dsolve('D(MPF) = (k3*Cyclin*Cdk)-(k2*-kwee)*MPF+(k25*preMPF)' , 'MPF(0)=0', 't');
preMPF= dsolve('D(preMPF)= ((kwee*MPF)-(k25=k2)*preMPF)','preMPF(0)= 0','t');
Cdc25P = dsolve('D(Cdc25P)=(kA*MPF*(TotCdc25-Cdc25P)/(KKA+(TotCdc25-Cdc25P)))-(kB*PPase*Cdc25P/(KKB+Cdc25P))','Cdc25P(0)=0','t');
Wee1P = dsolve('D(Wee1P)=(kE*MPF*(TotWee1-Wee1P)/(KKE+(TotWee1-Wee1P)))-(kF*PPase*Wee1P/(KKF+Wee1P))','Wee1P(0)=0','t');
IEP = dsolve('D(IEP)=(kG*MPF*(TotIE-IEP)/(KKG+(TotIE-IEP)))-(kH*PPase*IEP/(KKH+IEP))','IEP(0)=1','t');
APC = dsolve('D(APC)=(kC*IEP*(TotAPC-APC)/(KKC+(TotAPC-APC)))-(kD*PPase*APC/(KKD+APC))','APC(0)=1','t');



My objective is to solve these equations and plot it on a single graph.
I would really appreciate it if someone could help me with this.

Thank you

Manjeera
From: Alan Weiss on 29 Jul 2010 13:21
On 7/29/2010 12:59 PM, Manjeera wrote:
> Hi I'm a Biology student. I am trying to solve a set of 7 ODE
> simultaneously. These are the signals in a cancerous cell. I think there
> is something wrong with the code but I am not able to figure it out. Can
> someone please help?
(SNIP)
>
> Manjeera

I am not sure why you are trying to solve these equations symbolically
instead of numerically. Take a look at how to solve differential
equations numerically in MATLAB:
http://www.mathworks.com/access/helpdesk/help/techdoc/math/f1-662913.html

Alan Weiss
MATLAB mathematical toolbox documentation
From: Manjeera on 30 Jul 2010 13:08
Alan Weiss <aweiss(a)mathworks.com> wrote in message <i2sdae$ke5$1(a)fred.mathworks.com>...

> I am not sure why you are trying to solve these equations symbolically
> instead of numerically. Take a look at how to solve differential
> equations numerically in MATLAB:
> http://www.mathworks.com/access/helpdesk/help/techdoc/math/f1-662913.html
>
> Alan Weiss
> MATLAB mathematical toolbox documentation

Thank You Alan,
The main reason i am trying to solve this graphically is to compare results graphically to another method used to solve these equations. The results of another method is just in form of graphs.
I will try solving it numerically and then try plot everything.
From: Alan Weiss on 30 Jul 2010 13:56
On 7/30/2010 1:08 PM, Manjeera wrote:
> Alan Weiss <aweiss(a)mathworks.com> wrote in message
> <i2sdae$ke5$1(a)fred.mathworks.com>...
>
>> I am not sure why you are trying to solve these equations symbolically
>> instead of numerically. Take a look at how to solve differential
>> equations numerically in MATLAB:
>> http://www.mathworks.com/access/helpdesk/help/techdoc/math/f1-662913.html
>>
>> Alan Weiss
>> MATLAB mathematical toolbox documentation
>
> Thank You Alan,
> The main reason i am trying to solve this graphically is to compare
> results graphically to another method used to solve these equations. The
> results of another method is just in form of graphs. I will try solving
> it numerically and then try plot everything.

You might want to look at the plot of a solution to a 3-D ODE towards
the end of this demo:
http://www.mathworks.com/products/optimization/demos.html?file=/products/demos/shipping/optim/goaldemo.html

Alan Weiss
MATLAB mathematical toolbox documentation
 | 
Pages: 1
Prev: Please Somebody Help! Digital hologram
Next: Best OS/platform for Matlab