Opened 2 years ago

Last modified 22 months ago

#21789 new defect

Maxima can solve this bvp why not sage?

Reported by: kcrisman Owned by:
Priority: major Milestone: sage-7.5
Component: calculus Keywords:
Cc: nbruin Merged in:
Authors: Reviewers:
Report Upstream: N/A Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Description

According to this sage-devel post

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

gives an error but Maxima can nonetheless solve it.

Change History (1)

comment:1 Changed 22 months ago by mforets

strangely enough, it works if you use e^1, as in:

sage: y(x) = function('y')(x); y
x |--> y(x)
sage: eq = x^2*diff(y,x,2)-x*diff(y,x,1)+y==0; eq
x |--> x^2*diff(y(x), x, x) - x*diff(y(x), x) + y(x) == 0
sage: h = desolve(eq, dvar=y, ivar=x,ics=[1,1,e^1,2*e^1]); h
x*(log(x) + 1)
Note: See TracTickets for help on using tickets.