1on ranum; 2 3f := (x**2-2)*x*(x-5); 4 5rl := isolate f; 6 7for each r in rl collect sub(x=r, f); 8 9r1 := refine(part(rl, 1), 10); 10 11on rarat; 12 13r1; 14 15refine(r1, 100); 16 17off rarat; 18 19x1 := ra(f, 1, 2); 20 21off nat; 22 23x1; 24 25on nat; 26 27x2 := ra(f, 4, 10); 28 29isolate f; 30 31x1-3; 32 33x1+1/7; 34 35x1/5; 36 37-x1; 38 39x1 + x2; 40 411/x1; 42 43x1 - x2; 44 45x1 * x2; 46 47x1^4; 48 49x1**10 - (x1*x1*x1*x1*x1*x1*x1*x1*x1*x1); 50 51wilk := 8388608*x**20 + 1761607681*x**19 + 172931153920*x**18 + 5210543221964800*x**17 + 447347234439168*x**16 + 5314028108264898560*x**15 + 336985244869591040*x**14 + 546342720331186176000*x**13 + 94877480085669019648*x**12 + 551137370949952460554240 *x**11 + 10968398649699241820160*x**10 + 5685079777790228273561600*x**9 + 528740774622641958944768 *x**8 + 572611655889692786813829120*x**7 + 10122095419974470210682880*x**6 + 5830198816984091441338777600*x**5 + 67426052557934862488567808*x**4 + 59107969196810523545855590400*x**3 + 115794329499468438700032000 *x**2 + 6073425049924762651852800000*x + 20408661249006627717120000$ 61 62sub(x=x1, wilk); 63 64sub(x=x2, wilk); 65 66isolatingivl wilk; 67 68on rahidepoly; 69 70wrl := isolate wilk; 71 72for each x in wrl collect refine(x, 10); 73 74off rahidepoly; 75 76on errcont; 77 78symbolic ra_wrappertest(6, 42); 79 80symbolic ra_wrappertest(7, 42); 81 82off msg; 83 84ra(f, 1, 10); 85 86ra(f, 1/2, 1); 87 88ra(f, 1, 5); 89 90ra(f, 5, 10); 91 92on msg; 93 94off errcont; 95 96end; 97