Robert Dodier <robert.dodier <at> gmail.com> writes:
>
> For the example at hand, assuming that the x(... is supposed to be x*(...
> I get the following.
>
> f(x) := (500*((0.0075 + log((x/15.5) + 1))^0.2)*x*((5.5^2) -
> (x^2))^0.5 - ((4.5^2) - (x^2))^0.5);
> quad_qags (f(x), x, 0, 4.5);
> => [15624.57336771343, 6.9011136083018057E-5, 273, 0]
>
With quad_qags I get the correct result, romberg did not match, thank you.
Copying the formula in two pieces I lost a bracket-pair (and the multiplier),
excuse me...
thx
Hans W. Hofmann