A new procedure that generates the transient solution of the first moment of the state of a Markovian queueing network with state-dependent arrivals, services, and routeing is developed. The procedure involves defining a partial differential equation that relates an approximate multivariate cumulant generating function to the intensity functions of the network. The partial differential equation then yields a set of ordinary differential equations which are numerically solved to obtain the first moment.