In the present work, numerical scheme has been devised and implemented for the solution of problems of complex one-dimensional flows involving simultaneous action of driving potentials, namely area variation, heat addition, friction and mass injection.