function[D,H] = linearize_new(x,y,f_i,M_t,n)

D_o = zeros(n,1);       % inizializzo il gradiente D (4x1)
H_o = zeros(n,n);       % inizializzo l'hessiano H (4x4)

for i = 1:M_t
    
    % gradiente
    D_o = D_o + (-y(i)*exp(-y(i)*f_i(i))/(1+exp(-y(i)*f_i(i)))*[x(i,:) 1]');
         
    % hessiano
    H_o = H_o + exp(-y(i)*f_i(i))/((1+exp(-y(i)*f_i(i)))^2)*[x(i,:) 1]'*[x(i,:) 1];
         
end
D = D_o;
H = H_o;
end
     