Programa SEL.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parte 1. Introduce un Sistema de Ecuaciones Lineales"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#importamos las bibliotecas que utilizaremos\n",
    "%matplotlib inline\n",
    "from matplotlib import pyplot as plt\n",
    "from IPython.display import display, Math, Latex, Markdown\n",
    "from sympy.interactive import printing\n",
    "from sympy import Eq, solve_linear_system, Matrix,Symbol, init_printing, init_session\n",
    "from numpy import linalg\n",
    "import numpy as np\n",
    "import sympy as sp\n",
    "printing.init_printing(use_latex=True)\n",
    "init_printing(use_latex='mathjax')\n",
    "\n",
    "#Definimos las variables que se utilizarán guardar datos\n",
    "\n",
    "def printmd(string):\n",
    "    display(Markdown(string))\n",
    "\n",
    "#Definimos las variables para los coeficientes\n",
    "\n",
    "x, y=sp.symbols ('x,y')\n",
    "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')\n",
    "\n",
    "to= r'$$Introduce\\ los\\ coeficientes\\ de\\ tu\\ primera\\ ecuación\\ lineal,\\ de\\ la\\ forma\\ ax+by=c$$'\n",
    "\n",
    "printmd(to)\n",
    "\n",
    "#Para que Python interprete la entrada como un número entero, se debe utilizar la función int() para decimal se usa float()\n",
    "\n",
    "a= int(input(\"Cuál es el coeficiente de x: \"))\n",
    "a\n",
    "b=int(input(\"Cuál es el coeficiente de y: \"))\n",
    "b\n",
    "c=int(input(\"Cuál es el término constante: \"))\n",
    "\n",
    "#\n",
    "printmd('$$Tu\\ ecuación\\ 1\\ es:$$')\n",
    "\n",
    "ec1= sp.Function('ec1')\n",
    "ec2= sp.Function('ec2')\n",
    "\n",
    "ec1= Eq(a*x + b*y, c)\n",
    "\n",
    "display(ec1)\n",
    "\n",
    "ta= r'$$Introduce\\ los\\ datos\\ de\\ la\\ segunda\\ ecuación\\ de\\ la\\ forma\\ dx+ey=f:$$'\n",
    "\n",
    "printmd(ta)\n",
    "\n",
    "d= int(input(\"Cuál es el coeficiente de x: \"))\n",
    "d\n",
    "e=int(input(\"Cuál es el coficiente de y: \"))\n",
    "e\n",
    "f=int(input(\"Cuál es el termino constante: \"))\n",
    "\n",
    "printmd('$$Tu\\ ecuación\\ 2\\ es:$$')\n",
    "\n",
    "ec2= Eq(d*x + e*y, f)\n",
    "\n",
    "display(ec2)\n",
    "\n",
    "printmd ('$$Por\\ lo\\ tanto\\ tu\\ sistema\\ es:$$')\n",
    "\n",
    "display(ec1, ec2)\n",
    "\n",
    "printmd('$$Por\\ lo\\ tanto\\ el\\ arreglo\\ rectangular\\ de\\ tu\\ sistema\\ es:$$')\n",
    "\n",
    "row1=[a,b,c]\n",
    "row2=[d,e,f]\n",
    "\n",
    "system = Matrix ((row1,row2))\n",
    "\n",
    "display (system)\n",
    "\n",
    "#Gráficamos el sistema\n",
    "\n",
    "printmd ('**La gráfica del sistema es:**')\n",
    "\n",
    "def y1(x):\n",
    "    return (c - a*x)/b \n",
    "# Función lineal.\n",
    "\n",
    "def y2(x):\n",
    "    return  (f - d*x)/f \n",
    "# Valores del eje X que toma el gráfico.\n",
    "x = range(-10, 15)\n",
    "# Graficar ambas funciones.\n",
    "plt.plot(x, [y1(i) for i in x])\n",
    "plt.plot(x, [y2(i) for i in x])\n",
    "# Establecer el color de los ejes.\n",
    "plt.axhline(0, color=\"black\")\n",
    "plt.axvline(0, color=\"black\")\n",
    "# Limitar los valores de los ejes.\n",
    "plt.xlim(-15, 15)\n",
    "plt.ylim(-15, 15)\n",
    "# Guardar gráfico como imágen PNG.\n",
    "plt.savefig(\"output.png\")\n",
    "# Mostrarlo.\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parte 2. Modifica un Sistema de Ecuaciones Lineales"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**-Tomando en cuenta el sistema que introdujiste anteriormente:**\n",
    "\n",
    "**-Elije una ecuación y multiplicala por un número entero**\n",
    "\n",
    "**-Observa cómo cambia el sistema respecto al original**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "x,y = sp.symbols('x,y')\n",
    "m=int(input('Qué ecuación quieres multiplicar por un número (introduce 1 o 2 según la ecuación que desees multiplicar):'))\n",
    "m\n",
    "\n",
    "if m == 1:\n",
    "    n=int(input('Por qué número deseas multiplicar: '))\n",
    "    n\n",
    "    \n",
    "    a1=a*n\n",
    "    b1=b*n\n",
    "    c1=c*n\n",
    "    \n",
    "    printmd('**Tu ecuación 1 ahora es:**')\n",
    "    \n",
    "    display(ec1*n)\n",
    "\n",
    "    eq3 = sp.Function('eq3')\n",
    "    \n",
    "    printmd('**Es decir:**')\n",
    "\n",
    "    eq3 = Eq(a1*x+b1*y, c1)\n",
    "\n",
    "    display(eq3)\n",
    "    \n",
    "    printmd('**Por lo tanto tu sistema es el siguiente:**')\n",
    "\n",
    "    display (eq3,ec2)\n",
    "    \n",
    "    printmd('$$Por\\ lo\\ tanto\\ el\\ arreglo\\ rectangular\\ de\\ tu\\ sistema\\ es:$$')\n",
    "\n",
    "    row3=[a1,b1,c1]\n",
    "    row4=[d,e,f]\n",
    "\n",
    "    system1 = Matrix ((row3,row4))\n",
    "\n",
    "    display (system1)\n",
    "\n",
    "    #Gráficamos el sistema\n",
    "\n",
    "    printmd ('**La gráfica del sistema es:**')\n",
    "\n",
    "    def y11(x):\n",
    "        return (c1 - a1*x)/b1 \n",
    "    # Función lineal.\n",
    "\n",
    "    def y2(x):\n",
    "        return  (f - d*x)/f \n",
    "    # Valores del eje X que toma el gráfico.\n",
    "    x = range(-10, 15)\n",
    "    # Graficar ambas funciones.\n",
    "    plt.plot(x, [y1(i) for i in x])\n",
    "    plt.plot(x, [y2(i) for i in x])\n",
    "    # Establecer el color de los ejes.\n",
    "    plt.axhline(0, color=\"black\")\n",
    "    plt.axvline(0, color=\"black\")\n",
    "    # Limitar los valores de los ejes.\n",
    "    plt.xlim(-15, 15)\n",
    "    plt.ylim(-15, 15)\n",
    "    # Guardar gráfico como imágen PNG.\n",
    "    plt.savefig(\"output.png\")\n",
    "    # Mostrarlo.\n",
    "    plt.show()\n",
    "\n",
    "elif m==2: \n",
    "    n=int(input('Por qué número deseas multiplicar: '))\n",
    "    n\n",
    "    \n",
    "    d1=d*n\n",
    "    e1=e*n\n",
    "    f1=f*n\n",
    "    \n",
    "    \n",
    "    printmd('**Tu ecuación 2 ahora es:**')\n",
    "    \n",
    "    display(ec2*n)\n",
    "        \n",
    "    eq4 = sp.Function('eq4')\n",
    "    \n",
    "    printmd('**Es decir:**')\n",
    "\n",
    "    eq4 = Eq(d1*x+e1*y, f1)\n",
    "\n",
    "    display(eq4)\n",
    "    \n",
    "    printmd('**Por lo tanto tu sistema es el siguiente:**')\n",
    "\n",
    "    display (ec1,eq4)\n",
    "    \n",
    "    printmd('$$Por\\ lo\\ tanto\\ el\\ arreglo\\ rectangular\\ de\\ tu\\ sistema\\ es:$$')\n",
    "\n",
    "    row4=[a,b,c]\n",
    "    row5=[d1,e1,f1]\n",
    "\n",
    "    system2 = Matrix ((row4,row5))\n",
    "\n",
    "    display (system2)\n",
    "\n",
    "    #Gráficamos el sistema\n",
    "\n",
    "    printmd ('**La gráfica del sistema es:**')\n",
    "\n",
    "    def y1(x):\n",
    "        return (c - a*x)/b \n",
    "    # Función lineal.\n",
    "\n",
    "    def y22(x):\n",
    "        return  (f1 - d1*x)/f1 \n",
    "    # Valores del eje X que toma el gráfico.\n",
    "    x = range(-10, 15)\n",
    "    # Graficar ambas funciones.\n",
    "    plt.plot(x, [y1(i) for i in x])\n",
    "    plt.plot(x, [y2(i) for i in x])\n",
    "    # Establecer el color de los ejes.\n",
    "    plt.axhline(0, color=\"black\")\n",
    "    plt.axvline(0, color=\"black\")\n",
    "    # Limitar los valores de los ejes.\n",
    "    plt.xlim(-15, 15)\n",
    "    plt.ylim(-15, 15)\n",
    "    # Guardar gráfico como imágen PNG.\n",
    "    plt.savefig(\"output.png\")\n",
    "    # Mostrarlo.\n",
    "    plt.show()\n",
    "\n",
    "else:\n",
    "    printmd('**Introdujiste un número distinto de 1 o 2, vuelve a intentarlo**')\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parte 3. Muestra la solución de un SEL"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parte 4. Introduce un Parámetro para resolver un SEL"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "hide_input": true,
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  },
  "latex_envs": {
   "LaTeX_envs_menu_present": true,
   "autoclose": true,
   "autocomplete": true,
   "bibliofile": "biblio.bib",
   "cite_by": "apalike",
   "current_citInitial": 1,
   "eqLabelWithNumbers": true,
   "eqNumInitial": 1,
   "hotkeys": {
    "equation": "Ctrl-E",
    "itemize": "Ctrl-I"
   },
   "labels_anchors": false,
   "latex_user_defs": false,
   "report_style_numbering": false,
   "user_envs_cfg": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}