
/* tcomp34: plot of diffusion of molecule in rod */

timinc = 1e-2;
ploti = timinc;
endexp =.1;
implicit = 1;
crit=1e-8;
relax = .3;

plmax=1;
plmin=0;

drm = 1e30;
vcl = 0;
molconv = 1e3;

/* D = 2e-5; */			/* cm2/sec */
D = 2e-6;			/* cm2/sec */
dri = molconv / D;

segdia = 10;			/* um */
seglen=10;
conn 1 to 2 cable dia segdia length seglen cm 1e-2*segdia*1e-5/4;
plot 1,2;

stimdur = endexp;
stim node 1 cclamp 1e-10 start 0 dur stimdur;

run;

