MATH 592: TOPICS IN APPLIED MATHEMATICS II

Mathematical Physiology 2004

Grading and Project Outline

Notes:

Computer Labs:

Figures

XPPAUT Code and documentation

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.

XTC Code and documentation

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).