Bibliography
63
for i=1:N
G4(i) = F4(i);
Flux12(i)=-D*dV1(i)*G4(i);
end
Y4 = diff(G4)/dx; Flux11=-D*Y4;
Flux1=(Flux11+Flux12);
%figure(2); plot(x,dV1);
figure(3); plot(x,G4);
%figure(4); plot(x1,Y4);
save F1.dat F1 -ascii;
save F2.dat F3 -ascii;
save F3.dat F4 -ascii;
save F4.dat F4 -ascii;
save F5.dat F5 -ascii;
save x.dat x -ascii;
%============================================================
%Flux2 IN STATE 2
%============================================================
H4=zeros(2*N,1);
for i=N+1:2*N
H4(i)=F4(i);
end
Z4 = diff(H4)/dx;
dV2=0.;
Flux21=-D*Z4;
%Flux22=-D*(dV2*H6’); Flux22=0.
Flux2 = Flux21;
%============================================================
Totalflow = Flux1 + Flux2;
TotalflowMean = sum(Totalflow)*dx;
disp(TotalflowMean) %Flux mean
%============================================================
%figure(5); plot(x,H4)
%figure(6); plot(x1,Z4)
figure(7); plot(x1,Totalflow)
Bibliography
1. Mogilner, A., Elston, T., Wang, H., Oster, G.: Molecular motors: examples. In: Fall, C.,
Marland, E., Tyson, J., Wagner, J. (eds.) Joel Keizer’s Computational Cell Biology, Chapter 12. Modeling Chemical Reactions vol. 12.3, p. 336 (2002). https://doi.org/10.1016/S03701573(01)00081-3
2. Wang, H., Peskin, C.S., Elston, T.C.: A robust numerical algorithm for studying biomolecular
transport processes. J. Theor. Biol. 221, 491–511 (2003)