I hope that's the right place to post the bug.
When solving this linear second order ODE boundary value problem:
> x^2*diff(y,x,2)-x*diff(y,x,1)+y == 0
> y(1)=1
> y(e)=2*e
Sage says Maxima can't solve it but when calling Maxima through Sage's
interface I get the right solution.
Steps to reproduce it
1. Try to solve the BVP in Sage:
y(x) = function('y')(x); y
> eq = x^2*diff(y,x,2)-x*diff(y,x,1)+y==0; eq
> h = desolve(eq, dvar=y, ivar=x,ics=[1,1,e,2*e]); h
Output:
Error in lines 3-3
Traceback (most recent call last):
File
"/projects/sage/sage-7.3/local/lib/python2.7/site-packages/smc_sagews/sage_server.py",
line 968, in execute
exec compile(block+'\n', '', 'single') in namespace, locals
File "", line 1, in <module>
File
"/projects/sage/sage-7.3/local/lib/python2.7/site-packages/sage/calculus/desolvers.py",
line 522, in desolve
raise NotImplementedError("Maxima was unable to solve this BVP. Remove
the initial condition to get the general solution.")
NotImplementedError: *Maxima was unable to solve this BVP. Remove the
initial condition to get the general solution.
*
2. Use Maxima's interface to solve the problem
maxima('depends(y,x)')
> maxima('eq: x^2*diff(y,x,2)-x*diff(y,x,1)+y=0')
> maxima('sol: ode2(eq, y,x)')
> maxima('exactSol: bc2(sol,x=1,y=1,x=%e,y=2*%e)')
Output:
...
y=x*(log(x)+1)
Tested using SageMathCloud on the 30/10/2016 so I guess Sage version 7.3
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
