new Jacobian implementation which runs in O(n*F) instead of O(n^2 *F) as the previous one!