Método de Jacobi en Lenguaje Python

Método de Jacobi en Lenguaje Python

Español | Online | Publicación:2013 | Free

Descripción:
Código que ejecuta el Método de Jacobi (para solución de sistemas de ecuaciones (éste código funciona para sistemas de [3x3]). Para garantizar la convergencia del método, hay que verificar que la matriz del sistema sea una matriz de diagonal estrictamente dominante.


El código es esl siguiente:
# -*- coding: UTF-8 -*- import os import decimal def Xval(a1, b1, c1, d1, y, z): # x = (-(b1*y)-(c1*z)+d1)/a1 # return x # # def Yval(a2, b2, c2, d2, x, z): # Formulas para recursividad y = (-(a2*x)-(c2*z)+d2)/b2 # return y # # def Zval(a3, b3, c3, d3, x, y): # z = (-(a3*x)-(b3*y)+d3)/c3 # return z # def main(): os.system("cls") vals = [[0.0,0.0,0.0],[1.0,1.0,1.0]] cons = [[0,0,0,0],[0,0,0,0],[0,0,0,0]] loop_val = 0 print "\n METODO DE JACOBI\n\nA1x + B1y + C1z = D1\n\ A2x + B2y + C2z = D2\nA3x + B3y + C3z = D3\n\n" cons[0][0] = input("\nIngrese A1 > ") # Constantes para ecuacion 1 cons[0][1] = input("Ingrese B1 > ") cons[0][2] = input("Ingrese C1 > ") cons[0][3] = input("Ingrese D1 > ") cons[1][0] = input("\nIngrese A2 > ") # Constantes para ecuacion 2 cons[1][1] = input("Ingrese B2 > ") cons[1][2] = input("Ingrese C2 > ") cons[1][3] = input("Ingrese D2 > ") cons[2][0] = input("\nIngrese A3 > ") # Constantes para ecuacion 3 cons[2][1] = input("Ingrese B3 > ") cons[2][2] = input("Ingrese C3 > ") cons[2][3] = input("Ingrese D3 > ") os.system("cls") print "Calculando..." while(vals[0] != vals [1]): print vals[0] vals[1][0] = Xval(cons[0][0], cons[0][1], cons[0][2], cons[0][3],\ vals[0][1], vals[0][2]) vals[1][1] = Yval(cons[1][0], cons[1][1], cons[1][2], cons[1][3],\ vals[0][0], vals[0][2]) vals[1][2] = Zval(cons[2][0], cons[2][1], cons[2][2], cons[2][3],\ vals[0][0], vals[0][1]) print vals[1] vals[0][0] = Xval(cons[0][0], cons[0][1], cons[0][2], cons[0][3],\ vals[1][1], vals[1][2]) vals[0][1] = Yval(cons[1][0], cons[1][1], cons[1][2], cons[1][3],\ vals[1][0], vals[1][2]) vals[0][2] = Zval(cons[2][0], cons[2][1], cons[2][2], cons[2][3],\ vals[1][0], vals[1][1]) raw_input("\n\nMetodo de Jacobi terminado...") main() exit()