7 October 2012 10:04:47.658 AM TEST_INTERP_1D_TEST FORTRAN77 version Test the TEST_INTERP_1D library. The R8LIB library is needed. TEST01 Print the title of each function. There are 8 functions available: Index Title 1 f(x) = steps -1/2/1 at [0,1/3], [1/3,4/5], [4/5,1]. 2 f(x) = (1-3x), x < 1/3 (6x-2) if 1/3 < x 3 f(x) = x (10*x-1) (5x-2) (5x-2) (4x-3.4) (x-1) 4 f(x) = atan ( 40 * x - 15 ) 5 f(x) = cos(7*x)+5*cos(11.2*x)-2*cos(14*x)+5*cos(31.5*x)+7*cos(63*x). 6 f(x) = exp ( - ( 4*x-1 )^2 ) 7 f(x) = exp ( 2 x ) if x < 0.5, 0 otherwise 8 f(x) = 1 / ( 1 + ( 10 * (x-1/2) )^2 ) TEST_INTERP_1D_TEST02 Use P00_F to sample each function. X, Y for problem 1 1: 0.00000 -1.00000 2: 0.100000 -1.00000 3: 0.200000 -1.00000 4: 0.300000 -1.00000 5: 0.400000 2.00000 6: 0.500000 2.00000 7: 0.600000 2.00000 8: 0.700000 2.00000 9: 0.800000 2.00000 10: 0.900000 1.00000 11: 1.00000 1.00000 X, Y for problem 2 1: 0.00000 1.00000 2: 0.100000 0.700000 3: 0.200000 0.400000 4: 0.300000 0.100000 5: 0.400000 0.400000 6: 0.500000 1.00000 7: 0.600000 1.60000 8: 0.700000 2.20000 9: 0.800000 2.80000 10: 0.900000 3.40000 11: 1.00000 4.00000 X, Y for problem 3 1: 0.00000 -0.00000 2: 0.100000 0.00000 3: 0.200000 0.416000 4: 0.300000 0.231000 5: 0.400000 0.00000 6: 0.500000 0.350000 7: 0.600000 1.20000 8: 0.700000 1.70100 9: 0.800000 0.896000 10: 0.900000 -0.900000 11: 1.00000 0.00000 X, Y for problem 4 1: 0.00000 -1.50423 2: 0.100000 -1.48014 3: 0.200000 -1.42890 4: 0.300000 -1.24905 5: 0.400000 0.785398 6: 0.500000 1.37340 7: 0.600000 1.46014 8: 0.700000 1.49402 9: 0.800000 1.51204 10: 0.900000 1.52321 11: 1.00000 1.53082 X, Y for problem 5 1: 0.00000 16.0000 2: 0.100000 4.60251 3: 0.200000 10.9479 4: 0.300000 -2.41286 5: 0.400000 8.33648 6: 0.500000 3.41326 7: 0.600000 17.0369 8: 0.700000 4.07749 9: 0.800000 7.82541 10: 0.900000 -3.02927 11: 1.00000 13.3791 X, Y for problem 6 1: 0.00000 0.367879 2: 0.100000 0.697676 3: 0.200000 0.960789 4: 0.300000 0.960789 5: 0.400000 0.697676 6: 0.500000 0.367879 7: 0.600000 0.140858 8: 0.700000 0.391639E-01 9: 0.800000 0.790705E-02 10: 0.900000 0.115923E-02 11: 1.00000 0.123410E-03 X, Y for problem 7 1: 0.00000 1.00000 2: 0.100000 1.49182 3: 0.200000 2.22554 4: 0.300000 3.32012 5: 0.400000 4.95303 6: 0.500000 0.00000 7: 0.600000 0.00000 8: 0.700000 0.00000 9: 0.800000 0.00000 10: 0.900000 0.00000 11: 1.00000 0.00000 X, Y for problem 8 1: 0.00000 0.384615E-01 2: 0.100000 0.588235E-01 3: 0.200000 0.100000 4: 0.300000 0.200000 5: 0.400000 0.500000 6: 0.500000 1.00000 7: 0.600000 0.500000 8: 0.700000 0.200000 9: 0.800000 0.100000 10: 0.900000 0.588235E-01 11: 1.00000 0.384615E-01 TEST_INTERP_1D_TEST Normal end of execution. 7 October 2012 10:04:47.659 AM