Symbolic derivation of finite difference approximations for the three‐dimensional Poisson equation