A new way of formulating a multipoint Padé approximant of a linear system transfer function is given. It is seen to be very flexible and computationally efficient. The expansion points can be a mixture of real, complex, purely imaginary and multiple points with no complex arithmetic required. The underlying idea of the method is to convert the rational approximation problem into one of polynomial approximation. This results in simple matrix multiplications and the solution of a set of linear equations to derive the reduced model. Examples are given to demonstrate the application of the technique.