Look into using Taylor polynomials for transcendental functions
The multivariate case seems interesting.
If we implement them it would also be interesting if we could make them parametric in the degree of approximation (i.e. n for an n-th order taylor polynomial)