# -*- coding: utf-8 -*-
"""
Created on Sun Nov 27 10:34:39 2022

@author: Sebas
"""

import numpy as np

def SND(x): #Dette er funksjonen vår, standard normalfordeling
    return 1/np.sqrt(2*np.pi) * np.exp(-(x**2)/2)

def trapesmetoden(f, a, b, N=1000): 
    #Dette er en funksjon for å ta trapesmetoden mellom a og b for en hvilken som helst funksjon "f" som gis som argument
    h = (b-a)/N #steglengde mellom a og b med N steg
    x = np.linspace(a,b,N+1) #punktene til x mellom a og b
    
    res = 0 #tellevariabel
    for i in range(N):
        res += (f(x[i])+f(x[i+1]))/2 * h #trapesmetodens definisjon
    
    return res

print(trapesmetoden(SND, 0, 1)) 
#Print standard normalfordelings sannsynlighet for et resultat mellom 0 og 1 (innenfor ett standardavvik over gjennomsnittet)
#Forventet resultat er ~34.1% (altså: 0.341...)