[Numpy-discussion] roots and high-order polynomial

Fabrice Silva silva@lma.cnrs-mrs...
Fri Jul 3 09:54:29 CDT 2009


Le vendredi 03 juillet 2009 à 14:43 +0200, Nils Wagner a écrit :
> Just curious - Can you provide us with the coefficients of 
> your polynomial ?



Working case :
Polynomial.c = 
[ -1.34100085e+57 +0.00000000e+00j  -2.28806781e+55 +0.00000000e+00j
  -4.34808480e+54 -3.27208577e+36j  -2.44499178e+52 -7.45966050e+34j
  -4.84319247e+51 -1.13783265e+34j   1.84249681e+47 +1.01074670e+33j
  -2.27761972e+48 +2.06997398e+31j   5.35043855e+45 +1.50239041e+30j
  -4.40346465e+44 +3.49745869e+27j   1.16055383e+42 -6.38037062e+25j
  -2.79370528e+40 -1.60550605e+23j   4.72851679e+37 +7.68429228e+20j
  -2.51340384e+35 +2.92105561e+17j]

A non-working case : pretty large range of values!
Polynomial.c = 
[ -1.34100085e+219 +0.00000000e+000j  -3.87620512e+217 +0.00000000e+000j
  -1.59644076e+218 -4.96044532e+200j  -4.36379116e+216 -1.78047353e+201j
  -8.66515071e+216 -7.57270655e+199j  -2.22899258e+215 +1.20695742e+200j
  -2.84415205e+215 +2.04437917e+198j  -6.84491536e+213 +4.11962630e+198j
  -6.31415522e+213 +8.60411394e+196j  -1.41164416e+212 +6.40341637e+195j
  -1.00484766e+212 +1.05851522e+194j  -2.06870108e+210 -6.90610706e+194j
  -1.18557253e+210 -7.04790806e+192j  -2.22304125e+208 +8.69684716e+191j
  -1.05816505e+208 +4.10351383e+190j  -1.78206407e+206 -3.01765476e+190j
  -7.22785804e+205 -3.75229199e+188j  -1.07357199e+204 +1.60517046e+188j
  -3.79912531e+203 +5.72547583e+186j  -4.85746209e+201 -6.62989846e+185j
  -1.53768394e+201 +6.29641862e+183j  -1.63653463e+199 -3.94069765e+183j
  -4.77535548e+198 +2.70460428e+180j  -4.02975379e+196 -8.22561006e+180j
  -1.12909884e+196 -3.28170925e+178j  -7.00288278e+193 +1.87782740e+178j
  -2.00766909e+193 +5.47118767e+176j  -7.99891319e+190 +5.56488189e+175j
  -2.63703112e+190 +1.09163655e+174j  -4.91486942e+187 +6.37859777e+172j
  -2.49451367e+187 +1.07180124e+171j   1.17174241e+183 +4.42130463e+169j
  -1.63881690e+184 -1.39641308e+167j   2.63207649e+181 -5.52568751e+165j
  -7.08593327e+180 +8.81583002e+163j   1.88293454e+178 +2.44846832e+162j
  -1.85373357e+177 +1.80169585e+160j   5.67531810e+174 +8.34051533e+158j
  -2.54119737e+173 +6.83930275e+156j   6.95713759e+170 -3.36080695e+154j
  -1.36763548e+169 -4.25101484e+151j   2.31484033e+166 +2.65804116e+149j
  -1.19894847e+164 +1.50542119e+146j]

-- 
Fabrice Silva <silva@lma.cnrs-mrs.fr>
LMA UPR CNRS 7051



More information about the NumPy-Discussion mailing list