File mttroot/mtt/bin/trans/m/dm2ir.m artifact 35a8da1aa8 part of check-in 67560661e3

function [t,ir] = dm2ir(A,B,C,D,E,tmax,tstep);

% Descriptor matrix to impulse response.
% NB At the moment - this assumes that E is unity .....

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


t=[0:tstep:tmax]';
ir = C*exp(A*t)*B;