“…For the sake of solving system (13), which is a linear system with nonconstant coefficients which is difficult to integrate, so the coefficient functions are approximated using Lagrange interpolation polynomial of degree 5, and hence system (13) will be reduced to: Then the solution is: π₯ 1 (π‘) = 3.176 π β 0.00292542π‘ 6 + 0.0969738 π‘ 5 β0.70755π‘ 4 +2.10439 π‘ 3 β2.49617π‘ 2 β2.49617t π₯ 2 (π‘) = 5.7217 π β 0.00478621 π‘ 6 + 0.0893614π‘ 5 β 0.558928π‘ 4 +1.44325π‘ 3 β0.798714π‘ 2 β3.91471t } The initial conditions for the third time-step interval [2,3] are supposed to be the previous step solutions which are given by equations ( 14), and similar approach is used as it is followed over the time-step interval [1,2], to get the next solutions: Similarly, we can proceed to evaluate the solutions and the control functions over the time-step intervals [3,4] and [4,5] which are sketched in Figures 7-10, respectively. Finally, the solution for the original system (7) for all conducted time step intervals are collected and sketched in Figure 11, as well as, the controller functions π’ 1 and π’ 2 are also collected in Figure 12.…”