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