Index: mttroot/mtt/bin/trans/m/dm2fr.m ================================================================== --- mttroot/mtt/bin/trans/m/dm2fr.m +++ mttroot/mtt/bin/trans/m/dm2fr.m @@ -8,10 +8,13 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ +// Revision 1.6 1996/08/24 14:22:23 peter +// Put in a ; to avoid excessive log output. +// % %% Revision 1.5 1996/08/16 14:26:37 peter % %% Check and fix size of u0. % %% % %% Revision 1.4 1996/08/15 12:50:51 peter % %% Put in a conj to undo effect of transpose. @@ -50,11 +53,11 @@ [n,m]=size(u0); if m>n u0=u0'; end; - +j = sqrt(-1); fr = zeros(N,Ny); i = 0; for w = W' i = i+1; FR = C*( (E*j*w - A) \ B*u0 ) + D*u0;