packages/hmatrix/examples/deriv.hs


1
2
3
4
5
6
7
8

-- Numerical differentiation

import Numeric.GSL

d :: (Double -> Double) -> (Double -> Double)
d f x = fst $ derivCentral 0.01 f x

main = print $ d (\x-> x * d (\y-> x+y) 1) 1