In our algorithm, we use Heun’s method to derive thenumerical solution of (8), (9), (10), (11), (12), (13), and (14)by the following reasons: Euler’s method is simple but hasunstable performance in terms of accuracy.