The introduction of the bicausality concept in the bond graph language has allowed new analytical methodologies, for instance in the context of model inversion, mechatronic system sizing and control. The bicausality concept is here applied for solving the equilibrium state of a mechatronic system. We propose a new method, which permits us to determine the size of the equilibrium set and the algebraic system to be solved. The proposed method is applied to linear systems in a first step, and a generalization is also given for some non-linear systems. Several examples are included in order to explain the method.