added finite differences approximation for Gauss-Newton matrix