using NumericFuns