Bueno-Orovio - Cherry - Fenton (BCF) model. Temporal dynamics

    du/dt = -(Jfi + Jso + Jsi ) ,
    dv/dt = [1 - H(u - Vv )](v - v)/τv - - H(u - Vv )v/τv ,
    dw/dt = [1 - H(u - Vw )](w - w)/τw- - H(u - Vw )w/τw ,        
    ds/dt = ((1 + tanh[ks (u - Us )])/2 - s)/τs ,

    Jfi = -vH(u - Vv)(u - Vv)(Uu - u)/τfi ,
    Jso = (u - Uo)[1 - H(u - Vw )]/τo + H(u - Vw )/τso ,
    Jsi = -H(u - Vw)ws/τsi ,

    τv- = [1 - H(u - Vv v1 + H(u - Vv v2,
    τs = [1 - H(u - Vw s1 + H(u - Vw s2 ,
    τo = [1 - H(u - Vo )]τo1 + H(u - Vo o2 ,
    τw- = τw1 + (τw2 - τw1 ) (1 + tanh[kw(u - Uw)])/2,
    τso = τso1 + (τso2 - τso1-) (1 + tanh[kso (u - Uso )])/2,
    v = 1 - H(u - Vv ) ,
    w = [1 - H(u - Vo )](1 - u/τw∞ ) + H(u - Vo )w*,

    V mV = 85.7u - 84 .

Uo = 0
Uu = 1.58
Vv = 0.3
Vw = 0.015
Vv- = 0.015
Vo = 0.006
τv1 = 60
τv2 = 1150
τv = 1.4506        
τw1 = 70
τw2 = 20
kw = 65
Uw = 0.03
τw = 280
τfi = 0.11
τo1 = 6
τo2 = 6
τso1 = 43
τso2 = 0.2
kso = 2
Uso = 0.65
τs1 = 2.7342
τs2 = 3
ks = 2.0994
Us = 0.9087
τsi = 2.8723
τw∞ = 0.07
w* = 0.94
Time is mesured in msec. Below you can explore the model for different parameter values.

Uo Vo Wo So Ymax dtms it
u, v, w, s variables are black, blue, red and green curves. The script makes 1000 it time steps dt.
Heart rhythms     updated 8 Dec 2011