MIBM Higher Order Central Differencing

[danieljvickers]

Problem it solves

Much of the moving immersed boundary code uses second-order central differences, which is obviously sub-ideal.

Proposed solution

Replace the hard-coded central differences with a method that uses the fd_coefficients in the code