10 February 2012 08:37:55 PM TOMS446_PRB C version Test the TOMS446 library. TEST01 Test CHEBY, which computes the Chebyshev series for several functions. Sin(x) Cos(x) Sin(2x) Cos(2x) X^5 0.000000 1.530395 0.000000 0.447782 0.000000 0.880101 0.000000 1.153450 -0.000000 0.625000 0.000000 -0.229807 0.000000 -0.705668 0.000000 -0.039127 0.000000 -0.257886 0.000000 0.312500 -0.000000 0.004953 -0.000000 0.067991 -0.000000 0.000500 -0.000000 0.014079 -0.000000 0.062500 0.000000 -0.000042 0.000000 -0.002405 0.000000 -0.000003 -0.000000 -0.000350 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000044 -0.000000 0.000000 0.000000 0.000010 0.000000 0.000000 TEST02 Test MLTPLY, which computes the product of two Chebyshev series. Multiply series for SIN(X) and COS(X) and compare with series for 1/2*SIN(2X). Sin(x) Cos(x) 1/2*Sin(2x) RESULT 4.93432e-17 1.5304 4.93432e-17 9.2625e-17 0.880101 9.86865e-17 0.576725 0.576725 1.23358e-16 -0.229807 1.11022e-16 1.74846e-16 -0.0391267 8.63507e-17 -0.128943 -0.128943 -4.93432e-17 0.00495328 -1.23358e-17 -3.31745e-17 0.000499515 -2.46716e-17 0.00703963 0.00703963 6.16791e-17 -4.18767e-05 4.31753e-17 2.93099e-17 -3.00468e-06 -2.46716e-17 -0.000174967 -0.000174946 0 1.87921e-07 1.85037e-17 -6.70784e-18 2.0997e-08 2.46716e-17 4.98469e-06 2.50038e-06 TEST03 Test ECHEB, which evaluates a Chebyshev series. Sin(x) -1 -0.841471 -0.841471 -0.6 -0.564642 -0.564642 -0.2 -0.198669 -0.198669 0.2 0.198669 0.198669 0.6 0.564642 0.564642 1 0.841471 0.841471 Cos(x) -1 0.540302 0.540302 -0.6 0.825336 0.825336 -0.2 0.980067 0.980067 0.2 0.980067 0.980067 0.6 0.825336 0.825336 1 0.540302 0.540302 Sin(2x) -1 -0.909297 -0.909302 -0.6 -0.932039 -0.932037 -0.2 -0.389418 -0.389423 0.2 0.389418 0.389423 0.6 0.932039 0.932037 1 0.909297 0.909302 Cos(2x) -1 -0.416147 -0.416147 -0.6 0.362358 0.362357 -0.2 0.921061 0.921061 0.2 0.921061 0.921061 0.6 0.362358 0.362357 1 -0.416147 -0.416147 x^5 -1 -1 -1 -0.6 -0.07776 -0.07776 -0.2 -0.00032 -0.00032 0.2 0.00032 0.00032 0.6 0.07776 0.07776 1 1 1 TEST04 Test EDCHEB, which evaluates the derivative of a Chebyshev series. Sin(x) -1 0.540302 0.540303 -0.6 0.825336 0.825336 -0.2 0.980067 0.980067 0.2 0.980067 0.980067 0.6 0.825336 0.825336 1 0.540302 0.540303 Cos(x) -1 0.841471 0.841471 -0.6 0.564642 0.564642 -0.2 0.198669 0.198669 0.2 -0.198669 -0.198669 0.6 -0.564642 -0.564642 1 -0.841471 -0.841471 Sin(2x) -1 -0.832294 -0.831887 -0.6 0.724716 0.724764 -0.2 1.84212 1.84211 0.2 1.84212 1.84211 0.6 0.724716 0.724764 1 -0.832294 -0.831887 Cos(2x) -1 1.81859 1.81858 -0.6 1.86408 1.86408 -0.2 0.778837 0.778828 0.2 -0.778837 -0.778828 0.6 -1.86408 -1.86408 1 -1.81859 -1.81858 x^5 -1 5 5 -0.6 0.648 0.648 -0.2 0.008 0.008 0.2 0.008 0.008 0.6 0.648 0.648 1 5 5 TEST05 Test DFRNT, which computes the Chebyshev series for the derivative of several functions. Chebyshev series for d/dx of: Sin(x) Cos(x) Sin(2x) Cos(2x) X^5 1.5304 5.67447e-16 0.895653 -3.45403e-16 3.75 8.38835e-16 -0.880101 2.31913e-15 -2.3069 1.97373e-16 -0.229807 3.70074e-16 -1.41125 -8.63507e-17 2.5 3.45403e-16 0.0391267 1.43095e-15 0.515775 -4.93432e-16 0.00495347 -1.4803e-16 0.136073 -1.23358e-16 0.625 7.40149e-16 -0.000499513 1.62833e-15 -0.0281566 -3.94746e-16 -4.16875e-05 9.86865e-17 -0.00471963 -9.86076e-32 -3.45403e-16 0 3.00673e-06 5.92119e-16 0.000701696 -3.94746e-16 3.77946e-07 4.44089e-16 0.000179449 7.77156e-16 0 0 0 0 0 0 TEST06 Test NTGRT, which computes the Chebyshev series for the indefinite integral of several functions. Chebyshev series for indefinite integral of: Sin(x) Cos(x) Sin(2x) Cos(2x) X^5 0 0 0 0 0 -3.70074e-17 0.880101 -6.16791e-17 0.576725 0 0.229807 3.08395e-18 0.352834 -3.39235e-17 0.078125 2.87836e-17 -0.0391267 4.11194e-17 -0.128943 3.08395e-17 -0.00495328 1.38778e-17 -0.0339957 2.31296e-18 0.03125 -1.11022e-17 0.000499515 -1.11022e-17 0.00703963 -1.23358e-18 4.18767e-05 0 0.00120243 3.59794e-18 0.00520833 4.40565e-18 -3.00461e-06 3.52452e-18 -0.000174908 1.76226e-18 -1.89105e-07 -3.08395e-18 -2.2494e-05 -6.16791e-18 -1.54198e-18 0 1.044e-08 2.05597e-18 2.43645e-06 -1.37065e-18 TOMS446_PRB Normal end of execution. 10 February 2012 08:37:55 PM