High-Order Compact Implicit Difference Methods For Parabolic Equations in Geodynamo Simulation

Figure 1

(a) Solutions of the benchmark problem 1. (b) Solution errors of the sixth- and fourth-order compact schemes. The total grid points were 𝑁=100 and the parameter was 𝜖=0.01.

(a)

(b)