/* tcomp22:  test of 2 rods */
/* same as tcomp21 but uses "rod.m" */



plmax = .0;
plmin = -.06;
endexp = 1.0;
ploti = endexp/200;
drm = 10000;
relax= 0.57;

rodrm = 5000;

 include "rod.m";	/* */

p1 = mrod(0,0,100);

conn [p1] to [120] cable length 10 dia .1;
conn [120] to [130] cable length 10 dia .5;

p2 = mrod(2,0,200);

conn [p2] to [220] cable length 10 dia .1;
conn [220] to [130] cable length 10 dia .5;

conn [130] to [140] cable length 10 dia 2;
at [140] sphere dia 30; 

conn [140] to [150] cable length 50 dia 5;
						/* */

stim rod 100 inten 100e3 start .009 dur .001;
stim rod 200 inten 100e3 start .009 dur .001;

/*stim node 1 vclamp -.06 start .005 dur 1;  */

/*plot I[1] max 40e-12 min 0; */

plot V[140],V[p2],V[p1],V[200],V[100];				 /* */

run;

