In a previous post I demonstrated how to use R’s simple built-in symbolic engine to generate Jacobian and (pseudo)-Hessian matrices that make non-linear optimization perform much more efficiently. Another related application is Gaussian error propagation. Say you have data from a set of measurements in variables x and y where you know the corresponding measurement errors (dx...