Type Help(); to begin. **************************************************** **************************************************** One trial of runTest(10,7,100) gives: 8.2100 One trial of runTestPGF(10,7,100,x) gives: 3 10 11 9 1 8 3 7 2 6 3 5 1 4 1 3 -- x + -- x + - x + -- x + -- x + --- x + --- x + --- x 10 50 5 20 25 100 100 100 and [seq(moms(f,x,i), i=1..6)] gives the moments: [8.330, 1.569, -0.8657, 3.382, -7.652, 23.90] This took 0.209000 seconds. **************************************************** **************************************************** One trial of runTest(100,70,100) gives: 97.690 One trial of runTestPGF(100,70,100,x) gives: 7 100 13 99 3 98 9 97 2 96 3 95 3 94 -- x + -- x + -- x + --- x + -- x + --- x + -- x 25 50 25 100 25 100 50 1 93 1 92 1 91 1 86 + -- x + --- x + -- x + --- x 25 100 50 100 and [seq(moms(f,x,i), i=1..6)] gives the moments: [97.72, 2.593, -1.636, 6.243, -22.82, 61.43] This took 6.990000 seconds. **************************************************** **************************************************** One trial of runTest(1000,900,100) gives: 991.92 One trial of runTestPGF(1000,900,100,x) gives: 3 986 3 1000 3 998 3 997 11 996 1 995 --- x + -- x + --- x + -- x + --- x + -- x 100 50 100 25 100 25 2 994 7 993 1 992 1 991 3 990 1 988 + -- x + --- x + -- x + -- x + -- x + -- x 25 100 50 20 50 50 1 987 1 985 1 984 1 982 1 980 1 979 + --- x + -- x + -- x + --- x + --- x + -- x 100 50 20 100 100 50 1 976 1 975 1 973 1 972 1 968 + --- x + --- x + --- x + -- x + --- x 100 100 100 50 100 1 958 7 999 1 989 + --- x + --- x + -- x 100 100 20 and [seq(moms(f,x,i), i=1..6)] gives the moments: [991.4, 7.763, -1.645, 6.230, -61.53, -12490.] This took 257.891000 seconds. **************************************************** ****************************************************