polar form
[Ro,Theta]=polar(A)
:A real or complex square matrix : :Ro, real matrix : :Theta, real or complex matrix :
[Ro,Theta]=polar(A) returns the polar form of A i.e. A=Ro*expm(%i*Theta) Ro symmetric >=0 and Theta hermitian >=0.
A=`rand`_(5,5);
[Ro,Theta]=polar(A);
`norm`_(A-Ro*`expm`_(%i*Theta),1)