Numerical computation of arbitrary order transfer maps and reconstructive correction of aberrations in the large acceptance spectrometer MAGNEX