Rosetta Code

Numerical integration

Approximate definite integrals with the trapezoidal rule over fixed functions.

Edit the source, then Run. Checks start after your first Run. Ctrl/⌘ + Enter runs; Tab indents. Preview services use sample data and never contact a real service.

Sourcerosettacode/popular/numerical_integration.vibe
def square(value: float) -> float
  value * value
end

def line(value: float) -> float
  (3 * value) + 2
end

def trapezoidal_square(start_value: float, finish_value: float, steps: int) -> float
  width = (finish_value - start_value) / steps
  total = (square(start_value) + square(finish_value)) / 2
  step = 1

  while step < steps
    total = total + square(start_value + (step * width))
    step = step + 1
  end

  total * width
end

def trapezoidal_line(start_value: float, finish_value: float, steps: int) -> float
  width = (finish_value - start_value) / steps
  total = (line(start_value) + line(finish_value)) / 2
  step = 1

  while step < steps
    total = total + line(start_value + (step * width))
    step = step + 1
  end

  total * width
end

def run -> { integral_of_three_x_plus_two: float, integral_of_x_squared: float }
  {
    integral_of_x_squared: trapezoidal_square(0.0, 1.0, 50),
    integral_of_three_x_plus_two: trapezoidal_line(0.0, 5.0, 50)
  }
end

Output

Click Run example to see the result.

Code checks

Run once to load the checker. Your code stays in this browser.