Maxima can solve this bvp why not sage?

[kcrisman]

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.

CC: @nbruin

Component: calculus

Issue created by migration from https://trac.sagemath.org/ticket/21789

[c22b6800-ec0b-4cbf-94c4-0a74aecc2093]

comment:1

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)

[7822f248-ba56-45f1-ab3d-4de7482bdf9f]

comment:2

May be related to #24102 (see this comment...)