JavaScript/Maths/MidpointIntegration.js at master · TheAlgorithms/JavaScript · GitHub

History

63 lines (54 loc) · 1.88 KB

Raw

/**

* @title Midpoint rule for definite integral evaluation

* @author [ggkogkou](https://github.com/ggkogkou)

* @brief Calculate definite integrals with midpoint method

* @details The idea is to split the interval in a number N of intervals and use as interpolation points the xi

* for which it applies that xi = x0 + i*h, where h is a step defined as h = (b-a)/N where a and b are the

* first and last points of the interval of the integration [a, b].

* We create a table of the xi and their corresponding f(xi) values and we evaluate the integral by the formula:

* I = h * {f(x0+h/2) + f(x1+h/2) + ... + f(xN-1+h/2)}

* N must be > 0 and a<b. By increasing N, we also increase precision

* [More info link](https://tutorial.math.lamar.edu/classes/calcii/approximatingdefintegrals.aspx)

function integralEvaluation(N, a, b, func) {

// Check if all restrictions are satisfied for the given N, a, b

if (!Number.isInteger(N) || Number.isNaN(a) || Number.isNaN(b)) {

throw new TypeError('Expected integer N and finite a, b')

}

if (N <= 0) {

throw Error('N has to be >= 2')

} // check if N > 0

if (a > b) {

throw Error('a must be less or equal than b')

} // Check if a < b

if (a === b) return 0 // If a === b integral is zero

// Calculate the step h

const h = (b - a) / N

// Find interpolation points

let xi = a // initialize xi = x0

const pointsArray = []

// Find the sum {f(x0+h/2) + f(x1+h/2) + ... + f(xN-1+h/2)}

let temp

for (let i = 0; i < N; i++) {

temp = func(xi + h / 2)

pointsArray.push(temp)

xi += h

}

// Calculate the integral

let result = h

temp = pointsArray.reduce((acc, currValue) => acc + currValue, 0)

result *= temp

if (Number.isNaN(result)) {

throw Error(

'Result is NaN. The input interval does not belong to the functions domain'

)

}

return result

}

export { integralEvaluation }

JavaScript/Maths/MidpointIntegration.js at master · TheAlgorithms/JavaScript · GitHub

Uh oh!

Uh oh!

JavaScript/Maths/MidpointIntegration.js at master · TheAlgorithms/JavaScript · GitHub

Uh oh!

Uh oh!

Expand file tree

Breadcrumbs

MidpointIntegration.js

Latest commit

History

Breadcrumbs

MidpointIntegration.js

File metadata and controls