This paper applies the Gear-4 algorithm to constitute the companion models lor the circuit elements and electromagnetic relays commonly used in relay protective devices. The interface computation, which describes the node equation of companion model and differential equation of mechanical moving part for the electromagnetic Jelay, is also presented. For the reason mentioned above, it makes possible* to solve more accurately the transient process of relay protective devices including electromagnetic relay using the companion model method. The simplicity and effectiveness of interface computation have been proved by the example calcultaion.