polar - forme polaire
renvoie la forme polaire deA c'est à dire :
A=Ro*expm(%i*Theta) Ro symétrique >=0 et Theta hermitienne >=0.
A=rand(5,5); [Ro,Theta]=polar(A); norm(A-Ro*expm(%i*Theta),1)