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