function pars(ethwv,tmitvl,dmpr)
ethwv1;
tmitvl1;
dmpr1;
% function pars(ethwv,tmitvl,dmpr)
% Purpose: To draw pseudo-acceleration response spectrum figure in 6 s
% variable names:
% ethwv acceleration of earthquake wave record in one column form
% tmitvl time interval of the earthquake wave record
% dmpr damping ratio
t = 0:tmitvl:(length(ethwv)-1)*tmitvl;
Prd = 0:0.01:6; Prd = Prd;
omgar = 2*pi./Prd;
for i = 1:length(omgar)
SYS = tf(1,[1,2*dmpr*omgar(i),omgar(i)*omgar(i)]);
disp(i) = max(abs(lsim(SYS,ethwv,t)));
spd(i) = 2*pi.*disp(i)/Prd(i);
acc(i) = 2*pi.*spd(i)/Prd(i);
end
plot(Prd,acc);