       disp(' ')
       disp(' Miles.m     ver 1.1    February 15, 2005')
	   disp(' by Tom Irvine ')
   	   disp(' Email: tomirvine@aol.com ')
	   disp(' ')
	   disp(' This program calculates the response of a single-degree-of-freedom ') 
	   disp(' system to a base input power spectral density input via Miles equation. ')
	   disp(' ')
       %
       clear miles_acc;
       clear miles_rd;
       clear f;
       clear a;
       clear q;
       clear om;
       clear om2;
       %
       disp(' Enter the natural frequency (Hz) ')
	   f=input(' ');
       %
       disp(' ')   
       disp(' Enter the PSD (G^2/Hz) value at the natural frequency ')
	   a=input(' ');
       disp(' ')
       disp(' Enter the Q value (typically Q=10) ');
	   q=input(' ');
       %
       miles_acc=sqrt((pi/2.)*a*f*q);
       %
	   out1=sprintf('\n Acceleration Response = %10.2f GRMS (1 sigma) ',miles_acc);
	   out2=sprintf('                       = %10.2f G    (3 sigma) ',3.*miles_acc);
	   disp(out1)
       disp(out2)
       %
       om = 2.*pi*f;
       om2 = (om^2.);
       miles_rd=miles_acc/om2;
	   miles_rd=miles_rd*386.;
       %
	   out3=sprintf('\n Relative Displacement = %10.3g inch (1 sigma) ',miles_rd);
	   out4=sprintf('                       = %10.3g inch (3 sigma) ',3.*miles_rd);
	   disp(out3)
       disp(out4)
       %
       miles_rd=miles_rd*25.4;
       %
	   out5=sprintf('\n Relative Displacement = %10.3g mm (1 sigma) ',miles_rd);
	   out6=sprintf('                       = %10.3g mm (3 sigma) ',3.*miles_rd);
       disp(out5)
       disp(out6)