Post new topic Reply to topic  [ 3 posts ] 
Author Message
 Post subject: Resultant
PostPosted: Thu Mar 04, 2010 9:10 pm 
I want to compute Resultant of two polynomial:

ring r = (0,t), (a(1..3),b(1..3),c(1..3),d(1..3),s) , lp;
poly x = s^3 + a(3)*s^2 + a(2)*s + a(3);
poly y = s^3 + b(3)*s^2 + b(2)*s + b(3);
poly z = s^3 + c(3)*s^2 + c(2)*s + c(3);
poly u = s^3 + d(3)*s^2 + d(2)*s + d(3);
poly f = x^2 - s*u^2 + s - (s^2*u^2-s)*y^2;
poly g = s*x^2 - 2*s^2*u^2 + 1 - s*(s^2*u^2 - s)*z^2;
resultant (f,g,s);

But it seems that Singular does not give me the answer, since there comes nothing. Can you help me please?
--------------
By the way i have a question, as above, but i want to have a(1..n), b(1..n), c(1..n), d(1..n), so i have to define

int n;

But then i have trouble with writting down polynomials, e.g. i want to define

x = s^n + \sum_{i=1}^{n} a_i * s^{n-i}

How can i represent it in SINGULAR?


Report this post
Top
  
Reply with quote  
 Post subject: Re: Resultant
PostPosted: Sat Mar 06, 2010 4:32 am 

Joined: Mon Aug 29, 2005 9:22 am
Posts: 41
Location: Kaiserslautern, Germany
The resultant computation is perhaps too long. To see that Singular indeed computes something, and to watch the progress of a computation, you should type option(prot); before the computation starts.

To your second question:
int n = 37; (n must be here a concete integer)
ring r = 0,(a(1..n),s),dp;
int ii;
poly x = s^n;
for( ii=n; ii>0; ii--)
{
x = x + (-1)^n*a(ii)*s^(n-ii);
}

GMG


Report this post
Top
 Profile  
Reply with quote  
 Post subject: Re: Resultant
PostPosted: Sat Mar 06, 2010 4:37 am 

Joined: Mon Aug 29, 2005 9:22 am
Posts: 41
Location: Kaiserslautern, Germany
Sorry, option(prot); does not show anything for resultant.
You can use option(mem); to see the memory usage.


Report this post
Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 3 posts ] 

You can post new topics in this forum
You can reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

It is currently Fri May 13, 2022 11:06 am
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group