I've used this code for a while and never had this problem before. So I compared what I did to the ode2() documentation, and as you can see, I get the same output as the help files. But the problem I've actually been working on has a denominator of (A - 2x) and the natural domain of the problem (a draining brine tank) is [0,A/2]. But for some reason when the solution is displayed, it's got (2x - A) so the power of 3/2 is creating problems. I verified (at the bottom) that the right numbers are hidden...