t is twice beta times xi*li</span>
alpha = 1;   <span class="comment">% conductance per segment is alpha times xi/li</span>
G0 = 1;      <span class="comment">% source conductance</span>
C0 = 10;     <span class="comment">% load capacitor</span>
wmax = 1;    <span class="comment">% upper bound on x</span>

<span class="comment">%</span>
<span class="comment">% Node positions</span>
<span class="comment">%</span>

xpos = [ 0   1   6   8  -4  -1 ;
         0  -1   4  -2   1   4 ] ;
X11 = repmat(xpos(1,:),n,1)