from math import exp, pi
from numpy import linspace
from time import time

def f(x, mu=0, sigma=1):
    return exp(-(x-mu)**2/(2*sigma**2))/((2*pi*sigma**2)**(1/2))

def midpoint(x, h):
    # TODO: Skriv ferdig funksjonen
    return 0

def integral(xmin, xmax, n):
    t0 = time()
    print()
    print("n = ", n)
    xs = linspace(xmin, xmax, n)
    h = (xmax - xmin)/(n - 1)
    print("h = ", h)
    sum = 0.0
    for x in xs[:-1]:
        sum += midpoint(x, h)
    t1 = time()
    t = t1 - t0
    return sum, t

n = 10000
numcalc, t = integral(-10, 10, n)
print("Numerical calculation:", numcalc)
print("Time (seconds):", round(t, 4))