## Parte 1. Introduce un Sistema de Ecuaciones Lineales


```python
#importamos las bibliotecas que utilizaremos
%matplotlib inline
from matplotlib import pyplot as plt
from IPython.display import display, Math, Latex, Markdown
from sympy.interactive import printing
from sympy import Eq, solve_linear_system, Matrix,Symbol, init_printing, init_session
from numpy import linalg
import numpy as np
import sympy as sp
printing.init_printing(use_latex=True)
init_printing(use_latex='mathjax')

#Definimos las variables que se utilizarán guardar datos

def printmd(string):
    display(Markdown(string))

#Definimos las variables para los coeficientes

x, y=sp.symbols ('x,y')
a,b,c,d,e,f,g,h,r,j,k,l,m,n,a1,b1,c1=sp.symbols ('a,b,c,d,e,f,g,h,r,j,k,l,m,n,a1,b1,c1')

to= r'$$Introduce\ los\ coeficientes\ de\ tu\ primera\ ecuación\ lineal,\ de\ la\ forma\ ax+by=c$$'

printmd(to)

#Para que Python interprete la entrada como un número entero, se debe utilizar la función int() para decimal se usa float()

a= int(input("Cuál es el coeficiente de x: "))
a
b=int(input("Cuál es el coeficiente de y: "))
b
c=int(input("Cuál es el término constante: "))

#
printmd('$$Tu\ ecuación\ 1\ es:$$')

ec1= sp.Function('ec1')
ec2= sp.Function('ec2')

ec1= Eq(a*x + b*y, c)

display(ec1)

ta= r'$$Introduce\ los\ datos\ de\ la\ segunda\ ecuación\ de\ la\ forma\ dx+ey=f:$$'

printmd(ta)

d= int(input("Cuál es el coeficiente de x: "))
d
e=int(input("Cuál es el coficiente de y: "))
e
f=int(input("Cuál es el termino constante: "))

printmd('$$Tu\ ecuación\ 2\ es:$$')

ec2= Eq(d*x + e*y, f)

display(ec2)

printmd ('$$Por\ lo\ tanto\ tu\ sistema\ es:$$')

display(ec1, ec2)

printmd('$$Por\ lo\ tanto\ el\ arreglo\ rectangular\ de\ tu\ sistema\ es:$$')

row1=[a,b,c]
row2=[d,e,f]

system = Matrix ((row1,row2))

display (system)

#Gráficamos el sistema

printmd ('**La gráfica del sistema es:**')

def y1(x):
    return (c - a*x)/b 
# Función lineal.

def y2(x):
    return  (f - d*x)/f 
# Valores del eje X que toma el gráfico.
x = range(-10, 15)
# Graficar ambas funciones.
plt.plot(x, [y1(i) for i in x])
plt.plot(x, [y2(i) for i in x])
# Establecer el color de los ejes.
plt.axhline(0, color="black")
plt.axvline(0, color="black")
# Limitar los valores de los ejes.
plt.xlim(-15, 15)
plt.ylim(-15, 15)
# Guardar gráfico como imágen PNG.
plt.savefig("output.png")
# Mostrarlo.
plt.show()
```

## Parte 2. Modifica un Sistema de Ecuaciones Lineales

**-Tomando en cuenta el sistema que introdujiste anteriormente:**

**-Elije una ecuación y multiplicala por un número entero**

**-Observa cómo cambia el sistema respecto al original**


```python
x,y = sp.symbols('x,y')
m=int(input('Qué ecuación quieres multiplicar por un número (introduce 1 o 2 según la ecuación que desees multiplicar):'))
m

if m == 1:
    n=int(input('Por qué número deseas multiplicar: '))
    n
    
    a1=a*n
    b1=b*n
    c1=c*n
    
    printmd('**Tu ecuación 1 ahora es:**')
    
    display(ec1*n)

    eq3 = sp.Function('eq3')
    
    printmd('**Es decir:**')

    eq3 = Eq(a1*x+b1*y, c1)

    display(eq3)
    
    printmd('**Por lo tanto tu sistema es el siguiente:**')

    display (eq3,ec2)
    
    printmd('$$Por\ lo\ tanto\ el\ arreglo\ rectangular\ de\ tu\ sistema\ es:$$')

    row3=[a1,b1,c1]
    row4=[d,e,f]

    system1 = Matrix ((row3,row4))

    display (system1)

    #Gráficamos el sistema

    printmd ('**La gráfica del sistema es:**')

    def y11(x):
        return (c1 - a1*x)/b1 
    # Función lineal.

    def y2(x):
        return  (f - d*x)/f 
    # Valores del eje X que toma el gráfico.
    x = range(-10, 15)
    # Graficar ambas funciones.
    plt.plot(x, [y1(i) for i in x])
    plt.plot(x, [y2(i) for i in x])
    # Establecer el color de los ejes.
    plt.axhline(0, color="black")
    plt.axvline(0, color="black")
    # Limitar los valores de los ejes.
    plt.xlim(-15, 15)
    plt.ylim(-15, 15)
    # Guardar gráfico como imágen PNG.
    plt.savefig("output.png")
    # Mostrarlo.
    plt.show()

elif m==2: 
    n=int(input('Por qué número deseas multiplicar: '))
    n
    
    d1=d*n
    e1=e*n
    f1=f*n
    
    
    printmd('**Tu ecuación 2 ahora es:**')
    
    display(ec2*n)
        
    eq4 = sp.Function('eq4')
    
    printmd('**Es decir:**')

    eq4 = Eq(d1*x+e1*y, f1)

    display(eq4)
    
    printmd('**Por lo tanto tu sistema es el siguiente:**')

    display (ec1,eq4)
    
    printmd('$$Por\ lo\ tanto\ el\ arreglo\ rectangular\ de\ tu\ sistema\ es:$$')

    row4=[a,b,c]
    row5=[d1,e1,f1]

    system2 = Matrix ((row4,row5))

    display (system2)

    #Gráficamos el sistema

    printmd ('**La gráfica del sistema es:**')

    def y1(x):
        return (c - a*x)/b 
    # Función lineal.

    def y22(x):
        return  (f1 - d1*x)/f1 
    # Valores del eje X que toma el gráfico.
    x = range(-10, 15)
    # Graficar ambas funciones.
    plt.plot(x, [y1(i) for i in x])
    plt.plot(x, [y2(i) for i in x])
    # Establecer el color de los ejes.
    plt.axhline(0, color="black")
    plt.axvline(0, color="black")
    # Limitar los valores de los ejes.
    plt.xlim(-15, 15)
    plt.ylim(-15, 15)
    # Guardar gráfico como imágen PNG.
    plt.savefig("output.png")
    # Mostrarlo.
    plt.show()

else:
    printmd('**Introdujiste un número distinto de 1 o 2, vuelve a intentarlo**')
    
```

## Parte 3. Muestra la solución de un SEL


```python

```

## Parte 4. Introduce un Parámetro para resolver un SEL


```python

```


```python

```

255023a7086fdddfc3559b836b9bd74a04dc3368
ipynb
Jupyter Notebook
Programa SEL.ipynb
LuisEnrique3149/SEL_soluciones
c615dbb9b2cd984193e5ad666f6d4f81705c948b
MIT
Programa SEL.ipynb
LuisEnrique3149/SEL_soluciones
c615dbb9b2cd984193e5ad666f6d4f81705c948b
MIT
Programa SEL.ipynb
LuisEnrique3149/SEL_soluciones
c615dbb9b2cd984193e5ad666f6d4f81705c948b
MIT
Qwen/Qwen-72B
1. YES
2. YES