Join Our Mystery Crochet Along!
def trapezoidal(f, a, b, n): h = (b - a) / n s = 0.5 * (f(a) + f(b)) for i in range(1, n): s += f(a + i*h) return h * s
# Exemple: int_0^1 e^{-x^2} dx f = lambda x: math.exp(-x*x) approx = trapezoidal(f, 0.0, 1.0, 1000) print("Approx:", approx) function T = trapezoidal(f, a, b, n) h = (b - a) / n; x = a:h:b; y = f(x); T = h*(0.5*y(1) + sum(y(2:end-1)) + 0.5*y(end)); end def trapezoidal(f, a, b, n): h = (b - a) / n s = 0
Join Our Mystery Crochet Along!
Join our free crochet along and unravel the mystery as you stitch along with other fellow crocheters.
Thank you for signing up to our crochet along.