?# Program som tester tre ulike numeriske derivasjonsmetoder

from math import *

def dsin1(x,h):
	return (sin(x+h)-sin(x))/h # Enklest mulige
	
def dsin2(x,h):
	return (sin(x+h)-sin(x-h))/(2*h) # Basert p? interpolasjon med parabel

def dsin3(x,h):
	return (sin(x-2*h)-8*sin(x-h)+8*sin(x+h)-sin(x+2*h))/(12*h)
	     # Basert p? interpolasjon med 4. grads polynomer
	
for i in range(0,20):
	print ('i: %2i    e1: % 3.6e     e2: % 3.6e    e3: % 3.6e'
	%(i, dsin1(1,10**(-i))-cos(1),
	     dsin2(1,10**(-i))-cos(1),
	     dsin3(1,10**(-i))-cos(1)))
	     
# Legg merke til hvordan den ?kte n?yaktigheten faktisk realiseres og ikke bare ?dlegges
# av avrundingsfeil.