new script for calculating a reference for a colleague's finite difference implementation