Re: [Maxima-discuss] solve() involving sqrt()

[Raymond T. <toy...@gm...>]

On Mon, Sep 14, 2020 at 4:31 AM Me Self <vla...@gm...> wrote:

> Hello Everyone
>
>
> This system of polynomials involves sqrt():
>
> eq1 : 0.97561*sqrt(x*y)*z-0.009879$
> eq2 : 0.097561*z-0.4878*sqrt(x*y)+0.2439*y+0.2439*x-0.98506$
> eq3 : 0.02439*y*z+0.97561*x*z-0.90247$
>
> But solve() seems to fail, to_poly_solve() takes ages without results,
> while minpack_solve aborts with:
>
>
What starting points did you use?

I got:

minpack_solve([eq1,eq2,eq3],[x,y,z], [1,1,1]);
Evaluation took 0.0100 seconds (0.0100 elapsed) using 942.305 KB.

[[5.740409864561979,-1.755503737103853,0.1624466869294288],
 0.01024886625263842,4]$

>
> I think it's because there's both sqrt(x*y) and x+y in the 3rd
> equation. Is there any workaround for this? The hompack version that
> Raymond Toy is working on seems to accept non-integer powers for the
> coefficients, but it's not in Maxima (yet?).
>

No only non-negative integer powers are allowed.  At least for polsys.

>
> Regards,
> Vlad
>
>
> _______________________________________________
> Maxima-discuss mailing list
> Max...@li...
> https://lists.sourceforge.net/lists/listinfo/maxima-discuss
>


-- 
Ray