Re: Is it possible to have several solutions in this wa y to this equation using Pari/ɢᴘ ?

Laël Cellier <[email protected]>
Newsgroups gmane.comp.mathematics.pari.user
Message-ID <[email protected]>
So yes, I want to solve a system of equations where alpha and beta can 
be set to any values as long as they’re valid for solving the first 
equation.

Are you meaning this/your code allows to get more than 2 solutions for 
the second equation ?

Cordialement,

Le 19/01/2025 à 21:39, Bill Allombert a écrit :
> On Sun, Jan 19, 2025 at 08:05:09PM +0100, Laël Cellier wrote:
>> first,
>> beta=-(V\W);
>> alpha=W*(V+W*beta);
>>
>> is just 1 way to find a suitable solution alpha == w (v + w beta). I’m
>> needing to find other ways in order to get different results.
>>
>> Let’s give a numerical example :
>> V=25 c=60 W=3 b=85 f=-1
>> give :
>> alpha=3
>> beta=-8
>> nfr=[-4, -1/9]~
> So you want this system of equation to have solutions:
>
> alpha^2*x^2+(2*alpha*beta-f*b)*x+(beta^2-c) = 0
> alpha=W*(V+W*beta);
>
> you can eliminate alpha and obtain
>
> Q = subst(P,alpha,W*(V+W*beta))
> Q=(W^4*beta^2+2*V*W^3*beta+V^2*W^2)*x^2+(2*W^2*beta^2+2*V*W*beta-b*f)*x+(beta^2-c)
> To have rational solution the discriminant of Q needs to be a square
> ? poldisc(Q)
> %3 = (-4*W^2*b*f+4*W^4*c)*beta^2+(-4*V*W*b*f+8*V*W^3*c)*beta+(b^2*f^2+4*V^2*W^2*c)
>
> so we need to solve for (beta, D)
>
> D^2= (-4*W^2*b*f+4*W^4*c)*beta^2+(-4*V*W*b*f+8*V*W^3*c)*beta+(b^2*f^2+4*V^2*W^2*c)
>
> This is in general a conic which can be solved with qfsolve/qfparam
> (or hyperellratpoints).
>
> fun(V,W,b,c,f)=
> {
>    my(A= (-4*W^2*b*f+4*W^4*c), B = (-4*V*W*b*f+8*V*W^3*c)/2, C = b^2*f^2+4*V^2*W^2*c);
>    my(M=[A, B, 0; B, C, 0; 0, 0, -1]);
>    my(S=qfsolve(M));
>    my(beta =S[1]/S[2]);
>    my(alpha= W*(V+W*beta));
>    my(X=nfroots(,alpha^2*x^2+(2*alpha*beta-f*b)*x+(beta^2-c)));
>    [alpha,beta,X];
> }
>
> However in your example " V=25 c=60 W=3 b=85 f=-1 ",
> the conic degenerates to
> D^2 = 25*(30*beta+233)^2
> so you can pick whatever value you want for beta...
> Maybe this is the actual mistery ?
>
> Cheers,
> Bill.
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.