Código de Python - Método de Gauss-Seidel

Método de Gauss-Seidelgráfica de visualizaciones


Python

estrellaestrellaestrellaestrellaestrella(1)
Publicado el 06 de Junio del 2013 por Abraham
10.750 visualizaciones desde el 06 de Junio del 2013. Una media de 63 por semana
Código de ejecuta el Método de Gauss-Seidel para resolución de sistemas de ecuaciones lineales de tres ecuaciones por tres variables. Para garantizar que el método funcione, el sistema de ecuaciones lineales a resolver debe tener una matriz de diagonal estrictamente dominante.
Desarrollado con Python 2.7

Versión 1
estrellaestrellaestrellaestrellaestrella(1)

Publicado el 06 de Junio del 2013gráfica de visualizaciones de la versión: Versión 1
10.751 visualizaciones desde el 06 de Junio del 2013. Una media de 63 por semana
estrellaestrellaestrellaestrellaestrella
estrellaestrellaestrellaestrella
estrellaestrellaestrella
estrellaestrella
estrella

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
# -*- coding: UTF-8 -*-
 
# Metodo de Gauss-Seidel
# Abraham Teran
 
# NOTA: Para garantizar la convergencia del Metodo de Gauss-Seidel
#       hay que verificar que la matriz del sistema de ecuaciones
#       sea una matriz de diagonal estrictamente dominante.
#       Este codigo funciona para sistemas de ecuaciones
#       de tres ecuaciones por tres variables.
 
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():
    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]]\
 
    print "\n\tMETODO DE GAUSS-SEIDEL\n\nA1x + B1y + C1z = D1\n\
A2x + B2y + C2z = D2\nA3x + B3y + C3z = D3\n\n"\
 
    for i in range(3):
        print " ** FILA NUMERO",i+1,"**"
        for j in range(4):
            print "Dato [",i+1,",",j+1,"] de Matriz"
            cons[i][j] = float(input("INGRESE DATO > "))
            j = j+1
        i = i+1
        print "\n"\
 
    while vals[0] != vals[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[0][0], vals[1][2])
        vals[0][2] = Zval(cons[2][0], cons[2][1], cons[2][2], cons[2][3],\
vals[0][0], vals[0][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[1][0], vals[0][2])
        vals[1][2] = Zval(cons[2][0], cons[2][1], cons[2][2], cons[2][3],\
vals[1][1], vals[1][2])\
 
    raw_input("\n\nMetodo de Gauss-Seidel terminado...")
 
main()



Comentarios sobre la versión: Versión 1 (1)

maycol
07 de Octubre del 2014
estrellaestrellaestrellaestrellaestrella
No ha dejado ningún comentario
Responder

Comentar la versión: Versión 1

Nombre
Correo (no se visualiza en la web)
Valoración
Comentarios

http://lwp-l.com/s2391