
function MS_S = ms(gp, kc, ti)

    [gp_n gp_d] = tfdata(gp,'v');
    
    delay = totaldelay(gp);

    w = logspace(-4,4,4000);

    Gol_n_s = conv([ti*kc kc],gp_n);
    Gol_d_s = conv([ti 0],gp_d);
    
    Gol_n_w = polyval(Gol_n_s,w*1i);
    Gol_d_w = polyval(Gol_d_s,w*1i);
    delay_w = exp(-delay*w*1i);    
    
    Gol = Gol_n_w./Gol_d_w.*delay_w;
    
    S = 1./abs(1+Gol);
    MS_S = max(S);
end