“…The exact solutions of these systems are difficult to find. In finding the approximate solutions, different approaches in the literature have been used such as homotopy perturbation method (HPM), 1 homotopy analysis method (HAM), 2 differential transformation method (DTM), 3 Genocchi polynomials (GP), 4 reproducing kernel Hilbert space method (RKHSM), 5 Chebyshev wavelets method (CWM), 6 fixed point method (FPM), 7 Sinc-collocation method (SCM), 8 Haar functions method (HFM), 9 Adomian decomposition method (ADM), 10 and relaxed Monte Carlo method (RMCM). 11 In the present attempt, optimal homotopy asymptotic method (OHAM) is used, which was introduced by Marinca et al 12,13 for the solution of differential equations.…”