Skip main navigation

Hands-on: Numerical integration

Vectorised numerical integration using the Riemann sum.
© CSC - IT Center for Science Ltd.

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

Reach your personal and professional goals

Unlock access to hundreds of expert online courses and degrees from top universities and educators to gain accredited qualifications and professional CV-building certificates.

Join over 18 million learners to launch, switch or build upon your career, all at your own pace, across a wide range of topic areas.

Start Learning now