An implementation of global error estimation algorithm in a transient analysis program has been described. This algorithm does not affect the efficiency of the program and may be useful in the analysis of a wide class of electronics circuits. The global error estimate may be used for controlling the step size and order of the backward differentiation formulas. An algorithm for changing the step size and order based on the global error estimate has been described and compared with two approaches commonly used in circuit analysis programs, i. e. constant error per unit step and constant error per step methods.