function [gk, Gk] = nonlinear_ok_simple_g(k, xk1) n = size(xk1, 1); if k == 1 Gk = zeros(n, n); gk = Gk * xk1 + 1; else Gk = eye(n); gk = Gk * xk1 + 1; end return end