maple 16 无法与内核连接

ttttmmmm

用maple 16 在win7下 求高次多项式零点 当最高次数较小时可以计算 怎加最高次数 就报错无法与内核连接
  求解~~~~
另外，maple 始终没有调用多核去运算，需要另外设置吗

feiyuzhen

我没有用win7，一些建议如下：
1.版本问题，maple16的版本是32还是64位的，win7的版本是否一致？
2.多项式零点的话应该没问题，符号解低于5次是可以的，高的话maple的结果是rootof提示符，数值解更不是问题

ttttmmmm

求助大神：在maple中使用fsolve能瞬间求出方程一中多项式的零点，所以应该不是问题1，而当稍微增加最高项次数，到方程二中，就完全无法计算出结果，感觉应该不是计算量的问题。
方程一：
-7.506556000*10^(-23)*omega^12+2.243366854*10^(-78)*omega^34+3.651941282*10^(-50)*omega^25+1.145341958*10^(-105)*omega^41-3.837716260*10^(-102)*omega^40+9.952750000*10^(-10)*omega^5+2.520327960*10^(-113)*omega^43-5.961864252*10^(-98)*omega^39+2.934254874*10^(-60)*omega^29+3.877088344*10^(-109)*omega^42-3.482296310*10^(-82)*omega^35+2.846747460*10^(-123)*omega^45-1.089582984*10^(-70)*omega^32-5.927854900*10^(-48)*omega^24+0.7371682860e-5*omega^4+1.150636782*10^(-52)*omega^26+1.000000000*10^(-174)*omega^56-3.421709139*10^(-86)*omega^36+3.115541000*10^(-149)*omega^51+0.4898830000e-5*omega^3-1.137646660*10^(-117)*omega^44+4.672450454*10^(-94)*omega^38+1.000000000*10^(-171)*omega^55+4.857583679*10^(-90)*omega^37-2.199194507*10^(-150)*omega^52-1.000000000*10^(-188)*omega^59-1.000000000*10^(-191)*omega^60+9.074004810*10^(-126)*omega^46+5.000000000*10^(-181)*omega^57+1.032629580*10^(-14)*omega^8+4.638000000*10^(-200)*omega^62-6.192226965*10^(-55)*omega^27-3.370200000*10^(-44)*omega^21-1.173960000*10^(-36)*omega^18+7.406933000*10^(-44)*omega^22+1.220000000*10^(-183)*omega^58+6.646947000*10^(-140)*omega^49-4.374587100*10^(-46)*omega^23-8.234527370*10^(-58)*omega^28-1.861180000*10^(-30)*omega^15-6.739449245*10^(-134)*omega^48+4.485566928*10^(-10)*omega^6-1.628871800*10^(-131)*omega^47-4.058700000*10^(-41)*omega^20-2.179900000*10^(-39)*omega^19+58.56452040-1.753122200*10^(-20)*omega^11+8.286312400*10^(-19)*omega^10-1.000000000*10^(-208)*omega^64+3.319370000*10^(-27)*omega^14+0.4456662984e-1*omega^2-7.350510339*10^(-67)*omega^31-1.559547000*10^(-156)*omega^53+0.6506950000e-2*omega+3.986744000*10^(-25)*omega^13+2.630858102*10^(-63)*omega^30-5.244560000*10^(-14)*omega^7+5.468958700*10^(-159)*omega^54+3.894350000*10^(-34)*omega^17+2.997800000*10^(-32)*omega^16+4.438707600*10^(-142)*omega^50+1.200000000*10^(-205)*omega^63+2.099570412*10^(-74)*omega^33+1.368262600*10^(-16)*omega^9 = 0

