In this thesis, we study the contractivity preserving, high order, time discretization methods for solving non-stiff ordinary differential equations. We construct a class of one-step, explicit, contractivity preserving, multi-stage, multi-derivative, Hermite-Birkhoff-Taylor methods of order p=5,6, ..., 15, that we denote by CPHBT, with nonnegative coefficients by casting s-stage Runge-Kutta methods of order 4 and 5 with Taylor methods of order p-3 and p-4, respectively. The constructed CPHBT methods are implemented using an efficient variable step algorithm and are compared to other well-known methods on a variety of initial value problems. The results show that CPHBT methods have larger regions of absolute stability, require less function evaluations and hence they require less CPU time to achieve the same accuracy requirements as other methods in the literature. Also, we show that the contractivity preserving property of CPHBT is very efficient in suppressing the effect of the propagation of discretization errors when a long-term integration of a standard N-body problem is considered. The formulae of 49 CPHBT methods of various orders are provided in Butcher form.