Annotation of OpenXM_contrib/pari/src/test/32/polyser, Revision 1.1
1.1 ! maekawa 1: realprecision = 38 significant digits
! 2: echo = 1 (on)
! 3: ? apol=x^3+5*x+1
! 4: x^3 + 5*x + 1
! 5: ? changevar(x+y,[z,t])
! 6: y + z
! 7: ? deriv((x+y)^5,y)
! 8: 5*x^4 + 20*y*x^3 + 30*y^2*x^2 + 20*y^3*x + 5*y^4
! 9: ? ((x+y)^5)'
! 10: 5*x^4 + 20*y*x^3 + 30*y^2*x^2 + 20*y^3*x + 5*y^4
! 11: ? dz=vector(30,k,1);dd=vector(30,k,k==1);dm=dirdiv(dd,dz)
! 12: [1, -1, -1, 0, -1, 1, -1, 0, 0, 1, -1, 0, -1, 1, 1, 0, -1, 0, -1, 0, 1, 1, -
! 13: 1, 0, 0, 1, 0, 0, -1, -1]
! 14: ? direuler(s=1,40,1+s*X+s^2*X)
! 15: [1, 6, 12, 0, 30, 72, 56, 0, 0, 180, 132, 0, 182, 336, 360, 0, 306, 0, 380,
! 16: 0, 672, 792, 552, 0, 0, 1092, 0, 0, 870, 2160, 992, 0, 1584, 1836, 1680, 0,
! 17: 1406, 2280, 2184, 0]
! 18: ? dirmul(abs(dm),dz)
! 19: [1, 2, 2, 2, 2, 4, 2, 2, 2, 4, 2, 4, 2, 4, 4, 2, 2, 4, 2, 4, 4, 4, 2, 4, 2,
! 20: 4, 2, 4, 2, 8]
! 21: ? zz=yy;yy=xx;eval(zz)
! 22: xx
! 23: ? factorpadic(apol,7,8)
! 24:
! 25: [(1 + O(7^8))*x + (6 + 2*7^2 + 2*7^3 + 3*7^4 + 2*7^5 + 6*7^6 + O(7^8)) 1]
! 26:
! 27: [(1 + O(7^8))*x^2 + (1 + 6*7 + 4*7^2 + 4*7^3 + 3*7^4 + 4*7^5 + 6*7^7 + O(7^8
! 28: ))*x + (6 + 5*7 + 3*7^2 + 6*7^3 + 7^4 + 3*7^5 + 2*7^6 + 5*7^7 + O(7^8)) 1]
! 29:
! 30: ? factorpadic(apol,7,8,1)
! 31:
! 32: [(1 + O(7^8))*x + (6 + 2*7^2 + 2*7^3 + 3*7^4 + 2*7^5 + 6*7^6 + O(7^8)) 1]
! 33:
! 34: [(1 + O(7^8))*x^2 + (1 + 6*7 + 4*7^2 + 4*7^3 + 3*7^4 + 4*7^5 + 6*7^7 + O(7^8
! 35: ))*x + (6 + 5*7 + 3*7^2 + 6*7^3 + 7^4 + 3*7^5 + 2*7^6 + 5*7^7 + O(7^8)) 1]
! 36:
! 37: ? intformal(sin(x),x)
! 38: 1/2*x^2 - 1/24*x^4 + 1/720*x^6 - 1/40320*x^8 + 1/3628800*x^10 - 1/479001600*
! 39: x^12 + 1/87178291200*x^14 - 1/20922789888000*x^16 + O(x^17)
! 40: ? intformal((-x^2-2*a*x+8*a)/(x^4-14*x^3+(2*a+49)*x^2-14*a*x+a^2),x)
! 41: (x + a)/(x^2 - 7*x + a)
! 42: ? newtonpoly(x^4+3*x^3+27*x^2+9*x+81,3)
! 43: [2, 2/3, 2/3, 2/3]
! 44: ? padicappr(apol,1+O(7^8))
! 45: [1 + 6*7 + 4*7^2 + 4*7^3 + 3*7^4 + 4*7^5 + 6*7^7 + O(7^8)]
! 46: ? padicappr(x^3+5*x+1,Mod(x*(1+O(7^8)),x^2+x-1))
! 47: [Mod((1 + 3*7 + 3*7^2 + 4*7^3 + 4*7^4 + 4*7^5 + 2*7^6 + 3*7^7 + O(7^8))*x +
! 48: (2*7 + 6*7^2 + 6*7^3 + 3*7^4 + 3*7^5 + 4*7^6 + 5*7^7 + O(7^8)), x^2 + x - 1)
! 49: ]~
! 50: ? Pol(sin(x),x)
! 51: -1/1307674368000*x^15 + 1/6227020800*x^13 - 1/39916800*x^11 + 1/362880*x^9 -
! 52: 1/5040*x^7 + 1/120*x^5 - 1/6*x^3 + x
! 53: ? Pol([1,2,3,4,5],x)
! 54: x^4 + 2*x^3 + 3*x^2 + 4*x + 5
! 55: ? Polrev([1,2,3,4,5],x)
! 56: 5*x^4 + 4*x^3 + 3*x^2 + 2*x + 1
! 57: ? polcoeff(sin(x),7)
! 58: -1/5040
! 59: ? polcyclo(105)
! 60: x^48 + x^47 + x^46 - x^43 - x^42 - 2*x^41 - x^40 - x^39 + x^36 + x^35 + x^34
! 61: + x^33 + x^32 + x^31 - x^28 - x^26 - x^24 - x^22 - x^20 + x^17 + x^16 + x^1
! 62: 5 + x^14 + x^13 + x^12 - x^9 - x^8 - 2*x^7 - x^6 - x^5 + x^2 + x + 1
! 63: ? pcy=polcyclo(405)
! 64: x^216 - x^189 + x^135 - x^108 + x^81 - x^27 + 1
! 65: ? pcy*pcy
! 66: x^432 - 2*x^405 + x^378 + 2*x^351 - 4*x^324 + 4*x^297 - x^270 - 4*x^243 + 7*
! 67: x^216 - 4*x^189 - x^162 + 4*x^135 - 4*x^108 + 2*x^81 + x^54 - 2*x^27 + 1
! 68: ? poldegree(x^3/(x-1))
! 69: 2
! 70: ? poldisc(x^3+4*x+12)
! 71: -4144
! 72: ? poldiscreduced(x^3+4*x+12)
! 73: [1036, 4, 1]
! 74: ? polinterpolate([0,2,3],[0,4,9],5)
! 75: 25
! 76: ? polisirreducible(x^5+3*x^3+5*x^2+15)
! 77: 0
! 78: ? pollegendre(10)
! 79: 46189/256*x^10 - 109395/256*x^8 + 45045/128*x^6 - 15015/128*x^4 + 3465/256*x
! 80: ^2 - 63/256
! 81: ? zpol=0.3+pollegendre(10)
! 82: 46189/256*x^10 - 109395/256*x^8 + 45045/128*x^6 - 15015/128*x^4 + 3465/256*x
! 83: ^2 + 0.053906249999999999999999999999999999999
! 84: ? polrecip(3*x^7-5*x^3+6*x-9)
! 85: -9*x^7 + 6*x^6 - 5*x^4 + 3
! 86: ? polresultant(x^3-1,x^3+1)
! 87: 8
! 88: ? polresultant(x^3-1.,x^3+1.,,1)
! 89: 8.0000000000000000000000000000000000000
! 90: ? polroots(x^5-5*x^2-5*x-5)
! 91: [2.0509134529831982130058170163696514536 + 0.E-38*I, -0.67063790319207539268
! 92: 663382582902335603 + 0.84813118358634026680538906224199030917*I, -0.67063790
! 93: 319207539268663382582902335603 - 0.84813118358634026680538906224199030917*I,
! 94: -0.35481882329952371381627468235580237077 + 1.39980287391035466982975228340
! 95: 62081964*I, -0.35481882329952371381627468235580237077 - 1.399802873910354669
! 96: 8297522834062081964*I]~
! 97: ? polroots(x^4-1000000000000000000000,1)
! 98: [-177827.94100389228012254211951926848447 + 0.E-38*I, 177827.941003892280122
! 99: 54211951926848447 + 0.E-38*I, 6.6530622500127354998594589316364200753 E-111
! 100: + 177827.94100389228012254211951926848447*I, 6.65306225001273549985945893163
! 101: 64200753 E-111 - 177827.94100389228012254211951926848447*I]~
! 102: ? polrootsmod(x^16-1,41)
! 103: [Mod(1, 41), Mod(3, 41), Mod(9, 41), Mod(14, 41), Mod(27, 41), Mod(32, 41),
! 104: Mod(38, 41), Mod(40, 41)]~
! 105: ? polrootspadic(x^4+1,41,6)
! 106: [3 + 22*41 + 27*41^2 + 15*41^3 + 27*41^4 + 33*41^5 + O(41^6), 14 + 20*41 + 2
! 107: 5*41^2 + 24*41^3 + 4*41^4 + 18*41^5 + O(41^6), 27 + 20*41 + 15*41^2 + 16*41^
! 108: 3 + 36*41^4 + 22*41^5 + O(41^6), 38 + 18*41 + 13*41^2 + 25*41^3 + 13*41^4 +
! 109: 7*41^5 + O(41^6)]~
! 110: ? polsturm(zpol)
! 111: 4
! 112: ? polsturm(zpol,0.91,1)
! 113: 1
! 114: ? polsylvestermatrix(a2*x^2+a1*x+a0,b1*x+b0)
! 115:
! 116: [a2 b1 0]
! 117:
! 118: [a1 b0 b1]
! 119:
! 120: [a0 0 b0]
! 121:
! 122: ? polsym(x^17-1,17)
! 123: [17, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 17]~
! 124: ? poltchebi(10)
! 125: 512*x^10 - 1280*x^8 + 1120*x^6 - 400*x^4 + 50*x^2 - 1
! 126: ? polzagier(6,3)
! 127: 4608*x^6 - 13824*x^5 + 46144/3*x^4 - 23168/3*x^3 + 5032/3*x^2 - 120*x + 1
! 128: ? serconvol(sin(x),x*cos(x))
! 129: x + 1/12*x^3 + 1/2880*x^5 + 1/3628800*x^7 + 1/14631321600*x^9 + 1/1448500838
! 130: 40000*x^11 + 1/2982752926433280000*x^13 + 1/114000816848279961600000*x^15 +
! 131: O(x^16)
! 132: ? serlaplace(x*exp(x*y)/(exp(x)-1))
! 133: 1 + (y - 1/2)*x + (y^2 - y + 1/6)*x^2 + (y^3 - 3/2*y^2 + 1/2*y)*x^3 + (y^4 -
! 134: 2*y^3 + y^2 - 1/30)*x^4 + (y^5 - 5/2*y^4 + 5/3*y^3 - 1/6*y)*x^5 + (y^6 - 3*
! 135: y^5 + 5/2*y^4 - 1/2*y^2 + 1/42)*x^6 + (y^7 - 7/2*y^6 + 7/2*y^5 - 7/6*y^3 + 1
! 136: /6*y)*x^7 + (y^8 - 4*y^7 + 14/3*y^6 - 7/3*y^4 + 2/3*y^2 - 1/30)*x^8 + (y^9 -
! 137: 9/2*y^8 + 6*y^7 - 21/5*y^5 + 2*y^3 - 3/10*y)*x^9 + (y^10 - 5*y^9 + 15/2*y^8
! 138: - 7*y^6 + 5*y^4 - 3/2*y^2 + 5/66)*x^10 + (y^11 - 11/2*y^10 + 55/6*y^9 - 11*
! 139: y^7 + 11*y^5 - 11/2*y^3 + 5/6*y)*x^11 + (y^12 - 6*y^11 + 11*y^10 - 33/2*y^8
! 140: + 22*y^6 - 33/2*y^4 + 5*y^2 - 691/2730)*x^12 + (y^13 - 13/2*y^12 + 13*y^11 -
! 141: 143/6*y^9 + 286/7*y^7 - 429/10*y^5 + 65/3*y^3 - 691/210*y)*x^13 + (y^14 - 7
! 142: *y^13 + 91/6*y^12 - 1001/30*y^10 + 143/2*y^8 - 1001/10*y^6 + 455/6*y^4 - 691
! 143: /30*y^2 + 7/6)*x^14 + O(x^15)
! 144: ? serreverse(tan(x))
! 145: x - 1/3*x^3 + 1/5*x^5 - 1/7*x^7 + 1/9*x^9 - 1/11*x^11 + 1/13*x^13 - 1/15*x^1
! 146: 5 + O(x^16)
! 147: ? subst(sin(x),x,y)
! 148: y - 1/6*y^3 + 1/120*y^5 - 1/5040*y^7 + 1/362880*y^9 - 1/39916800*y^11 + 1/62
! 149: 27020800*y^13 - 1/1307674368000*y^15 + O(y^16)
! 150: ? subst(sin(x),x,x+x^2)
! 151: x + x^2 - 1/6*x^3 - 1/2*x^4 - 59/120*x^5 - 1/8*x^6 + 419/5040*x^7 + 59/720*x
! 152: ^8 + 13609/362880*x^9 + 19/13440*x^10 - 273241/39916800*x^11 - 14281/3628800
! 153: *x^12 - 6495059/6227020800*x^13 + 69301/479001600*x^14 + 26537089/1188794880
! 154: 00*x^15 + O(x^16)
! 155: ? taylor(y/(x-y),y)
! 156: (O(y^16)*x^15 + y*x^14 + y^2*x^13 + y^3*x^12 + y^4*x^11 + y^5*x^10 + y^6*x^9
! 157: + y^7*x^8 + y^8*x^7 + y^9*x^6 + y^10*x^5 + y^11*x^4 + y^12*x^3 + y^13*x^2 +
! 158: y^14*x + y^15)/x^15
! 159: ? variable(name^4-other)
! 160: name
! 161: ? getheap
! 162: [61, 7111]
! 163: ? print("Total time spent: ",gettime);
! 164: Total time spent: 258
! 165: ? \q
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to polyser CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / pari / src / test / 32