Resolución de sistemas de ecuaciones con Sympy

Hasta ahora se tienen los siguientes artículos sobre el uso de la librería Sympy:

• Cálculo de límites con la librería Sympy
• Cálculo de derivadas con la librería Sympy
• Cálculo de integrales con la librería Sympy

En este caso se explicará el uso de la librería Sympy en la resolución de sistemas de ecuaciones, se utilizará el ejemplo de un artículo anterior (Resolución de sistemas de ecuaciones lineales por descomposición QR usando Numpy).

Se inicia ipython notebook (si desean usar el archivo lo pueden bajar del siguiente enlace): 

A continuación se muestra el código del script:

#!/usr/bin/env python







#Se importa sympy

from sympy import *

#Se define las variables simbolicas x,y,z

x = Symbol('x')

y = Symbol('y')

z= Symbol('z')










#Resolver el sistema de ecuaciones

#3x+9y-10z  =   24

#x-6y+4z      =   -4

#10x-2y+8z  =  20

resultado =solve([3*x+9*y-10*z-24,x-6*y+4*z+4,10*x-2*y+8*z-20],[x,y,z])




print resultado

Al ejecutar el script se obtiene lo siguiente:

{x: 308/103, z: -117/103, y: 42/103}

Se nota que el resultado obtenido es el mismo del artículo ya mencionado.

La imagen de la utilización de notebook se muestra a continuación: