Skip main navigation

Get 30% off one whole year of Unlimited learning. Subscribe for just £249.99 £174.99. T&Cs apply

Hands-on: Numerical integration

Vectorised numerical integration using the Riemann sum.
We continue studying vectorized operations with numerical integration.

Source code for this exercise is located in numpy/integration/

A simple method for evaluating integrals numerically is by the middle Riemann
sum

[S = int_a^b f(x) dx = sum_{i=1}^n f(x’_i) Delta x]

with

[x’_i = (x_i + x_{i-1}) / 2 ; x_0 = a, x_n = b]

Calculate the integral in the interval [0,π/2] and investigate how much the
Riemann sum of sin differs from 1.0. Avoid for loops. Investigate also how the results changes with the choice of Δx.

© CSC - IT Center for Science Ltd.
This article is from the free online

Python in High Performance Computing

Created by
FutureLearn - Learning For Life