An Efficient Approximation Method for Phase-Type Distributions