C 代码 使用以下方法评估基本函数和特殊函数 切比雪夫多项式.rar

03 March 2022 10:59:32 PM FN_TEST: C version FN evaluates special functions. I4_MACH_TEST: I4_MACH evaluates integer machine numbers. I4_MACH(1) = the standard input unit. I4_MACH(2) = the standard output unit. I4_MACH(3) = the standard punch unit. I4_MACH(4) = the standard error message unit. I4_MACH(5) = the number of bits per integer storage unit. I4_MACH(6) = the number of characters per integer storage unit. I4_MACH(7) = A, the base. I4_MACH(8) = S, the number of base A digits. I4_MACH(9) = A^S-1, the largest integer. I4_MACH(10) = B, the base. I4_MACH(11) = T, the number of single precision base B digits. I4_MACH(12) = EMIN, the smallest single precision exponent E. I4_MACH(13) = EMAX, the largest single precision exponent E. I4_MACH(14) = T, the number of double precision base B digits. I4_MACH(15) = EMIN, the smallest double precision exponent E. I4_MACH(16) = EMAX, the largest double precision exponent E. I I4_MACH(I) 1 5 2 6 3 7 4 6 5 32 6 4 7 2 8 31 9 2147483647 10 2 11 24 12 -125 13 128 14 53 15 -1021 16 1024 R8_ACOS_TEST: R8_ACOS computes the arc cosine function. X ARCCOS(X) R8_ACOS(X) Diff -0.1 1.67096 1.67096 2.22045e-16 0 1.5708 1.5708 0 0.1 1.47063 1.47063 0 0.2 1.36944 1.36944 2.22045e-16 0.3 1.2661 1.2661 2.22045e-16 0.4 1.15928 1.15928 0 0.5 1.0472 1.0472 0 0.6 0.927295 0.927295 0 0.7 0.795399 0.795399 0 0.8 0.643501 0.643501 1.11022e-16 0.9 0.451027 0.451027 5.55112e-17 1 0 0 0 R8_ACOSH_TEST: R8_ACOSH evaluates the arc hyperbolic cosine X ARCCOSH(X) R8_ACOSH(X) Diff 1 0 0 0 1.01 0.141304 0.141304 8.32667e-17 1.1 0.443568 0.443568 1.11022e-16 1.2 0.622363 0.622363 1.11022e-16 1.3 0.756433 0.756433 0 1.4 0.867015 0.867015 1.11022e-16 1.5 0.962424 0.962424 0 2 1.31696 1.31696 2.22045e-16 3 1.76275 1.76275 0 3.14159 1.81153 1.81153 0 4 2.06344 2.06344 4.44089e-16 5 2.29243 2.29243 0 10 2.99322 2.99322 0 100 5.29829 5.29829 0 1000 7.6009 7.6009 0 R8_AI_TEST: R8_AI evaluates the Airy function Ai(x) X AIRY_AI(X) R8_AI(X) Diff 0 0.355028 0.355028 0 0.1 0.329203 0.329203 0 0.2 0.303703 0.303703 0 0.3 0.278806 0.278806 0 0.4 0.254742 0.254742 0 0.5 0.231694 0.231694 0 0.6 0.2098 0.2098 2.77556e-17 0.7 0.189162 0.189162 5.55112e-17 0.8 0.169846 0.169846 5.55112e-17 0.9 0.151887 0.151887 5.55112e-17 1 0.135292 0.135292 2.77556e-17 R8_AID_TEST: R8_AID evaluates the derivative of the Airy Ai(x) function. X AIRY_AID(X) R8_AID(X) Diff 0 -0.258819 -0.258819 0 0.1 -0.25713 -0.25713 5.55112e-17 0.2 -0.252405 -0.252405 5.55112e-17 0.3 -0.245146 -0.245146 0 0.4 -0.235832 -0.235832 2.77556e-17 0.5 -0.224911 -0.224911 0 0.6 -0.212793 -0.212793 2.77556e-17 0.7 -0.199851 -0.199851 2.77556e-17 0.8 -0.186413 -0.186413 0 0.9 -0.172764 -0.172764 0 1 -0.159147 -0.159147 0 R8_AINT_TEST: R8_AINT rounds an R8 towards 0. X AINT(X) R8_AINT(X) Diff -2.01 -2 -2 0 -1.99 -1 -1 0 -1.5 -1 -1 0 -1.1 -1 -1 0 -1.01 -1 -1 0 -1 -1 -1 0 -0.99 0 -0 0 -0.9 0 -0 0 -0.51 0 -0 0 -0.5 0 -0 0 -0.49 0 -0 0 -0.01 0 -0 0 0 0 0 0 0.01 0 0 0 0.49 0 0 0 0.5 0 0 0 0.51 0 0 0 0.9 0 0 0 0.99 0 0 0 1 1 1 0 1.01 1 1 0 1.1 1 1 0 1.5 1 1 0 1.99 1 1 0 2.01 2 2 0 R8_ASIN_TEST: R8_ASIN evaluates the arcsine. X ARCSIN(X) R8_ASIN(X) Diff -0.1 -0.100167 -0.100167 0 0 0 0 0 0.1 0.100167 0.100167 0 0.2 0.201358 0.201358 0 0.3 0.304693 0.304693 0 0.4 0.411517 0.411517 5.55112e-17 0.5 0.523599 0.523599 1.11022e-16 0.6 0.643501 0.643501 0 0.7 0.775397 0.775397 1.11022e-16 0.8 0.927295 0.927295 1.11022e-16 0.9 1.11977 1.11977 0 1 1.5708 1.5708 0 R8_ASINH_TEST: R8_ASINH evaluates the hyperbolic arcsine X ARCSINH(X) R8_ASINH(X) Diff -5 -2.31244 -2.31244 0 -1 -0.881374 -0.881374 0 0 0 0 0 0.1 0.0998341 0.0998341 1.38778e-17 0.2