## 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