/* tcomp21:  test of rod */
/* same as tcomp15 but uses "rod.m" */



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

pois = 0;

 include "rod.m";	/* */

n = 100;

n = mrod(0,0,n+1);

conn n to 6  cable length 10 dia .1;

conn 6 to 7 cable length 10 dia .5;
conn 7 to 8 cable length 10 dia 2;
at 8 sphere dia 30; 

conn 8 to 9 cable length 50 dia 5;
						/* */

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

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

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

plot V[8],V[7],V[101];		/* */

 run; 

