I found a bug in Maple V-4: when you use dsolve with initial conditions, and when there are two or more values of the constants that can be found, Maple gives just one solution. Example
But there are two solutions.
The workaround for this bug is to patch the end of `dsolve/consts`:
But to use this patch, you must unprotect dsolve.
The bug is removed with Maple V Release 5. (U. Klein)
There is indeed a bug in dsolve (R4), but the problem is not that only one solution is reported for the initial value problem given by Alain Leroux. The problem is that the wrong solution is reported!
The initial value problem has a unique solution (locally, - by the general theorem on uniqueness).
Consider the following two simpliﬁed versions:
The last solution is obviously incorrect, as equ2 is equivalent to equ1. That the correct solution can be joined with the constant solution \(x=1\) at \(t=2\) is irrelevant in this context.
Try this one
\[ y(x)''+2/x y(x)' - m^2 y(x)=0 \]
the solution in R4 is mixed exponentials and hyperbolics. Go ﬁgure! Why not just return hyperbolics or exponentials as R2 and R3 did.