function imp=respd(impact,G1,nt)
%function imp=respd(impact,G1,nt)
[nv,ns]=size(impact);
imp=zeros(nv,ns,nt);
imp(:,:,1)=impact;
for i=2:nt
   imp(:,:,i)=G1*imp(:,:,i-1);
end