fsolve.m: Fix undefined output error when using Output function (bug #60144).
Maxima is returning the answer for a different integral. E.g. in MATLAB: (1-exp(-pi/2))/2 ans = 0.3961 quad(@(x) exp(-acos(x)),0,1) ans = 0.3961 Another way to look at it is that Maxima has chosen a different branch. Ordinarily, one would say that 0 < acos(x) < pi/2 if 0 < x < 1. That's the branch that MATLAB uses: acos(.9) ans = 0.4510 cos(.4510) ans = 0.9000 However, cos(x) is of course an even function: cos(-.4510) ans = 0.9000 ...which means one could claim that acos(.9) = -.4510 That's essentially...
I would like to second the original bug report. I've recently been dealing with a number of nested radicals, and was disappointed that Maxima couldn't simplify them. I wrote an awkward workaround. I certainly didn't find that loadable package. And I don't consider that a good solution anyway. Half the reason for using Maxima is automatic simplification. I should not have to rummage through the toolbox looking for something appropriate - or to rely on the kindness of strangers to point the way. I'll...