          % id of the nodes of the edge
                        idNds = [NdBorders(j); NdBorders(j+1)];
            
                        % Equivalent nodal load
                        if quadMesh == false
                            feq = [0.5*p*l;0.5*p*l];
                        else
                            feq = [p*l;p*l;4.0*p*l]/6.0;
                            idNds = [idNds; this.ELEM{el}(j+nLinNodes)];
                        end
            
                        % Add contribution to the