The fluxes on the faces of the control region are approximated by difference quotients in the same way as in the one dimensional case. That is,

If we use to denote the FDE approximation to , we then obtain the following finite difference methods.