Numerical method for the solution of the transient behavior of bipolar semiconductor devices