Code des valeurs propres d'une matrice carrée en Fortran 90 ou algo

Bonjour à tous, je dois coder un programme prenant en entrée une matrice (donc je demande à l'utilisateur de rentrer la taille, souvent de 2 ou 3, ainsi que les composantes une par une), et en sortie ses valeurs propres.

Pour l'entrée de la matrice c'est ok, mais je ne vois pas du tout le code pour avoir les valeurs propres ensuite (sachant qu'il faut avoir le polynôme caractéristique etc, ou y a t il une façon plus simple?)

... Quelqu'un pourrait t-il m'aider svp? Même le langage algorithmique m'aiderait ...