bifurcation diagram

[Giuseppe Ventura]

Hi

Can I create with Scilab a bifurcation diagram?

I have this system of differential equations (in dt):

 

Glucose'=k0*Glucose-k1*Glucose*ATP/(1+(ATP/ki)^n)

ATP'=-2*k1*Glucose*ATP/(1+(ATP/ki)^n)+((k31*k32*TGP*NAD*(atot-ATP)-k33*k34*P
G*ATP*(NADH))/(k33*(NADH)+k32*(atot-ATP)))+(k4*PG*(atot-ATP))-(k7*ATP)

F16P'=k1*Glucose*ATP/(1+(ATP/ki)^n)-k2*F16P

TGP'=2*k2*F16P-((k31*k32*TGP*NAD*(atot-ATP)-k33*k34*PG*ATP*(NADH))/(k33*(NAD
H)+k32*(atot-ATP)))-(k8*TGP*(NADH))

NAD'=(k8*TGP*(NADH))+((k31*k32*TGP*NAD*(atot-ATP)-k33*k34*PG*ATP*(NADH))/(k3
3*(NADH)+k32*(atot-ATP)))+(k6*Acetaldehyde*(NADH))

PG'=((k31*k32*TGP*NAD*(atot-ATP)-k33*k34*PG*ATP*(ntot-NAD))/(k33*(ntot-NAD)+
k32*(atot-ATP)))-(k4*PG*(atot-ATP))

Pyruvate'=k4*PG*(atot-ATP)-(k5*Pyruvate)

Acetaldehyde'=(k5*Pyruvate)-(k6*Acetaldehyde*(NADH))-(k10*(Acetaldehyde-extr
acellularacetaldehyde))

extracellularacetaldehyde'=0.1*(k10*(Acetaldehyde-extracellularacetaldehyde)
)-k9*extracellularacetaldehyde

NADH'=((k31*k32*TGP*NAD*(atot-ATP)-k33*k34*PG*ATP*(NADH))/(k33*(NADH)+k32*(a
tot-ATP)))-(k6*Acetaldehyde*(NADH))-(k8*TGP*(NADH))

 

And I want a diagram Nadh vs k9, Is possible?

I also would like to know what the Hopf point.

 

There is an example?

thanks

Regards

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20101129/7ff60036/attachment.htm>