An automatic algorithm for accurate numerical gradient calculations has been developed. The algorithm is based on both finite differences and Chebyshev interpolation approximations. The novelty of the method is an automated ...