Index: mttroot/mtt/lib/comp/simple/R_eqn.m
==================================================================
--- mttroot/mtt/lib/comp/simple/R_eqn.m
+++ mttroot/mtt/lib/comp/simple/R_eqn.m
@@ -13,10 +13,14 @@
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % %% Version control history
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % %% $Id$
 % %% $Log$
+% %% Revision 1.3  1996/09/12 16:40:45  peter
+% %% Now does a default (when no constitutive relationship) with a unit
+% %% constitutive relationship on each port.
+% %%
 % %% Revision 1.2  1996/09/12 13:12:04  peter
 % %% Removed an fprintf that was a mistake.
 % %%
 % %% Revision 1.1  1996/09/12 11:00:56  peter
 % %% Initial revision
@@ -38,12 +42,12 @@
 for outport = 1:ports
   outcause = -bonds(outport,1);
   outnumber = bond_number(outport);
   if length(cr)>0
     incause = bonds(:,1);
-    inports = 1:ports;
-    inbonds = bond_number;
+    inport = 1:ports;
+    inbond = bond_number;
   else % default if no cr: unit cr on each port
     incause = bonds(outport,1);
     inport = outport;
     inbond = outnumber;
   end;