#1302 Another solve() problem

None
closed
5
2016-10-02
2007-12-17
Jan Fricke
No

Maxima 5.13.0 http://maxima.sourceforge.net
Using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (aka GCL)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
This is a development version of Maxima. The function bug_report()
provides bug reporting information.
(%i1) solve([x^2+y^2=2^2,(x-1)^2+(y-1)^2=2^2]);
(%o1) []
(%i2)

Related

Bugs: #3208

Discussion

  • Dieter Kaiser

    Dieter Kaiser - 2010-10-12

    Copying the example from the bug report ID: 2037060 - Symbolic Center of Circle Calculation into this thread. It is the same problem:

    -------------------------------------------
    (%i1) s: solve([(x-x0)^2+(y-y0)^2=r^2,(x-x1)^2+(y-y1)^2=r^2],[x,y]);
    (%o1) []
    (%i2)

    Maxima uses up a lot of processing power and doesn't find the correct answer.

    ftp://ftp.deltacompsys.com/public/NG/Mathcad%20-%20CenterOfCircleCalculations.pdf

    Peter Nachtwey
    peter@deltamotion.com
    -------------------------------------------

    Dieter Kaiser

     
  • David Billinghurst

    • status: open --> closed
    • assigned_to: David Billinghurst
    • Group: --> None
     
  • David Billinghurst

    Fixed by commit [2986fc] See [bugs:#3208] for details

    (%i1) solve([x^2+y^2=2^2,(x-1)^2+(y-1)^2=2^2]);
    
                                            sqrt(7)      1
                                   sqrt(2) (------- - -------) + 2
                  sqrt(7) - 1               sqrt(2)   sqrt(2)
    (%o1) [[y = - -----------, x = -------------------------------],
                       2                          2
                                                           sqrt(7)      1
                                                  sqrt(2) (------- + -------) - 2
                               sqrt(7) + 1                 sqrt(2)   sqrt(2)
                          [y = -----------, x = - -------------------------------]]
                                    2                            2
    (%i2) time(%o1);
    
    (%o2)                               [0.016]
    
     

    Related

    Bugs: #3208
    Commit: [2986fc]


    Last edit: David Billinghurst 2016-10-02

Log in to post a comment.

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:





No, thanks