1#!/usr/bin/env python3
2
3from math import sin, pi
4
5
6def f(x):
7    return sin(x ** 2)
8
9
10def integrate_f(a, b, N):
11    s = 0
12    dx = (b - a) / N
13    for i in range(N):
14        s += f(a + i * dx)
15    return s * dx
16
17
18if __name__ == '__main__':
19    print(integrate_f(0, pi, 1000))
20