результант

Форум пользователей пакета Mathematica

Модератор: Admin

suntanka
Сообщения: 16
Зарегистрирован: Ср дек 22, 2010 9:03 pm

результант

Сообщение suntanka » Ср дек 22, 2010 9:09 pm

Подскажите, пожалуйста, в чем проблема? Может несовершенство алгоритма, а может не хватает ресурсов. Вычисляю результант для двух многочленов, первый код работает, второй нет

работает:

Код: Выделить всё

P1 = x^3 + a*x*y - 1;
P2 = y^3 + b*x*y + c*x - 1;
J = 9*x^2*y^2 + 3*a*y^3 + 3*b*x^3 - a*c*x;
Resultant[J, u1*P1 + u2*P2, y]


результат:

Код: Выделить всё

-27 a^3 u1^3 - 81 a^3 u1^2 u2 - 81 a^3 u1 u2^2 - 27 a^3 u2^3 +
 108 a^3 c u1^2 u2 x + 216 a^3 c u1 u2^2 x + 108 a^3 c u2^3 x -
 144 a^3 c^2 u1 u2^2 x^2 - 144 a^3 c^2 u2^3 x^2 -
 162 a^2 b u1^2 u2 x^3 - 324 a^2 b u1 u2^2 x^3 - 162 a^2 b u2^3 x^3 +
 64 a^3 c^3 u2^3 x^3 + 9 a^6 c u1^3 x^4 - 81 a^3 c u1^2 u2 x^4 +
 27 a^5 b c u1^2 u2 x^4 - 81 a^3 c u1 u2^2 x^4 +
 351 a^2 b c u1 u2^2 x^4 + 27 a^4 b^2 c u1 u2^2 x^4 +
 351 a^2 b c u2^3 x^4 + 9 a^3 b^3 c u2^3 x^4 +
 108 a^3 c^2 u1 u2^2 x^5 - 180 a^2 b c^2 u2^3 x^5 + 81 a^3 u1^3 x^6 -
 27 a^5 b u1^3 x^6 - 729 u1^2 u2 x^6 + 81 a^3 u1^2 u2 x^6 +
 405 a^2 b u1^2 u2 x^6 - 81 a^4 b^2 u1^2 u2 x^6 - 1458 u1 u2^2 x^6 +
 405 a^2 b u1 u2^2 x^6 - 81 a^3 b^3 u1 u2^2 x^6 - 729 u2^3 x^6 -
 27 a^2 b^4 u2^3 x^6 + 54 a^3 c u1^2 u2 x^7 + 1458 c u1 u2^2 x^7 -
 378 a^2 b c u1 u2^2 x^7 + 1458 c u2^3 x^7 - 729 c^2 u2^3 x^8 -
 54 a^3 u1^3 x^9 + 1458 u1^2 u2 x^9 - 486 a^2 b u1^2 u2 x^9 +
 1458 u1 u2^2 x^9 - 324 a b^2 u1 u2^2 x^9 - 108 b^3 u2^3 x^9 -
 1458 c u1 u2^2 x^10 - 729 u1^2 u2 x^12


второй код:

Код: Выделить всё

r03 = -162*x^3*a^2*b*u2^3 + 64*a^3*c^3*x^3*u2^3 - 27*a^3*u2^3 -
   108*b^3*x^9*u2^3 + 108*a^3*c*x*u2^3 - 729*x^8*c^2*u2^3 -
   27*x^6*b^4*a^2*u2^3 - 180*x^5*c^2*a^2*b*u2^3 - 729*x^6*u2^3 +
   1458*x^7*c*u2^3 + 351*x^4*a^2*b*c*u2^3 - 144*a^3*c^2*x^2*u2^3 +
   9*x^4*b^3*a^3*c*u2^3;
r30 = -27*a^5*x^6*b*u1^3 - 27*a^3*u1^3 + 9*a^6*x^4*c*u1^3 +
   81*x^6*a^3*u1^3 - 54*x^9*a^3*u1^3;
r21 = 54*u2*x^7*a^3*c*u1^2 + 1458*u2*x^9*u1^2 + 108*a^3*u2*c*x*u1^2 +
  81*u2*x^6*a^3*u1^2 + 405*u2*x^6*a^2*b*u1^2 - 729*u2*x^12*u1^2 -
  729*u2*x^6*u1^2 + 27*u2*x^4*a^5*c*b*u1^2 - 81*u2*x^4*a^3*c*u1^2 -
  486*a^2*x^9*u2*b*u1^2 - 81*a^3*u2*u1^2 - 162*u2*x^3*a^2*b*u1^2 -
  81*u2*x^6*b^2*a^4*u1^2; r12 = -81*u2^2*x^4*a^3*c*u1 -
  1458*u2^2*x^10*c*u1 - 1458*u2^2*x^6*u1 - 324*u2^2*x^9*b^2*a*u1 +
  216*a^3*u2^2*c*x*u1 - 378*u2^2*x^7*a^2*c*b*u1 +
  27*u2^2*x^4*b^2*a^4*c*u1 + 405*u2^2*x^6*a^2*b*u1 +
  351*u2^2*x^4*a^2*c*b*u1 - 144*a^3*u2^2*c^2*x^2*u1 -
  81*a^3*u2^2*u1 + 1458*u2^2*x^9*u1 - 81*a^3*u2^2*b^3*x^6*u1 -
  324*u2^2*x^3*a^2*b*u1 + 108*u2^2*x^5*c^2*a^3*u1 + 1458*u2^2*x^7*c*u1;
Resultant[r30, v1*r03 + v2*r21 + v3*r12, x]


результат нет :(

TrushkovVV
Сообщения: 45
Зарегистрирован: Сб ноя 14, 2009 2:55 pm
Контактная информация:

Re: результант

Сообщение TrushkovVV » Чт дек 23, 2010 7:42 am

Умножение - это не звездочка, а пробел.