 - edgeY(1);
                        l = sqrt(dx*dx + dy*dy);

                        <span class="comment">% id of the nodes of the edge</span>
                        idNds = [NdBorders(j); NdBorders(j+1)];

                        <span class="comment">% Equivalent nodal load</span>
                        <span class="keyword">if</span> quadMesh == false
                            feq = [0.5*p*l;0.5*p*l];
                        <span class="keyword">else</span>
                            feq = [p*l;