File Exchange

image thumbnail

M-File Application OF Model Pcemaker VI1 of Cardiac Cells With Runge Kutta 4

version 1.0 (1.47 KB) by

it is M-file of runge kutta 4 of model pacemaker VI1 for all the cardiac cells parametrized

0 Downloads

Updated

No License

It is M-file Runge Kutta 4 for model pacemaker VI1 of all the cardiac cells

Comments and Ratings (3)

Urs (us) Schwarz

simply said: mental SCHMUTZ...
to be removed at your earliest convenience
us

John D'Errico

Again, more dreck. The author could at least have spelled MATLAB correctly. Its not "math lab". While this is actually an m-file, its essentially uncommented, and so useless to anyone but the author.

Dimitri Shvorob

First of all, it doesn't run: the m-file name and function name differ. What this in-essence-script does is:

for n=1:1:1000
A1(1,n)=(a/C)*(V(n)^3)+(b/C)*(V(n)^2)+(do/C)*V(n)+d1/C -(1/C)*I1(n)+(A*sin(2*(3.14)*(n*h)*F+K))/C;
B1(1,n)=(1/L)*V(n)-(R/L)*I1(n)+E/L;
A2(1,n)=(a/C)*(V(n)+(h/2)*A1(n)).^3+(b/C)*(V(n)+(h/2)*A1(n)).^2+(do/C)*(V(n)+(h/2)*A1(n))+d1/C-(1/C)*(I1(n)+(h/2)*B1(n))+(A*sin((2*3.14)*F*(n*h)+K))/C;
B2(1,n)=(1/L)*(V(n)+(h/2)*A1(n))-(R/L)*(I1(n)+(h/2)*B1(n))+E/L;
A3(1,n)=(a/C)*(V(n)+(h/2)*A2(n)).^3 +(b/C)*(V(n)+(h/2)*A2(n)).^2 +(do/C)*(V(n)+(h/2)*A2(n)) +d1/C -(1/C)*(I1(n)+(h/2)*B2(n))+(A*sin((2*3.14)*F*(n*h)+K))/C;
B3(1,n)=(1/L)*(V(n)+(h/2)*A2(n))-(R/L)*(I1(n)+(h/2)*B2(n))+E/L;
A4(1,n)=(a/C)*(V(n)+h*A3(n)).^3+(b/C)*(V(n)+h*A3(n)).^2+(do/C)*(V(n)+h*A3(n))+d1/C-(1/C)*(I1(n)+h*B3(n))+(A*sin((2*3.14)*F*(n*h)+K))/C;
B4(1,n)=(1/L)*(V(n)+h*A3(n))-(R/L)*(I1(n)+h*B3(n))+E/L;
V(1,n+1)=V(n)+(h/6)*(A1(n)+A2(n)+2*A3(n)+A4(n));
I1(1,n+1)=I1(n)+(h/6)*(B1(n)+B2(n)+2*B3(n)+B4(n));
end

MATLAB Release
MATLAB 7.0.1 (R14SP1)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video