These two issues probably came up before, but I can't find where and when searching Maple primes.
So I thought to summarize the issues I see with odetest in one post, in the hope to get clarification on current status on these from the powers who know.
The first issue
The order in which odetest returns the answer. When odetest is called to check the ode and IC, as in
the_status := odetest(sol,[ode,IC])
One would expect the_status to be a list, where the first entry tells if sol verifies the ode, and the second entry tells if sol verifies IC.
i.e. the order is the same as in the input. right? Since ode is first and IC is second in the input list.
But Maple sometimes mixes the order. See example 1 below. This makes it impossible to determine if the solution verifies the ode or IC, when one of the entries in the_status is zero and the other is not, since order can be reversed sometimes.
Second issue:
When the solution is implicit, Maple gives invalid odetest result on the IC, unless one rewrites the solution using (lhs-rhs)(sol)=0.
i.e. move everything to the left side of the equation with RHS zero. This happens sometimes and when the solution is implicit.
I have thought this was fixed in current Maple, but it is not. I remember this came up before, but can't find when and where.
Example 2 below shows an example.
Will these two issues hopefully be fixed in Maple 2025? Sometimes one forgets to rewrite the solution using (lhs-rhs)(sol)=0 and this results in false negative.
Please see worksheet below. ps. I hope forum manager does not delete this question.
Download issues_with_odetest_dec_16_2024.mw