function[fig_num]=two_dof_frf_apply_rba(m,k,fig_num,iu);
%
disp(' ');
disp(' Apply rigid-body acceleration? ');
disp(' 1=yes  2=no ');
irba = input(' ');
%
if(irba==1)
        clear x;
        clear aaa;
%        
        disp(' ');
        disp(' Enter acceleration (G) ');
        accel=input(' ');
%
        aaa=zeros(2,1);
%
        aaa(1)=accel;
        aaa(2)=accel;
%
        aaa=aaa*386;
        mm=m*aaa;
%
        x=k\mm;
%
        disp(' ');
        disp('   Displacement  ');
        disp('    x1        x2       ');
        disp('   (in)      (in)      ');
        out1=sprintf(' %6.2e  %6.2e ',x(1),x(2));
        disp(out1);
%
end   