lcmdiag
=======

least common multiple diagonal factorization



Calling Sequence
~~~~~~~~~~~~~~~~


::

    [N,D]=lcmdiag(H)
    [N,D]=lcmdiag(H,flag)




Arguments
~~~~~~~~~

:H rational matrix
: :N polynomial matrix
: :D diagonal polynomial matrix
: :flag character string: `'row'` or `'col'` (default)
:



Description
~~~~~~~~~~~

`[N,D]=lcmdiag(H,'row')` computes a factorization `D*H=N`, i.e.
`H=D^(-1)*N` where D is a diagonal matrix with D(k,k)=lcm of kth row
of H('den').

`[N,D]=lcmdiag(H)` or `[N,D]=lcmdiag(H,'col)` returns `H=N*D^(-1)`
with diagonal D and D(k,k)=lcm of kth col of H('den')



Examples
~~~~~~~~


::

    s=`poly`_(0,'s');
    H=[1/s,(s+2)/s/(s+1)^2;1/(s^2*(s+2)),2/(s+2)];
    [N,D]=lcmdiag(H);
    N/D-H




See Also
~~~~~~~~


+ `lcm`_ least common multiple
+ `gcd`_ gcd calculation
+ `bezout`_ Bezout equation for polynomials or integers


.. _lcm: lcm.html
.. _gcd: gcd.html
.. _bezout: bezout.html