Question: Should dsolve return this?

given

ode:=2*x^(1/2)*diff(y(x),x)-y(x) = -sin(x^(1/2))-cos(x^(1/2)); 
ic:=y(infinity) = y__0; 
sol:=dsolve([ode,ic]);

It gives  

This solution satisfies the ode itself. Now cos(sqrt(x)) when x=infinity is  -1..+1

But IC says y(infinity)=y0  so odetest do not verify the IC and gives this

odetest(sol,[ode,ic]);

I think dsolve should not have returned a solution at all. 

What do the experts here think of this result?

Maple 2026.1 on windows 10

Please Wait...