Reconstruction of MIMO nonlinear differential equation models from the generalised frequency response function matrix