Buenas, estamos trancados en el obligatorio en la parte 2 de realizar el programa en octave que resuelva el sistema lineal. El problema que tenemos es que hay dos maneras de aproximar la resolución, una es resolverlo con la matriz esparsa estocástica por columnas y otra es con la matriz ya perturbada. El beneficio de usar la estocástica sería que si se encuentra un método eficiente se puede explotar el hecho de que tiene gran cantidad de ceros, y de hecho encontramos varios métodos para matrices esparsas pero estos requieren que la matriz sea tridiagonal y simétrica, requisitos que no cumple la nuestra. Tampoco logramos encontrar métodos para llegar a una matriz de este tipo a partir de la matriz nuestra.
Por otro lado con la matriz perturbada no podemos aplicar Gauss Seidel o Jacobi porque el radio espectral es mayor que 1 así que no podemos asegurar la convergencia del método y de hecho hicimos un programa para comprobar esto.
Necestiaríamos saber si hay algún método para resolver matrices esparsas no tridiagonales ni siméticas o uno general que pueda asegurar convergencia en nuestras condiciones.
Desde ya muchas gracias