“…Some of these models do not always have theoretical solutions, thus numerical methods are often employed to solve them. Researchers in most cases always use method of reduction of higher order ODEs into system of first order ODEs to solve (1). This technique though quite good, is bedeviled with many problems such tediousness, complexity of the method, waste of time, and the need for large computer storage memories because of too many auxiliary functions, etc.…”