1
2exp1:=441*x^7+780*x^6-2861*x^5+4085*x^4+7695*x^3+3713*x^2-43253*x+24500;
3
4exp2:=9*x^6+6*x^5-65*x^4+20*x^3+135*x^2-154*x+49;
5
6aa:=7*x^13+10*x^8+4*x^7-7*x^6-4*x^3-4*x^2+3*x+3;
7bb:=x^14-2*x^8-2*x^7-2*x^4-4*x^3-x^2+2*x+1;
8
9% example 2.14
10
11a:=48*x^3-84*x^2+42*x-36;
12
13b:=-4*x^3-10*x^2+44*x-30;
14
15
16% square free ex
17
18c:=x^8-2*x^6+2*x^2-1;
19
20%square_free(a,x);
21
22%eval_sq_free(ws,x);
23
24
25% hr example
26
27pu:=441*x^7+780*x^6-2861*x^5+4085*x^4+7695*x^3+3713*x^2-43253*x+24500;
28
29qu:=9*x^6+6*x^5-65*x^4+20*x^3+135*x^2-154*x+49;
30
31%makemon(p,q,x);
32
33%pf(pu/qu,x);
34%trst hr;
35
36
37k:=36*x^6+126*x^5+183*x^4+(13807/6)*x^3-407*x^2-(3242/5)*x+(3044/15);
38l:=(x^2+(7/6)*x+(1/3))^2*(x-(2/5))^3;
39%trst hr;
40
41ratint(k,l,x);
42
43ratint(exp1,exp2,x);
44
45ratint(1,x^2,x);
46
47ratint(1,x^2+1,x);
48
49ratint(1,x^3+x+1,x);
50
51ratint(1,x^5+1,x);
52
53end;
54