方程二
-5.328930497*10^(-142)*omega^53+1.237165060*10^(-121)*omega^48+3.132963202*10^(-134)*omega^51-8.845398500*10^(-114)*omega^46-58.56451500-1.585807425*10^(-129)*omega^50+6.284390000*10^(-67)*omega^30-4.354200000*10^(-64)*omega^29-9.063300000*10^(-28)*omega^13-0.6505470000e-2*omega-7.113926804*10^(-153)*omega^56+1.967398000*10^(-31)*omega^15+2.223385000*10^(-29)*omega^14-3.005981008*10^(-145)*omega^54-1.602789610*10^(-149)*omega^55-2.670972550*10^(-126)*omega^49-2.062702871*10^(-156)*omega^57-7.290678396*10^(-160)*omega^58+1.017808700*10^(-37)*omega^18-6.599905000*10^(-40)*omega^19+3.161508000*10^(-42)*omega^20-9.774340000*10^(-45)*omega^21-2.920000000*10^(-239)*omega^75-8.800000000*10^(-235)*omega^74-1.608093182*10^(-281)*omega^85-2.000000000*10^(-277)*omega^84-6.300000000*10^(-269)*omega^82-9.680000000*10^(-273)*omega^83-1.249110000*10^(-49)*omega^23+3.707220000*10^(-47)*omega^22+5.600000000*10^(-243)*omega^76+3.227228700*10^(-166)*omega^59-4.977817000*10^(-207)*omega^69-2.671297000*10^(-174)*omega^61+1.779049231*10^(-169)*omega^60+1.231013995*10^(-215)*omega^71-2.703286505*10^(-210)*omega^70+5.559268853*10^(-202)*omega^68+1.010944300*10^(-198)*omega^67-1.475405321*10^(-177)*omega^62+1.165031374*10^(-185)*omega^64+2.106651400*10^(-182)*omega^63-0.4949890000e-5*omega^3-0.7370082900e-5*omega^4-1.548250097*10^(-190)*omega^65-0.4456682430e-1*omega^2+2.426000000*10^(-247)*omega^77+6.597038000*10^(-219)*omega^72-8.561061300*10^(-194)*omega^66-1.000000000*10^(-255)*omega^79-8.184670000*10^(-10)*omega^5+5.378000000*10^(-264)*omega^81-1.430850000*10^(-18)*omega^9-1.374700620*10^(-14)*omega^8-5.031350000*10^(-14)*omega^7-4.514233800*10^(-10)*omega^6-2.418953400*10^(-19)*omega^10-2.566360000*10^(-23)*omega^11-1.686868900*10^(-35)*omega^17-1.311732500*10^(-33)*omega^16+8.019500000*10^(-63)*omega^28+2.542000000*10^(-60)*omega^27-3.140226000*10^(-24)*omega^12-1.166250000*10^(-54)*omega^25+3.979610000*10^(-52)*omega^24-1.922907520*10^(-98)*omega^42-1.155606000*10^(-68)*omega^31+2.354260000*10^(-57)*omega^26+2.068130150*10^(-118)*omega^47-1.462791588*10^(-110)*omega^45+5.331826784*10^(-106)*omega^44+8.647300300*10^(-103)*omega^43+5.961174800*10^(-89)*omega^39-1.133247470*10^(-85)*omega^38-2.471905100*10^(-83)*omega^37-2.903800750*10^(-95)*omega^41+1.353674471*10^(-91)*omega^40+3.839986430*10^(-80)*omega^36+2.777638500*10^(-78)*omega^35-5.151661900*10^(-75)*omega^34+1.047186000*10^(-73)*omega^33+2.039300800*10^(-70)*omega^32+1.855519426*10^(-137)*omega^52 = 0

feiyuzhen

用fsolve计算的方程一结果如下：
    -94114.207656771, -15380.844049232, -9216.9135698340,

      -942.60732312467, -397.06220442507, -272.66293589432,

      -148.91774928318, 4246.4437268739, 5565.6653185140,

      94307.066789625，
方程二的结果如下：
    -18038.891523284, -3728.3576262443, -443.57025649626,

      -313.51529859293, -193.57490351760, 366.42808282128,

      2220.1018753893, 18151.886658634, 63712.287490560
因为omega在复数域的根太多，时间太长，上述计算只是在实数域上进行。

计算环境：xp+maple14

多核计算需要用maple[Threads]编程，我没有用过

btw，如果已知该多项式为某矩阵的特征多项式的话，用求特征值的方法获取omega的值是很快的。