Yo no me acuerdo de donde saqué el código si de numerical recipes o de otro lado. Para el médodo de Householder hay varias rutinas diferentes y también para jacobi.
Yo al código lo probé y comparé los valores y andaba bien, la primera vez lo hice con una rutina que traduje de fortran a c y todos los cálculos los hacía en una dll, el programa estaba hecho en visual basic porque en esa época no sabía nada de delphi.
Después usé el código de nuemrical recipes que está en pascal pero no creo que haya errores.
Si tenés tempo un día de estos vemos que está pasando. Podemos ahcer pruebas para ver que margen de error hay.
Si tenés dudas el código de jean pierre moreau para el cálculo de vectores y valores propios, da el código en c, pascal, fortran y basic.
Es fácil implementarlo en un programa delphi en consola o se puede hacer en forma gráfica.
Saludos
|