Allow functions for Field.coordinate_transform

This lets you quickly switch to polar coordinate systems:

    field.coordinate_transform = lambda x: sympy.Matrix((x.norm(), sympy.atan2(*x) / (2 * sympy.pi) * field.shape[1]))