C2DPulse:

--------------------------------------------------------------------------
   Create a discrete time system assuming a pulse at the input.

   Create a discrete time system from a continuous system
   assuming a pulse at the input.

   Given
   .
   x = ax + bu

   Find f and g where

   x(k+1) = fx(k) + gu(k)

   Since version 1.
--------------------------------------------------------------------------
   Form:
   [f, g] = C2DPulse( a, b, T, tau )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   a            (n,n)  Continuous plant matrix
   b            (n,m)  Input matrix
   T            (1,1)  Time step

   -------
   Outputs
   -------
   f            (n,n)  Discrete plant matrix
   g            (n,m)  Discrete input matrix

--------------------------------------------------------------------------
	References: Van Loan, C.F., Computing Integrals Involving the Matrix
               Exponential, IEEE Transactions on Automatic Control
               Vol. AC-23, No. 3, June 1978, pp. 395-404.
--------------------------------------------------------------------------