Décomposition d'une matrice en valeurs singuileres

Bonjour,
je cherche à mettre en place un algorithme pour obtenir (numeriquement) les valeurs singulieres d'une matrice carree, methode connue aussi sous le nom de SVD.
Quelles sont les methodes aisement programmables pour ca? Si quelqu un a un code "exemple" je suis preneur (de toute facon je ne le reprendrais pas tel quel, mais je le reecrirais en OCaml).
Merci pour votre aide.

ha et tu ne pourras evidemment obtenir au mieux que des valeurs approchées de ces valeurs propres car comme chacun sait si la taille de ta matrice est supérieur a 5, ses valeurs propres qui sont solutions d'un polynome de degré >=5 ne se mettent pas en général sous la forme de radicaux, cf notre ami galois et la non resolubilité des groupes qui portent son nom pour n>=5

Merci pour toutes ces indications, je ne cherche en effet que des valeurs approchees, donc le probleme ne se pose pas.