fixed definition of function f in divdiff.hs