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.2 1996/08/15 10:24:28 peter +% %% Includes u0 argument. +% %% % %% Revision 1.1 1996/08/10 14:11:28 peter % %% Initial revision % %% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -24,19 +27,23 @@ u0(1) = 1; end; fr = zeros(N,Ny); i = 0; -for w = W +for w = W' i = i+1; FR = C*( (E*j*w - A) \ B ) + D; fr(i,:) = FR'; end; + + + +