ADDED   mttroot/mtt/bin/trans/cse2csm_r
Index: mttroot/mtt/bin/trans/cse2csm_r
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/cse2csm_r
@@ -0,0 +1,130 @@
+#! /bin/sh
+
+     ###################################### 
+     ##### Model Transformation Tools #####
+     ######################################
+
+# Bourne shell script: cse2csm_r
+# Constrained-state equation to linear constrained-state matrices conversion
+# P.J.Gawthrop  6th September 1991, May 1994
+# Copyright (c) P.J.Gawthrop, 1991, 1994.
+
+###############################################################
+## Version control history
+###############################################################
+## $Id$
+## $Log$
+###############################################################
+
+# Inform user
+echo Creating $1_csm.r
+
+# Remove the old log file
+rm -f cse2csm_r.log
+
+# Use reduce to accomplish the transformation
+reduce >cse2csm_r.log << EOF
+
+in "$1_def.r";
+in "$1_cse.r";
+in "$1_cr.r";
+in "$1_sympar.r";
+
+OFF Echo;
+OFF Nat;
+
+% Find MTTA : the A matrix
+matrix MTTA(MTTNx,MTTNx);
+
+  FOR j := 1:MTTNx DO
+    BEGIN
+    xj := MTTX(j,1);
+    FOR i := 1:MTTNx DO
+      MTTA(i,j) := df(MTTEdx(i,1), xj, 1);
+    END;
+
+% Find MTTB : the B matrix
+matrix MTTB(MTTNx,MTTNu);
+
+  FOR j := 1:MTTNu DO
+    BEGIN
+    uj := MTTU(j,1);
+    FOR i := 1:MTTNx DO
+      MTTB(i,j) := df(MTTEdx(i,1), uj, 1);
+    END;
+
+% Find MTTC : the C matrix
+matrix MTTC(MTTNy,MTTNx);
+
+  FOR i := 1:MTTNy DO
+    FOR j := 1:MTTNx DO
+    BEGIN
+      xj := MTTX(j,1);
+      MTTC(i,j) := df(MTTY(i,1), xj, 1);
+    END;
+
+% Find MTTD : the D matrix
+matrix MTTD(MTTNy,MTTNu);
+
+  FOR i := 1:MTTNy DO
+    FOR j := 1:MTTNu DO
+    BEGIN
+      xj := MTTU(j,1);
+      MTTD(i,j) := df(MTTY(i,1), xj, 1);
+    END;
+
+%Substitute the ss values
+in "$1_ss.r";
+
+%Create the output file
+OUT "$1_csm.r";
+
+%Write out the matrices.
+IF MTTNx>0 THEN
+BEGIN
+  write "matrix MTTE(", MTTNx, ",", MTTNx, ");";
+  FOR i := 1:MTTNx DO
+   FOR j := 1:MTTNx DO IF MTTE(i,j) NEQ 0 THEN
+      write "MTTE(", i, ",", j, ") := ", MTTE(i,j);
+
+  write "matrix MTTA(", MTTNx, ",", MTTNx, ");";
+  FOR i := 1:MTTNx DO
+   FOR j := 1:MTTNx DO IF MTTA(i,j) NEQ 0 THEN
+      write "MTTA(", i, ",", j, ") := ", MTTA(i,j);
+END;
+
+IF MTTNx>0 THEN
+IF MTTNu>0 THEN
+BEGIN
+  write "matrix MTTB(", MTTNx, ",", MTTNu, ");";
+  FOR i := 1:MTTNx DO
+   FOR j := 1:MTTNu DO IF MTTB(i,j) NEQ 0 THEN
+      write "MTTB(", i, ",", j, ") := ", MTTB(i,j);
+END;
+
+%Write it out
+IF MTTNy>0 THEN
+IF MTTNx>0 THEN
+BEGIN
+  write "matrix MTTC(", MTTNy, ",", MTTNx, ");";
+  FOR i := 1:MTTNy DO
+    FOR j := 1:MTTNx DO IF MTTC(i,j) NEQ 0 THEN
+      write "MTTC(", i, ",", j, ") := ", MTTC(i,j);
+END;
+
+IF MTTNy>0 THEN IF MTTNu>0 THEN 
+BEGIN
+  write "matrix MTTD(", MTTNy, ",", MTTNu, ");";
+  FOR i := 1:MTTNy DO
+    FOR j := 1:MTTNu DO IF MTTD(i,j) NEQ 0 THEN
+      write "MTTD(", i, ",", j, ") := ", MTTD(i,j);
+END;
+
+write "END;";
+
+SHUT "$1_csm.r";
+
+EOF
+
+# Now invoke the standard error handling.
+mtt_error_r cse2csm_r.log

ADDED   mttroot/mtt/bin/trans/dm2tf_r
Index: mttroot/mtt/bin/trans/dm2tf_r
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/dm2tf_r
@@ -0,0 +1,82 @@
+#! /bin/sh
+
+     ###################################### 
+     ##### Model Transformation Tools #####
+     ######################################
+
+# Bourne shell script: dm2tf_r
+# Reduce  descriptor matrices to transfer function
+# P.J.Gawthrop 8th May 1991, Dec 1993, April 1994.
+# Copyright (c) P.J.Gawthrop, 1991, 1993, 1994.
+
+###############################################################
+## Version control history
+###############################################################
+## $Id$
+## $Log$
+###############################################################
+
+#Inform user
+echo Creating $1_tf.r
+
+# Remove the old log file
+rm -f dm2tf_r.log
+
+# Use reduce to accomplish the transformation
+reduce >dm2tf_r.log << EOF
+
+%ON FLOAT;
+
+IN "$1_def.r";
+IN "$1_dm.r";
+
+OFF Echo;
+OFF Nat;
+
+%create sE-A
+%MTT_SEA := s*MTTE-MTTA;
+
+%Find the denominator of the TF - det(sE-A);
+%comden := det(MTT_SEA);
+
+%Find the Adjoint transpose.
+%matrix AdjT(MTTNx,MTTNx);
+%FOR i := 1:MTTNx DO
+%  BEGIN
+%  FOR j := 1:MTTNx DO
+%    AdjT(i,j) := cofactor(MTT_SEA,i,j);
+%  END;
+
+%Adj := TP(AdjT);
+
+%Find the numerator matrix
+%Num := MTTC*Adj*MTTB + MTTD*comden;
+
+
+%Create the transfer function matrix
+MTTTF := MTTD;
+IF MTTNy>0 THEN MTTTF := MTTTF + (MTTC * ((s*MTTE-MTTA)^-1) * MTTB);
+
+%MTTTF := Num/comden;
+
+
+OUT "$1_tf.r";
+
+%Declare the transfer function matrix
+write "matrix MTTTF(", MTTNy, ",", MTTNu, ")$"$
+
+%And write it.
+%MTTTF := MTTTF;
+FOR i := 1:MTTNy DO
+  BEGIN
+  FOR j := 1:MTTNu DO
+    IF MTTTF(i,j) NEQ 0 THEN write "MTTTF(", i, ",", j, ") := ", MTTTF(i,j)$
+  END;
+
+write ";END;"$
+SHUT "$1_tf.r";
+
+EOF
+
+# Now invoke the standard error handling.
+mtt_error_r dm2tf_r.log