Index: mttroot/mtt/bin/trans/m/dm2sr.m
==================================================================
--- mttroot/mtt/bin/trans/m/dm2sr.m
+++ mttroot/mtt/bin/trans/m/dm2sr.m
@@ -8,10 +8,13 @@
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % %% Version control history
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % %% $Id$
 % %% $Log$
+% %% Revision 1.1  1996/08/11 09:42:40  peter
+% %% Initial revision
+% %%
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 
 [N,M] = size(T);
 if M>N
@@ -27,8 +30,9 @@
 
 sr = zeros(N,NN);
 i = 0;
 for t = T'
   i=i+1;
-  sr(i,:) = reshape(C*(1-exp(A*t))*B + D*ones(size(t)), 1,NN);
+  SR = C*( A\(exp(A*t)-1) )*B + D*ones(size(t));
+  sr(i,:) = reshape(SR, 1,NN);
 end;