function numEquations = computeEquationNumber(nDof, element_data)

nElements = length(element_data); 
numEquations = cell(nElements, 1); 

for i = 1:nElements
    elementConnectivity = element_data{i}{2};
    numEquationsElem = [];
    for j = 1:length(elementConnectivity)
        for k = 1:nDof
            % Compute the equation number
            numEquationsElem = [numEquationsElem, (elementConnectivity(j)-1) * nDof + (k - 1)+1];
        end
    end
    numEquations{i} = numEquationsElem;
end

end