MATH 592: TOPICS IN APPLIED MATHEMATICS II
Mathematical Physiology 2004
Documentation PDF documentation for xppaut
Hodgkin-Huxley Model xppaut code for Hodgkin -Huxley model of squid axon electrical activity.
FitzHugh-Nagumo Model xppaut code for a FitzHugh Nagumo model with nonlinearity f(u)=a*u*(u^3-u)/3. Set for AUTO calculations using applied current i as parameter.
Two Pool Calcium Model xppaut code for a 2-pool model of intracellular calcium. In the model u=cytosolic calcium and v=calcium concentration in an internal store with calcium inactivated calcium release (CICR). The calcium in the IP3 sensitive store is assumed constant. Initial conditions are set for an AUTO calculation.
Chay Cook beta-cell model xppaut code for a model of beta cell electrical activity. For different parameter values the model exhibits a variety of qualitatively different bursting behaviors (See Bull. Math. Biol. Vol 57, 00. 413-439, 1995 for a detailed account). The fast subsystem for the model with s=sinf(v,c) and c a bifurcation parameter can be downloaded by clicking here.
Polynomial model for Bursting For different parameter values model exhibits a variety of bursting patterns. Initial parameters are set for bursting patterns analogous to those in pancreatic beta-cells. (see: SIAM J. Appl. Math. Vol 54, pp. 814-832, 1994). The fast subsystem for the model can be obtained by clicking here.
Documentation Latex documentation for xtc - long version with examples.
Summary Latex documentation for xtc - summary.
Facilitated Diffusion of Oxygen XTC code for a dimensionless model of facilitated diffusion of oxygen in muscle fibers. Simple model uses S=oxygen, E=myoglobin and C=oxymyoglobin
Enzyme Capsule Model XTC code for numerical solution of enzyme capsule model in question 1 of Assignment 2 above.
Wavefront in FitzHugh Nagumo model #1 Numerically demonstrates wavefronts for a FHN model of voltage v(x,t) where x is on a finite interval. The kinetic portion of the model is set to the oscillatory mode and the stimulus is introduced by imposing v_x = ia2 at the right end.
Wavefront in FitzHugh Nagumo model #2 Numerically demonstrates travelling pulses for a FHN model of voltage v(x,t) where x is on a finite interval. The kinetic portion of the model is set to the stable steady state mode and the stimulus is introduced into the v eqn as follows: v_t = d v_xx + f(v) - w + I(x,t) where I is a pulse function equal to ia for x in (x0,x0+dx), t in (t0,t0+dt).
Wavefront in FitzHugh Nagumo model #3 Numerically demonstrates travelling pulses for a FHN model of voltage v(x,t) where x is on a finite interval. The kinetic portion of the model is set to the stable steady state mode and the stimulus is introduced as a short pulse at x=1 through the boundary condition v_x=I(t) where I is approximately equal to ia=1 for a duration t in (t0,t0+dt).