ment">% calculate j*omega</span>
0057 jw=j*2*pi*f;
0058 
0059 <span class="comment">% Allocate H</span>
0060 N=length(f);
0061 D=length(outdof);
0062 R=length(indof);
0063 Nm=length(p);
0064 H=zeros(N,D,R);
0065 <span class="comment">% Loop through all inputs</span>
0066 <span class="keyword">for</span> inno = 1:R
0067     <span class="comment">% Loop through response dofs</span>
0068     <span class="keyword">for</span> outno = 1:D
0069         <span class="comment">% Loop through modes</span>
0070        