                        <span class="comment">% number of material points &amp; dimensions</span>
0032 fext = zeros(nodes*nD,1);                                                   <span class="comment">% zero the external force vector</span>
0033 grav = zeros(nD,1); grav(nD) = -g;                                          <span class="comment">% gavity vector</span>
0034 <span class="keyword">for</span> mp = 1:nmp
0035    nIN = mpData(mp).nIN;                                                    <span class="co