# -*- coding: utf-8 -*-
"""
Created on Tue Nov  8 11:25:09 2022

@author: Ylva
"""

import numpy as np

a = 0.00001 #start-x-verdi
b = np.pi #slutt-x-verdi
N = 10000 #antall x-punkter

dX = (b - a)/N #lengden mellom punktene

def sinc(x):
    return np.sin(x)/(x) #sinc kan også defineres som sin(pi*x)/(pi*x). Livet er trist

def trapes(f): #MERK: denne funksjonen tar inn en funksjon! KULT!
    area = 0
    for i in range(N):
        area += 1/2 * (f(a + dX * i) + f(a + dX * (i + 1))) * dX
        #Her går det ann å skrive x = a + dX*i og legge det inn 
        #i funksjonen for å være enda tydeligere.
    return area

print(trapes(sinc))