¿Cómo calcular integrales usando el método de Simpson?

Daré un ejemplo de integración, como (θ^3*(1-θ)^4)dθ

El código es el siguiente:

Opción explícita

Private Sub form_Click()

Atenuar A como entero, B como entero, paso como entero, simple como simple

A = 0

B = 1

Paso = 10000

Simp = simps(A, B. Paso)

Sub formulario privado_Click()

Sub formulario privado_Click () Paso)

Imprimir Simp

End Sub

Función privada Diversión(X como única) Como única

Diversión = X ^ 3 * (1 - X) ^ 4

Función final

Función privada simps(A como entero, B como entero, paso como entero) como simple

Dim H como único, cuatro como único, S como único, X como único

Dim H como único, cuatro como único, dos como único, S como único, X como único, I como un número entero<. /p>

H = (B - A) / Paso

X = A

S = Diversión((A)) + Diversión((B))

Para I = 1 Al Paso / 2

X = X + H

Cuatro = Cuatro + Diversión(X)

X = X + H

Dos = Dos + Diversión(X)

Siguiente I

Dos = Dos - Diversión(X)