# -*- coding: utf-8 -*-
"""
Spyder Editor

This is a temporary script file.
"""

import numpy as np

def f(x):
    return np.exp(- x**2) #funksjonen vår

a = 0 #start-x-verdi
b = 1 #slutt-x-verdi
n = 1000 #antall punkter på x

res = 0 #variabel som lagrer summen "res" for "resultat"

dx = (b-a)/n #Delta x
x = np.zeros(n) #Våre x-verdier
#alternativ: x = np.linspace(a,b,n) #Trenger da ikke x i for-løkke


for i in range(0,n):
    x[i] = a + i*dx #MERK: starter i a
    res += f(x[i]) #Legger til rektangelets høyde
    #Skal egentlig også gange med dx, men tar det senere
    #for å spare prosessorkraft

res *= dx #ganger med dx
print("Riemann:",res)


#Trapesmetoden på Seb-måten:
res = 0
dx = (b-a)/n #Delta x
x = np.linspace(a,b,n)
for i in range(0,n-1):
    #antar vi brukte linspace tidligere
    res += (f(x[i])+f(x[i+1]))/2
res *= dx
print("trapes:",res)