More explicit form of Eq. We use fourth order predictor-corrector. To solve Eq. The time step is chosen to be small enough that the corrector part need only be applied once. The 4th order predictor-corrector method requires us to keep zeroth though 4th derivatives of and .
Similar equations apply for and .
20 continue Figure 4. Similarly for and . Similarly for and . So before applying the corrector cycle, we need to compute , and . 80 continue Figure 6.