Fortran Fortran mapa de henon

[AngelAT]

que algoritmo se podría usar para calcular los exponentes de liapunov en el mapa de Henon en Fortran 90, que teoría necesito conocer antes de programarlo
(se que tiene que ver con la matriz jacobiana y sus autovalores)? .....para mapas unidimensionales es sencillo pero en este caso me enrede...