In this paper, a numerical solution is presented for the generalized BLT equation with inhomogeneous transmission lines. In particular, a fully automatized method is developed for the computation of the transmission line propagation matrices and the junction scattering matrices from network structural specification and from transmission line characteristic parameters. Two numerical examples are presented, one for a tree-shaped network, and the other for a network involving a circuit loop.