Singular :: View topic - polynomial solving with Singular or Maple

> LIB "primdec.lib";
> ring r=0,(s1,s12,s13,s123,lambda),lp;
> ideal i = 1402*s13*s1+9*s1^2-20*s12*s1*s13+66*s12*s1*s13^2+66*s12*s1^2*s13
. -2822*s1^3+5608*s1^4-8430*s1^2*s13-5608*s13^2*s1+16824*s1^3*s13
. +16824*s13^2*s1^2+5608*s1*s13^3-9*s12*s1^2-11*s12*s13^2+22*s12*s1^3
. +22*s12*s13^3,
. 1691*s12*s1+1691*s12*s13-300*s1+600*s1^2+600*s13*s1,
. -1091*s13*s1+302*s1^2-602*s12*s1*s13+1812*s12*s1*s13^2+1812*s12*s1^2*s13
. -1208*s1^3+1208*s1^4+370*s1^2*s13+1578*s13^2*s1+3624*s1^3*s13
. +3624*s13^2*s1^2+1208*s1*s13^3-302*s12*s1^2-300*s12*s13^2+604*s12*s1^3
. +604*s12*s13^3;
> primdecGTZ(i);
[1]:
[1]:
_[1]=s13^2
_[2]=s1*s13
_[3]=1691*s1*s12-300*s1+1691*s12*s13
_[4]=s1^2
[2]:
_[1]=s13
_[2]=s1
[2]:
[1]:
_[1]=s12
_[2]=2*s1+2*s13-1
[2]:
_[1]=s12
_[2]=2*s1+2*s13-1
[3]:
[1]:
_[1]=s12
_[2]=s1
[2]:
_[1]=s12
_[2]=s1
[4]:
[1]:
_[1]=4016016*s13-46961
_[2]=282961794*s12-36018575
_[3]=4016016*s1-264661
[2]:
_[1]=4016016*s13-46961
_[2]=282961794*s12-36018575
_[3]=4016016*s1-264661

Author:	Guest [ Thu Aug 11, 2005 5:31 pm ]
Post subject:	polynomial solving with Singular or Maple
I have a system of polynomial equations, when I let Singular solve it it produces the message "Ideal i has to be 0-dimentional". Somehow, Singular does not allow to solve it, while Maple does allow. Do you know why? email: [email protected] Posted in old Singular Forum on: 2002-02-03 10:01:23+01

Author:	greuel [ Thu Sep 08, 2005 11:41 am ]
Post subject:
Singular can solve numerically only 0-dim systems by using the library solve.lib. This is just redisigned and will work much better in near future. Symbolically, Singular can treat systems of any dimension. For that you should consult the library primdec.lib. So, for the first system you can use primdecGTZ, it gives all irreducible components ([2]:) together with the primary components ([1]:). Setting the polynomials to 0 gives the solutions (as you see, the solution set is 1-dimensional): Code: > LIB "primdec.lib"; > ring r=0,(s1,s12,s13,s123,lambda),lp; > ideal i = 1402s13s1+9s1^2-20s12s1s13+66s12s1s13^2+66s12s1^2s13 . -2822s1^3+5608s1^4-8430s1^2s13-5608s13^2s1+16824s1^3s13 . +16824s13^2s1^2+5608s1s13^3-9s12s1^2-11s12s13^2+22s12s1^3 . +22s12s13^3, . 1691s12s1+1691s12s13-300s1+600s1^2+600s13s1, . -1091s13s1+302s1^2-602s12s1s13+1812s12s1s13^2+1812s12s1^2s13 . -1208s1^3+1208s1^4+370s1^2s13+1578s13^2s1+3624s1^3s13 . +3624s13^2s1^2+1208s1s13^3-302s12s1^2-300s12s13^2+604s12s1^3 . +604s12s13^3; > primdecGTZ(i); [1]: [1]: _[1]=s13^2 _[2]=s1s13 _[3]=1691s1s12-300s1+1691s12s13 _[4]=s1^2 [2]: _[1]=s13 _[2]=s1 [2]: [1]: _[1]=s12 _[2]=2s1+2s13-1 [2]: _[1]=s12 _[2]=2s1+2s13-1 [3]: [1]: _[1]=s12 _[2]=s1 [2]: _[1]=s12 _[2]=s1 [4]: [1]: _[1]=4016016s13-46961 _[2]=282961794s12-36018575 _[3]=4016016s1-264661 [2]: _[1]=4016016s13-46961 _[2]=282961794s12-36018575 _[3]=4016016s1-264661 email: [email protected] Posted in old Singular Forum on: 2002-02-03 10:03:27+01

Singular https://www.singular.uni-kl.de/forum/

polynomial solving with Singular or Maple https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=1314	Page 1 of 1

Page 1 of 1	All times are UTC + 1 hour [ DST ]
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group http://www.phpbb.com/