This paper shows how the multivariate finite time ruin probability function, in a phase-type environment, inherits the phase-type structure and can be efficiently approximated with only one Laplace transform inversion. From a theoretical point of view, we also provide below a generalization of Thorin’s formula (1971) for the double Laplace transform of the finite time ruin probability, by considering also the deficit at ruin; the model is that of a Sparre Andersen (renewal) risk process with phase-type interarrival times . In the case when the claims distribution is of phase-type as well, we obtain also an alternative formula for the single Laplace transform in time (or “exponentially killed probability’’), in terms of the roots with positive real part of the Lundberg’s equations, which complements Asmussen’s representation (1992) in terms of the roots with negative real part.