Since most ordinary differential equations cannot be solved in terms of explicit formulas, numerical methods are required to solve them approximately. Runge-Kutta methods are the standard for this. we present a method based on quadratic Taylor expansion of the differential equation. The motivation is to improve accuracy in certain ecological models. In particular, logistic models with harvesting are solved error-free by our methods