Index: mttroot/mtt/bin/trans/m/dm2fr.m
==================================================================
--- mttroot/mtt/bin/trans/m/dm2fr.m
+++ mttroot/mtt/bin/trans/m/dm2fr.m
@@ -8,13 +8,16 @@
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % %% 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.7  1996/11/06  16:40:38  peterg
+% %% Explicit definition of j
+% %%
+% %% 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.
@@ -58,11 +61,12 @@
 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;
+%%  FR = C*( (E*j*w - A) \ B*u0 ) + D*u0;
+  FR = C*( inv(E*j*w - A)*B*u0 ) + D*u0;
   fr(i,:) = conj(FR');
 end;