Numerical strategies for the system of IVPs using extended fourth order RK methods based on a variety of means