#9835: Make desolve more informative when solving BVP
----------------------------+-----------------------------------------------
   Reporter:  robert.marik  |       Owner:  burcin      
       Type:  defect        |      Status:  needs_review
   Priority:  minor         |   Milestone:  sage-4.6.1  
  Component:  calculus      |    Keywords:              
     Author:                |    Upstream:  N/A         
   Reviewer:                |      Merged:              
Work_issues:                |  
----------------------------+-----------------------------------------------

Comment(by robert.marik):

 Thanks, I will update the patch as time permits.

 The equation from sage-support still fails, since Maxima fails to solve
 it.  Anyway. Sage and Maxima now know that n is integer.
 {{{
 ma...@um-bc107:/opt/sage-4.5.3-Debian_Lenny_64-x86_64-Linux$ ./sage
 --maxima
 ;;; Loading
 #P"/opt/sage-4.5.3-Debian_Lenny_64-x86_64-Linux/local/lib/ecl/defsys
 tem.fas"
 ;;; Loading
 #P"/opt/sage-4.5.3-Debian_Lenny_64-x86_64-Linux/local/lib/ecl/cmp.fa s"
 ;;; Loading
 #P"/opt/sage-4.5.3-Debian_Lenny_64-x86_64-Linux/local/lib/ecl/sysfun .lsp"
 Maxima 5.20.1 http://maxima.sourceforge.net
 using Lisp ECL 10.2.1
 Distributed under the GNU Public License. See the file COPYING.
 Dedicated to the memory of William Schelter.
 The function bug_report() provides bug reporting information.
 (%i1) declare(n,integer);
 (%o1)                                done
 (%i2) eq: 'diff(x,t)^2=alpha-beta*abs(x)^n;
                           dx 2                      n
 (%o2)                    (--)  = alpha - beta abs(x)
                           dt
 (%i3) load('contrib_ode)$

 (%i4) contrib_ode(eq,x,t);
                           dx 2                      n
 (%t4)                    (--)  = alpha - beta abs(x)
                           dt

                      first order equation not linear in y'

 Improper argument to ratcoeff:
 0
 #0: linear2(expr=x,x=0)(ode2.mac line 75)
 #1: ode1_a(phi=-sqrt(alpha-beta*abs(x)^n),y=x,x=t)(ode1_lie.mac line 176)
 #2: ode1_lie(phi=-sqrt(alpha-beta*abs(x)^n),y=x,x=t)(ode1_lie.mac line 54)
  -- an error. To debug this try: debugmode(true);

 }}}

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9835#comment:4>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica, 
and MATLAB

-- 
You received this message because you are subscribed to the Google Groups 
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to 
[email protected].
For more options, visit this group at 
http://groups.google.com/group/sage-trac?hl=en.

Reply via email to