function [AMP] = SMDRFsweep(k, m, w, e, ti, tf)     % ‘SDSF’ for ‘Spring-Mass Stationary Frame’

% AMP = SMDRFsweep(104000,1/16,[1;5;10;15;20;21;22;23;24;25;30;35;40;45;50],1,0,300)

[i,~] = size(w);

for ii=1:i

    [t, y] = SMDRF(k,m,w(ii),e,ti,tf);
    
    AMP(ii) = max(y(:,1).*(t(:,1)>0.95*tf));

end

plot(w,AMP,'-or')

end