Differential equations are involved in many felds. Often these cannot be solved exactly, so they must be approximated. To do this, we derive a fourth order finite difference scheme for non-uniform meshes. This is paired ...