Runge Kutta perturbation method in one and two dimensions