simple numerical differentiation function with simple functions for now, will have to tune for multivariate system of odes