En orden para que todo esto sea valiosa, usted necesita para hacer una de dos cosas: 1) aumentar el número de puntos para estar seguro de que algún tipo de convergencia que está sucediendo, y/o 2) calcular el valor exacto de la integral, si es posible.
Para 2), resulta que a ser posible. Considere la posibilidad de
$$\int_0^a dx \: x \sqrt{a^3-x^3} = a^{7/2} \int_0^1 dx \: x \sqrt{1-x^3}$$
Sustituto $u=1-x^3$, $x=(1-u)^{1/3}$, $dx = -(1/3) (1-u)^{-2/3} du$:
$$\int_0^a dx\: x \sqrt{a^3-x^3} = a^{7/2} \int_0^1 du \: u^{1/2} (1-u)^{-1/3}$$
La última integral es una función beta, por lo que la integral es igual a
$$\int_0^a dx\: x \sqrt{a^3-x^3} = a^{7/2} \frac{\Gamma{\left ( \frac{3}{2}\right )}\Gamma{\left ( \frac{2}{3}\right )}}{\Gamma{\left ( \frac{13}{6}\right )}}$$
Al $a=5$, entonces la integral es $\approx 103.303$. El original de la integral es $100$ veces esta, o acerca de la $10330.3$.
Ahora que se realiza este paso, aumentar el número de puntos en su trapezoidal y de Simpson regla de aproximaciones. Están convergiendo para este resultado? Recuerde que hay algunos divertidos comportamiento de el integrando en $x=5$ (tangente vertical), por lo que incluso si aumenta, usted no puede aumentar la exactitud. Recordemos también que sólo porque hay una gran diferencia en la trapezoidal y de Simpson aproximaciones, esto no significa que usted está equivocado.
Ahora, para las aproximaciones. Vamos a considerar la regla trapezoidal para $n$ puntos de:
$$\text{trap}(n) = \frac{1}{n} \sum_{k=1}^{5 n-1} f(k/n)$$
donde $f(x)=100 x \sqrt{125-x^3}$. He hecho una tabla de valores de$1$$10$:
$$\left(
\begin{array}{cc}
1 & 9370.83 \\
2 & 9994.64 \\
3 & 10148.8 \\
4 & 10213. \\
5 & 10246.6 \\
6 & 10266.8 \\
7 & 10280. \\
8 & 10289.2 \\
9 & 10295.9 \\
10 & 10301. \\
\end{array}
\right)$$
Tenga en cuenta que, tal como hemos establecido en nuestra línea de base con el resultado exacto, nos puede sentirse bien acerca de este resultado. Su aproximación corresponde a $n=1$ y, por tanto, era un poco off, pero correcto para el nivel de aproximación buscada.
Usted puede hacer lo mismo para la regla de Simpson, excepto que usted debe dividir los pares y los impares puntos numerados:
$$\text{simp}(n) = \frac{2}{3 n} \sum_{k=1}^{(5 n-1)/2 - 1} f(2 k/n) + \frac{4}{3 n} \sum_{k=1}^{(5 n-1)/2} f((2 k-1)/n)$$
El resultado de $n=1$ $10$es
$$\left(
\begin{array}{cc}
1. & 8969.49 \\
2. & 8456.51 \\
3. & 10056.9 \\
4. & 9617.21 \\
5. & 10202.2 \\
6. & 9932.71 \\
7. & 10252.7 \\
8. & 10069. \\
9. & 10277. \\
10. & 10142. \\
\end{array}
\right)$$
Nota el comportamiento oscilatorio. De nuevo, porque hemos establecido una línea de base, sabemos lo bueno que la aproximación es. El resultado corresponde a $n=1$, y es correcto, sólo un muy bajo nivel de aproximación.