{VERSION 4 0 "SUN SPARC SOLARIS" "4.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 
0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 264 "signchanges:=proc(p
)\n  local c,d,i,s,x;\n  c:=0;\n  x:=indets(p)[1];\n  d:=degree(p,x);
\n  s:=lcoeff(p,x);\n  for i from d-1 to 0 by -1\n  do\n    if coeff(p
,x,i) <> 0\n    then\n      if s*coeff(p,x,i) < 0 then c:=c+1; s:=coef
f(p,x,i) fi;\n    fi;\n   od;\n  RETURN(c);\nend; " }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 655 "vincent:=proc(p)\n  local i,c,d,p1,p2,q,x,
\n        L,R,Z;\n  x:=indets(p)[1];\n  d:=degree(p,x);\n  q:=simplify
((y+1)^d*subs(x=1/(y+1),p));\n  c:=signchanges(q);\n  if c = 0 then RE
TURN([]) fi;\n  if c = 1 then RETURN([[0,1]]) fi;\n  p1:=simplify(2^d*
subs(x=y/2,p));\n  L:=vincent(p1);\n  for i from 1 to nops(L)\n  do\n \+
   L[i]:=[L[i][1]/2,L[i][2]/2];\n  od;\n  #p2:=simplify(2^d*subs(x=(y+
1)/2,p));\n  R:=vincent(q);\n  for i from 1 to nops(R)\n  do\n    R[i]
:=[1/(R[i][2]+1),1/(R[i][1]+1)];\n  od;\n  Z:=[];\n  for i from 1 to n
ops(R)\n  do\n    Z:=[R[i],op(Z)];\n  od;\n  if simplify(subs(x=1/2,p)
) = 0 then RETURN([op(L),[1/2,1/2],op(Z)]); fi;\n  RETURN([op(L),op(Z)
]);\n end;\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "f:=1;\nfor \+
i from 1 to 100\ndo\n  f:=f*(x-1/i):\nod:\nf:=numer(expand(f));" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "#solve(f,x);" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "Z:=vincent(f);" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 49 "for i from 1 to nops(Z) do\n  map(evalf,Z[i]
);\nod;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "map(evalf,Z[23])
;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "evalf(1/(100-23+1));" 
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "9" 0 }
{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }