runjaj, Spanish Calcular una integral numérica usando la regla de Simpson con #julialang es tan fácil como usar esta función:
∫(x, y) = (x[2]-x[1])/3*(y[1] + 4*sum(y[2:2:end-1]) + 2*sum(y[3:2:end-2]) + y[end])
Por ejemplo, la integral del seno entre 0 y π es:
x = LinRange(0, π, 7) # 7 puntos entre 0 y π y = sin.(x) # Cálculo del seno para los 7 puntos ∫(x, y) # Integración 2.000863189673536
El valor realizando la integral simbólicamente es 2. Todo bien sencillo.
Add comment