?# Program som tester en enkel metode for numerisk derivasjon

from math import *

def f(x):
	return sin(x)

def df(a,h):
	return (f(a+h)-f(a))/h 
	
	
for i in range(0,20):
	fderapp = df(0.5,10**(-i))
	fder = cos(0.5)
	print ('h=%e    df=%f   e=% 3.6e'
	%(10.**(-i), fderapp, fder-fderapp))