A numerical approximation method for the Fisher-Rao distance between multivariate normal distributions