open Impact
let f(x) = x**x
 
a <- 3.0
print (f(a))
 
h <- 1.0
print (h, f(a+h), (f(a+h) - f(a))/h)
h <- 0.1
print (h, f(a+h), (f(a+h) - f(a))/h)
h <- 0.01
print (h, f(a+h), (f(a+h) - f(a))/h)
h <- 0.001
print (h, f(a+h), (f(a+h) - f(a))/h)
h <- 0.0001
print (h, f(a+h), (f(a+h) - f(a))/h)