Finite element equations for electromagnetic fields are examined, in particular nodal elements using scalar potential formulation and edge elements for vector potential formulation. It is shown how the equations usually obtained via variational approach may be more conveniently derived using integral methods employing a geometrical description of the interpolating functions of edge and facet elements. Moreover, the resultant equations describe the equivalent multi-branch circuit models