Código de Python - Regla de Cramer

Regla de Cramer

Python

(3)

Publicado el 4 de Marzo del 2013 por Abraham (21 códigos)

24.198 visualizaciones desde el 4 de Marzo del 2013

Código fuente que resuelve sistemas de ecuaciones de 3x3 por medio de la Regla de Cramer.

Versión 1
(3)

Publicado el 4 de Marzo del 2013

gráfica de visualizaciones de la versión: Versión 1

24.199 visualizaciones desde el 4 de Marzo del 2013

# -*- coding: UTF-8 -*-

# Regla de Cramer

# Abraham Teran

import os

def sarrus(A):

    val = ((A[0][0]*A[1][1]*A[2][2])+\

           (A[0][1]*A[1][2]*A[2][0])+\

           (A[0][2]*A[1][0]*A[2][1]))-\

          ((A[2][0]*A[1][1]*A[0][2])+\

           (A[2][1]*A[1][2]*A[0][0])+\

           (A[2][2]*A[1][0]*A[0][1]))

    return val

def main():

    os.system("cls")

    sismat = [[0.0,0.0,0.0,0.0],\

              [0.0,0.0,0.0,0.0],\

              [0.0,0.0,0.0,0.0]]

    res = [0.0,0.0,0.0]

    print "\n\tREGLA DE CRAMER\n\n\nA1x + B1y + C1z = D1"

    print "A2x + B2y + C2z = D2\nA3x + B3y + C3z = D3\n\nIngresar valores:\n"

    sismat[0][0] = input("Ingrese A1 > ")

    sismat[0][1] = input("Ingrese B1 > ")

    sismat[0][2] = input("Ingrese C1 > ")

    sismat[0][3] = input("Ingrese D1 > ")

    sismat[1][0] = input("\nIngrese A2 > ")

    sismat[1][1] = input("Ingrese B2 > ")

    sismat[1][2] = input("Ingrese C2 > ")

    sismat[1][3] = input("Ingrese D2 > ")

    sismat[2][0] = input("\nIngrese A3 > ")

    sismat[2][1] = input("Ingrese B3 > ")

    sismat[2][2] = input("Ingrese C3 > ")

    sismat[2][3] = input("Ingrese D3 > ")

    mat_x = [[sismat[0][3],sismat[0][1],sismat[0][2]],\

             [sismat[1][3],sismat[1][1],sismat[1][2]],\

             [sismat[2][3],sismat[2][1],sismat[2][2]]]

    mat_y = [[sismat[0][0],sismat[0][3],sismat[0][2]],\

             [sismat[1][0],sismat[1][3],sismat[1][2]],\

             [sismat[2][0],sismat[2][3],sismat[2][2]]]

    mat_z = [[sismat[0][0],sismat[0][1],sismat[0][3]],\

             [sismat[1][0],sismat[1][1],sismat[1][3]],\

             [sismat[2][0],sismat[2][1],sismat[2][3]]]

    det_mat = sarrus(sismat)

    if (det_mat == 0):

        raw_input("Determinante de A nulo...")

    else:

        det_matx = sarrus(mat_x)

        det_maty = sarrus(mat_y)

        det_matz = sarrus(mat_z)

        res[0] = det_matx/det_mat

        res[1] = det_maty/det_mat

        res[2] = det_matz/det_mat

        print "\nP => ",res

        raw_input("\n\nRegla de Cramer terminado...")

main()

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

jobran

18 de Octubre del 2016

hola tiene un error y como estoy comenzando a usar python no veo el error me puedes ayudar

Responder

Abraham Terán

20 de Octubre del 2016

¿Qué versión de Python estás utilizando? Creo que lo programé usando Python 2.6

Responder

kevin

14 de Agosto del 2022

import os

def sarrus(A):
val = ((A[0][0]*A[1][1]*A[2][2])+\
(A[0][1]*A[1][2]*A[2][0])+\
(A[0][2]*A[1][0]*A[2][1]))-\
((A[2][0]*A[1][1]*A[0][2])+\
(A[2][1]*A[1][2]*A[0][0])+\
(A[2][2]*A[1][0]*A[0][1]))
return val

def main():
os.system("cls")
sismat = [[0.0,0.0,0.0,0.0],\
[0.0,0.0,0.0,0.0],\
[0.0,0.0,0.0,0.0]]
res = [0.0,0.0,0.0]
print ("\n\tREGLA DE CRAMER\n\n\nA1x + B1y + C1z = D1")
print ("A2x + B2y + C2z = D2\nA3x + B3y + C3z = D3\n\nIngresar valores:\n")
sismat[0][0] = int(input("Ingrese A1 > "))
sismat[0][1] = int(input("Ingrese B1 > "))
sismat[0][2] = int(input("Ingrese C1 > "))
sismat[0][3] = int(input("Ingrese D1 > "))
sismat[1][0] = int(input("\nIngrese A2 > "))
sismat[1][1] = int(input("Ingrese B2 > "))
sismat[1][2] = int(input("Ingrese C2 > "))
sismat[1][3] = int(input("Ingrese D2 > "))
sismat[2][0] = int(input("\nIngrese A3 > "))
sismat[2][1] = int(input("Ingrese B3 > "))
sismat[2][2] = int(input("Ingrese C3 > "))
sismat[2][3] = int(input("Ingrese D3 > "))
mat_x = [[sismat[0][3],sismat[0][1],sismat[0][2]],\
[sismat[1][3],sismat[1][1],sismat[1][2]],\
[sismat[2][3],sismat[2][1],sismat[2][2]]]
mat_y = [[sismat[0][0],sismat[0][3],sismat[0][2]],\
[sismat[1][0],sismat[1][3],sismat[1][2]],\
[sismat[2][0],sismat[2][3],sismat[2][2]]]
mat_z = [[sismat[0][0],sismat[0][1],sismat[0][3]],\
[sismat[1][0],sismat[1][1],sismat[1][3]],\
[sismat[2][0],sismat[2][1],sismat[2][3]]]
det_mat = sarrus(sismat)
if (det_mat == 0):
print("Determinante de A nulo...")
else:
det_matx = sarrus(mat_x)
det_maty = sarrus(mat_y)
det_matz = sarrus(mat_z)
res[0] = det_matx/det_mat
res[1] = det_maty/det_mat
res[2] = det_matz/det_mat
print ("\nP => ",res)
print("\n\nRegla de Cramer terminado...")

main()

Responder