Precise Numerical Differentiation of Thermodynamic Functions with Multicomplex Variables