disp(' ');
disp(' peak_response_random.m  ver 1.1  December 3, 2012 ');
disp(' by Tom Irvine ');
disp(' ');
disp(' Consider a single-degree-of-freedom system subjected ');
disp(' to a random vibration base input.  This script calculated ');
disp(' the maximum expected peak response in terms of the sigma value. ');
disp(' ');
disp(' This is based on the Rayleigh distribution of the peaks. ');
disp(' ');
%
disp(' Enter the natural frequency(Hz) ');
fn=input(' ');
disp(' ');
disp(' Enter the input duration(sec) ');
T=input(' ');
%
c=sqrt(2*log(fn*T));
ps=c + 0.5772/c;
%
disp(' ');
out1=sprintf(' The maximum expected peak response is: %4.3g sigma ',ps);
disp(out1);