{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"jupyter": {
"outputs_hidden": true
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\begin{array}{l}\n",
"\\begin{array}{lcrcrcrcrcl}\n",
" \\min \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 2 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 3 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 5 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 15 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 14 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 16 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 12 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 700 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 19 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 20 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 15 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 8 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 900 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 x_{1} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 5 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 6 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 10 \\\\\n",
"\\end{array} \\\\\n",
"x_{1}, x_{2}, x_{3}, x_{4} \\geq 0\n",
"\\end{array}$$"
],
"text/plain": [
"LP problem (use 'view(...)' or '%display typeset' for details)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\frac{42950}{219}$$"
],
"text/plain": [
"42950/219"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\frac{1750}{73},\\,0,\\,\\frac{6490}{219},\\,0\\right)$$"
],
"text/plain": [
"(1750/73, 0, 6490/219, 0)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\begin{array}{l}\n",
"\\begin{array}{lcrcrcrcrcl}\n",
" \\min \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 2 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 3 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 5 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 15 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 14 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 16 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 12 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 700 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 19 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 20 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 15 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 8 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 900 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 x_{1} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 5 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 6 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 10 \\\\\n",
"\\end{array} \\\\\n",
"x_{1}, x_{2}, x_{3}, x_{4} \\geq 0\n",
"\\end{array}$$"
],
"text/plain": [
"LP problem (use 'view(...)' or '%display typeset' for details)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\frac{42950}{219}$$"
],
"text/plain": [
"42950/219"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\frac{1750}{73},\\,0,\\,\\frac{6490}{219},\\,0\\right)$$"
],
"text/plain": [
"(1750/73, 0, 6490/219, 0)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{}$$"
],
"text/plain": [
"salto"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{\\texttt{El dual es:}}$$"
],
"text/plain": [
"dual"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\begin{array}{l}\n",
"\\begin{array}{lcrcrcrcl}\n",
" \\max \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 700 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 900 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 10 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 15 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 19 y_{2} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 2 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 14 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 20 y_{2} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 5 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 3 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 16 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 15 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 6 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 5 \\\\\n",
" \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 12 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 8 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 4 \\\\\n",
"\\end{array} \\\\\n",
"y_{1}, y_{2}, y_{3} \\geq 0\n",
"\\end{array}$$"
],
"text/plain": [
"LP problem (use 'view(...)' or '%display typeset' for details)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\frac{42950}{219}$$"
],
"text/plain": [
"42950/219"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(0,\\,\\frac{47}{219},\\,\\frac{65}{219}\\right)$$"
],
"text/plain": [
"(0, 47/219, 65/219)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{}$$"
],
"text/plain": [
"salto"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{\\ \\ \\ Fase\\ I}$$"
],
"text/plain": [
"FI"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{red}{t}}_{1} & {\\color{red}{t}}_{2} & {\\color{red}{t}}_{3} & {b} \\\\\n",
"\\hline\n",
" {\\color{red}{t}}_{1} & 15 & 14 & 16 & 12 & -1 & 0 & 0 & 1 & 0 & 0 & 700 \\\\\n",
"{\\color{red}{t}}_{2} & 19 & 20 & 15 & 8 & 0 & -1 & 0 & 0 & 1 & 0 & 900 \\\\\n",
"{\\color{red}{t}}_{3} & -7 & -5 & 6 & 4 & 0 & 0 & -1 & 0 & 0 & 1 & 10 \\\\\n",
"\\hline\n",
" {W_{ind}} & -27 & -29 & -37 & -24 & 1 & 1 & 1 & 0 & 0 & 0 & 1610\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 t1 t2 t3| rhs]\n",
"[----+-------------------------------------------------+----]\n",
"[ t1| 15 14 16 12 -1 0 0 1 0 0| 700]\n",
"[ t2| 19 20 15 8 0 -1 0 0 1 0| 900]\n",
"[ t3| -7 -5 6 4 0 0 -1 0 0 1| 10]\n",
"[----+-------------------------------------------------+----]\n",
"[ W| -27 -29 -37 -24 1 1 1 0 0 0|1610]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" \n",
" \n",
"cambio( 3 , 3 )\n",
"cambio( 1 , 1 )\n",
"cambio( 1 , 2 )\n",
"cambio( 2 , 5 )\n"
]
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{red}{t}}_{1} & {\\color{red}{t}}_{2} & {\\color{red}{t}}_{3} & {b} \\\\\n",
"\\hline\n",
" {x}_{2} & \\frac{73}{65} & 1 & 0 & -\\frac{4}{65} & 0 & -\\frac{2}{65} & \\frac{1}{13} & 0 & \\frac{2}{65} & -\\frac{1}{13} & \\frac{350}{13} \\\\\n",
"{\\color{green}{s}}_{1} & -\\frac{193}{65} & 0 & 0 & -\\frac{196}{65} & 1 & -\\frac{164}{195} & -\\frac{22}{39} & -1 & \\frac{164}{195} & \\frac{22}{39} & \\frac{2440}{39} \\\\\n",
"{x}_{3} & -\\frac{3}{13} & 0 & 1 & \\frac{8}{13} & 0 & -\\frac{1}{39} & -\\frac{4}{39} & 0 & \\frac{1}{39} & \\frac{4}{39} & \\frac{940}{39} \\\\\n",
"\\hline\n",
" {W_{ind}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 1 & 1 & 1 & 0\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 t1 t2 t3| rhs]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ x2| 73/65 1 0 -4/65 0 -2/65 1/13 0 2/65 -1/13| 350/13]\n",
"[ s1| -193/65 0 0 -196/65 1 -164/195 -22/39 -1 164/195 22/39| 2440/39]\n",
"[ x3| -3/13 0 1 8/13 0 -1/39 -4/39 0 1/39 4/39| 940/39]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ W| 0 0 0 0 0 0 0 1 1 1| 0]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"El problema de la fase I es óptimo y Wopt=0. Por tanto, el original es factible\n",
"Preparemos el cuadro inicial de la fase II\n"
]
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{\\ \\ \\ Fase\\ II}$$"
],
"text/plain": [
"FII"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}{===================================}$$"
],
"text/plain": [
"ry"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {b} \\\\\n",
"\\hline\n",
" {x}_{2} & \\frac{73}{65} & 1 & 0 & -\\frac{4}{65} & 0 & -\\frac{2}{65} & \\frac{1}{13} & \\frac{350}{13} \\\\\n",
"{\\color{green}{s}}_{1} & -\\frac{193}{65} & 0 & 0 & -\\frac{196}{65} & 1 & -\\frac{164}{195} & -\\frac{22}{39} & \\frac{2440}{39} \\\\\n",
"{x}_{3} & -\\frac{3}{13} & 0 & 1 & \\frac{8}{13} & 0 & -\\frac{1}{39} & -\\frac{4}{39} & \\frac{940}{39} \\\\\n",
"\\hline\n",
" {Z_{ind}} & -\\frac{14}{65} & 0 & 0 & \\frac{72}{65} & 0 & \\frac{43}{195} & \\frac{11}{39} & \\frac{7850}{39}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3| rhs]\n",
"[--------+--------------------------------------------------------------+--------]\n",
"[ x2| 73/65 1 0 -4/65 0 -2/65 1/13| 350/13]\n",
"[ s1| -193/65 0 0 -196/65 1 -164/195 -22/39| 2440/39]\n",
"[ x3| -3/13 0 1 8/13 0 -1/39 -4/39| 940/39]\n",
"[--------+--------------------------------------------------------------+--------]\n",
"[ Z| -14/65 0 0 72/65 0 43/195 11/39| 7850/39]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"cambio( 1 , 1 )\n"
]
},
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & \\frac{65}{73} & 0 & -\\frac{4}{73} & 0 & -\\frac{2}{73} & \\frac{5}{73} & \\frac{1750}{73} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & \\frac{193}{73} & 0 & -\\frac{232}{73} & 1 & -\\frac{202}{219} & -\\frac{79}{219} & \\frac{29290}{219} \\\\\n",
"{x}_{3} & 0 & \\frac{15}{73} & 1 & \\frac{44}{73} & 0 & -\\frac{7}{219} & -\\frac{19}{219} & \\frac{6490}{219} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & \\frac{14}{73} & 0 & \\frac{80}{73} & 0 & \\frac{47}{219} & \\frac{65}{219} & \\frac{42950}{219}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3| rhs]\n",
"[---------+---------------------------------------------------------------------+---------]\n",
"[ x1| 1 65/73 0 -4/73 0 -2/73 5/73| 1750/73]\n",
"[ s1| 0 193/73 0 -232/73 1 -202/219 -79/219|29290/219]\n",
"[ x3| 0 15/73 1 44/73 0 -7/219 -19/219| 6490/219]\n",
"[---------+---------------------------------------------------------------------+---------]\n",
"[ Z| 0 14/73 0 80/73 0 47/219 65/219|42950/219]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"El problema es óptimo\n"
]
}
],
"source": [
"tipo='min' #max o min\n",
"n=4 # número de variables x's\n",
"m=3 # número de restricciones (distintas de las de signo)\n",
"\n",
"c=matrix(QQ,1,n,[2,3,5,4] ) # vector de costos\n",
"\n",
"A=matrix(QQ,m,n,[ # coeficientes de las variables x's\n",
"[15,14,16,12],\n",
"[19,20,15,8],\n",
"[-7,-5,6,4]\n",
"])\n",
"\n",
"b=matrix(QQ,m,1,[700,900,10]) # vector de la derecha\n",
"\n",
"\n",
"######################################################\n",
"############### VARIABLES DE HOLGURA ################\n",
"\n",
"\n",
"h=matrix(ZZ,1,m,[-1,-1,-1]) # variables de holgura: \n",
" # Pon 1 si va sumando, \n",
" # -1 si va restando y \n",
" # 0 si no hay variable de holgura.\n",
"\n",
"#######################################################\n",
"\n",
"\n",
"######################################################\n",
"######### Problemas Paramétricos ###########\n",
" \n",
"prmtr_b=0 # 1 si quieres usar un parámetro en b\n",
" # 0 en caso contrario\n",
"\n",
"bp=matrix(QQ,m,1,[1,-1,2]) # vector de la derecha para problemas paramétricos.\n",
"\n",
"###################################################### \n",
" \n",
"prmtr_c=0 # 1 si quieres usar un parámetro en c\n",
" # 0 en caso contrario\n",
"\n",
" \n",
"cp=matrix(QQ,1,n,[-1,1,0,2] ) # vector de costos adicional para problemas paramétricos.\n",
"\n",
"######################################################\n",
"\n",
"\n",
"##############################################################################\n",
"####### Usando el paquete de programación lineal Pinteractive de Sage ########\n",
"##############################################################################\n",
"signos=[]\n",
"for i in range(m):\n",
" if h[0,i] == -1:\n",
" signos.append(\">=\")\n",
" elif h[0,i] == 1:\n",
" signos.append(\"<=\")\n",
" elif h[0,i] == 0:\n",
" signos.append(\"==\")\n",
" \n",
"var_nam = []\n",
"for i in range(n):\n",
" k=i+1\n",
" var_nam.append('x_%s' %k)\n",
" \n",
"PPL = InteractiveLPProblem(A, b, c, var_nam, \n",
"constraint_type=signos, variable_type=\">=\", problem_type=\"min\") \n",
"\n",
"##############################################################################\n",
"##############################################################################\n",
"\n",
"\n",
"####################################################################\n",
"############ HASTA AQUÍ LOS DATOS ################################\n",
"####################################################################\n",
"\n",
"\n",
"raya=var('ry', latex_name='===================================')\n",
"show(ry)\n",
"show(ry)\n",
"show(PPL)\n",
"show(PPL.optimal_value())\n",
"show(PPL.optimal_solution())\n",
"show(ry)\n",
"show(ry)\n",
"\n",
"#load('https://sage.unex.es/501708/simplex2023Auto.sage')\n",
"load('https://sage.unex.es/501708/simplex2023Auto_MC.sage')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"71/73"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"1750/73 -23\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & \\frac{65}{73} & 0 & -\\frac{4}{73} & 0 & -\\frac{2}{73} & \\frac{5}{73} & 0 & \\frac{1750}{73} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & \\frac{193}{73} & 0 & -\\frac{232}{73} & 1 & -\\frac{202}{219} & -\\frac{79}{219} & 0 & \\frac{29290}{219} \\\\\n",
"{x}_{3} & 0 & \\frac{15}{73} & 1 & \\frac{44}{73} & 0 & -\\frac{7}{219} & -\\frac{19}{219} & 0 & \\frac{6490}{219} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & -65 & 0 & -69 & 0 & -71 & -5 & 1 & -71 \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & \\frac{14}{73} & 0 & \\frac{80}{73} & 0 & \\frac{47}{219} & \\frac{65}{219} & 0 & \\frac{42950}{219}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4| rhs]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ x1| 1 65/73 0 -4/73 0 -2/73 5/73 0| 1750/73]\n",
"[ s1| 0 193/73 0 -232/73 1 -202/219 -79/219 0|29290/219]\n",
"[ x3| 0 15/73 1 44/73 0 -7/219 -19/219 0| 6490/219]\n",
"[ s4| 0 -65 0 -69 0 -71 -5 1| -71]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 14/73 0 80/73 0 47/219 65/219 0|42950/219]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"v=matrix(QQ,1,nv,[0,-65,0,-69,0,-71,-5] ); R=-71;\n",
"añadeRes(v,R)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ NO -0.0029505 NO -0.015882 NO -0.0030227 -0.059361 NO]"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(4)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -1 & 0 & -1 & 0 & \\frac{1}{73} & 23 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{389}{65} & 1 & -\\frac{743}{195} & -\\frac{22}{39} & \\frac{193}{4745} & \\frac{25517}{195} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{13} & 0 & -\\frac{10}{39} & -\\frac{4}{39} & \\frac{3}{949} & \\frac{1147}{39} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{69}{65} & 0 & \\frac{71}{65} & \\frac{1}{13} & -\\frac{1}{65} & \\frac{71}{65} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{58}{65} & 0 & \\frac{1}{195} & \\frac{11}{39} & \\frac{14}{4745} & \\frac{38284}{195}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4| rhs]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ x1| 1 0 0 -1 0 -1 0 1/73| 23]\n",
"[ s1| 0 0 0 -389/65 1 -743/195 -22/39 193/4745|25517/195]\n",
"[ x3| 0 0 1 5/13 0 -10/39 -4/39 3/949| 1147/39]\n",
"[ x2| 0 1 0 69/65 0 71/65 1/13 -1/65| 71/65]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 0 0 58/65 0 1/195 11/39 14/4745|38284/195]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(4,2)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"16/39"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"1147/39 -29"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -1 & 0 & -1 & 0 & \\frac{1}{73} & 0 & 23 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{389}{65} & 1 & -\\frac{743}{195} & -\\frac{22}{39} & \\frac{193}{4745} & 0 & \\frac{25517}{195} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{13} & 0 & -\\frac{10}{39} & -\\frac{4}{39} & \\frac{3}{949} & 0 & \\frac{1147}{39} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{69}{65} & 0 & \\frac{71}{65} & \\frac{1}{13} & -\\frac{1}{65} & 0 & \\frac{71}{65} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & -4 & 0 & -6 & -5 & -\\frac{6}{59} & 1 & -6 \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{58}{65} & 0 & \\frac{1}{195} & \\frac{11}{39} & \\frac{14}{4745} & 0 & \\frac{38284}{195}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5| rhs]\n",
"[---------+-----------------------------------------------------------------------------------------+---------]\n",
"[ x1| 1 0 0 -1 0 -1 0 1/73 0| 23]\n",
"[ s1| 0 0 0 -389/65 1 -743/195 -22/39 193/4745 0|25517/195]\n",
"[ x3| 0 0 1 5/13 0 -10/39 -4/39 3/949 0| 1147/39]\n",
"[ x2| 0 1 0 69/65 0 71/65 1/13 -1/65 0| 71/65]\n",
"[ s5| 0 0 0 -4 0 -6 -5 -6/59 1| -6]\n",
"[---------+-----------------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 0 0 58/65 0 1/195 11/39 14/4745 0|38284/195]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"v=matrix(QQ,1,nv,[0,0,0,-4,0,-6,-65/13,-6/59] ); R=-6;\n",
"añadeRes(v,R)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -0.22308 NO -0.00085470 -0.056410 -0.029013 NO]"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(5)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{1}{3} & 0 & 0 & \\frac{5}{6} & \\frac{132}{4307} & -\\frac{1}{6} & 24 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{31}{9} & 1 & 0 & \\frac{47}{18} & \\frac{1360}{12921} & -\\frac{743}{1170} & \\frac{404}{3} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{9} & 0 & 0 & \\frac{1}{9} & \\frac{97}{12921} & -\\frac{5}{117} & \\frac{89}{3} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{1}{3} & 0 & 0 & -\\frac{5}{6} & -\\frac{2}{59} & \\frac{71}{390} & 0 \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & \\frac{2}{3} & 0 & 1 & \\frac{5}{6} & \\frac{1}{59} & -\\frac{1}{6} & 1 \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{8}{9} & 0 & 0 & \\frac{5}{18} & \\frac{37}{12921} & \\frac{1}{1170} & \\frac{589}{3}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5| rhs]\n",
"[----------+--------------------------------------------------------------------------------------------------+----------]\n",
"[ x1| 1 0 0 -1/3 0 0 5/6 132/4307 -1/6| 24]\n",
"[ s1| 0 0 0 -31/9 1 0 47/18 1360/12921 -743/1170| 404/3]\n",
"[ x3| 0 0 1 5/9 0 0 1/9 97/12921 -5/117| 89/3]\n",
"[ x2| 0 1 0 1/3 0 0 -5/6 -2/59 71/390| 0]\n",
"[ s2| 0 0 0 2/3 0 1 5/6 1/59 -1/6| 1]\n",
"[----------+--------------------------------------------------------------------------------------------------+----------]\n",
"[ Z| 0 0 0 8/9 0 0 5/18 37/12921 1/1170| 589/3]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(5,6)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"39"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"117/3"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{1}{3} & 0 & 0 & \\frac{5}{6} & \\frac{132}{4307} & -\\frac{1}{6} & 0 & 24 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{31}{9} & 1 & 0 & \\frac{47}{18} & \\frac{1360}{12921} & -\\frac{743}{1170} & 0 & \\frac{404}{3} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{9} & 0 & 0 & \\frac{1}{9} & \\frac{97}{12921} & -\\frac{5}{117} & 0 & \\frac{89}{3} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{1}{3} & 0 & 0 & -\\frac{5}{6} & -\\frac{2}{59} & \\frac{71}{390} & 0 & 0 \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & \\frac{2}{3} & 0 & 1 & \\frac{5}{6} & \\frac{1}{59} & -\\frac{1}{6} & 0 & 1 \\\\\n",
"{\\color{green}{s}}_{6} & 0 & 0 & 0 & -5 & 0 & 0 & -1 & -\\frac{291}{4307} & -\\frac{10}{13} & 1 & -6 \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{8}{9} & 0 & 0 & \\frac{5}{18} & \\frac{37}{12921} & \\frac{1}{1170} & 0 & \\frac{589}{3}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6| rhs]\n",
"[----------+-------------------------------------------------------------------------------------------------------------+----------]\n",
"[ x1| 1 0 0 -1/3 0 0 5/6 132/4307 -1/6 0| 24]\n",
"[ s1| 0 0 0 -31/9 1 0 47/18 1360/12921 -743/1170 0| 404/3]\n",
"[ x3| 0 0 1 5/9 0 0 1/9 97/12921 -5/117 0| 89/3]\n",
"[ x2| 0 1 0 1/3 0 0 -5/6 -2/59 71/390 0| 0]\n",
"[ s2| 0 0 0 2/3 0 1 5/6 1/59 -1/6 0| 1]\n",
"[ s6| 0 0 0 -5 0 0 -1 -291/4307 -10/13 1| -6]\n",
"[----------+-------------------------------------------------------------------------------------------------------------+----------]\n",
"[ Z| 0 0 0 8/9 0 0 5/18 37/12921 1/1170 0| 589/3]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"v=matrix(QQ,1,nv,[0,0,0,-5,0,0,-1,-(97*3)/4307,-10/13] ); R=-6;\n",
"añadeRes(v,R)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -0.17778 NO NO -0.27778 -0.042383 -0.0011111 NO]"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(6)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{9}{2} & 0 & 0 & 0 & -\\frac{221}{8614} & -\\frac{21}{26} & \\frac{5}{6} & 19 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{33}{2} & 1 & 0 & 0 & -\\frac{613}{8614} & -\\frac{1031}{390} & \\frac{47}{18} & 119 \\\\\n",
"{x}_{3} & 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & -\\frac{5}{39} & \\frac{1}{9} & 29 \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{9}{2} & 0 & 0 & 0 & \\frac{193}{8614} & \\frac{107}{130} & -\\frac{5}{6} & 5 \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & -\\frac{7}{2} & 0 & 1 & 0 & -\\frac{339}{8614} & -\\frac{21}{26} & \\frac{5}{6} & -4 \\\\\n",
"{\\color{green}{s}}_{3} & 0 & 0 & 0 & 5 & 0 & 0 & 1 & \\frac{291}{4307} & \\frac{10}{13} & -1 & 6 \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & -\\frac{1}{2} & 0 & 0 & 0 & -\\frac{137}{8614} & -\\frac{83}{390} & \\frac{5}{18} & 198\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6| rhs]\n",
"[---------+---------------------------------------------------------------------------------------------------+---------]\n",
"[ x1| 1 0 0 -9/2 0 0 0 -221/8614 -21/26 5/6| 19]\n",
"[ s1| 0 0 0 -33/2 1 0 0 -613/8614 -1031/390 47/18| 119]\n",
"[ x3| 0 0 1 0 0 0 0 0 -5/39 1/9| 29]\n",
"[ x2| 0 1 0 9/2 0 0 0 193/8614 107/130 -5/6| 5]\n",
"[ s2| 0 0 0 -7/2 0 1 0 -339/8614 -21/26 5/6| -4]\n",
"[ s3| 0 0 0 5 0 0 1 291/4307 10/13 -1| 6]\n",
"[---------+---------------------------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 0 0 -1/2 0 0 0 -137/8614 -83/390 5/18| 198]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(6,7)"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"8/15"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"443/15 -29"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {\\color{green}{s}}_{7} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{9}{2} & 0 & 0 & 0 & -\\frac{221}{8614} & -\\frac{21}{26} & \\frac{5}{6} & 0 & 19 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{33}{2} & 1 & 0 & 0 & -\\frac{613}{8614} & -\\frac{1031}{390} & \\frac{47}{18} & 0 & 119 \\\\\n",
"{x}_{3} & 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & -\\frac{5}{39} & \\frac{1}{9} & 0 & 29 \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{9}{2} & 0 & 0 & 0 & \\frac{193}{8614} & \\frac{107}{130} & -\\frac{5}{6} & 0 & 5 \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & -\\frac{7}{2} & 0 & 1 & 0 & -\\frac{339}{8614} & -\\frac{21}{26} & \\frac{5}{6} & 0 & -4 \\\\\n",
"{\\color{green}{s}}_{3} & 0 & 0 & 0 & 5 & 0 & 0 & 1 & \\frac{291}{4307} & \\frac{10}{13} & -1 & 0 & 6 \\\\\n",
"{\\color{green}{s}}_{7} & 0 & 0 & 0 & -\\frac{7}{15} & 0 & 0 & 0 & -\\frac{137}{21535} & -\\frac{16}{325} & -\\frac{2}{15} & 1 & -\\frac{8}{15} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & -\\frac{1}{2} & 0 & 0 & 0 & -\\frac{137}{8614} & -\\frac{83}{390} & \\frac{5}{18} & 0 & 198\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7| rhs]\n",
"[----------+------------------------------------------------------------------------------------------------------------------------+----------]\n",
"[ x1| 1 0 0 -9/2 0 0 0 -221/8614 -21/26 5/6 0| 19]\n",
"[ s1| 0 0 0 -33/2 1 0 0 -613/8614 -1031/390 47/18 0| 119]\n",
"[ x3| 0 0 1 0 0 0 0 0 -5/39 1/9 0| 29]\n",
"[ x2| 0 1 0 9/2 0 0 0 193/8614 107/130 -5/6 0| 5]\n",
"[ s2| 0 0 0 -7/2 0 1 0 -339/8614 -21/26 5/6 0| -4]\n",
"[ s3| 0 0 0 5 0 0 1 291/4307 10/13 -1 0| 6]\n",
"[ s7| 0 0 0 -7/15 0 0 0 -137/21535 -16/325 -2/15 1| -8/15]\n",
"[----------+------------------------------------------------------------------------------------------------------------------------+----------]\n",
"[ Z| 0 0 0 -1/2 0 0 0 -137/8614 -83/390 5/18 0| 198]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"v=matrix(QQ,1,nv,[0,0,0,-7/15,0,0,0,-137/21535,-16/325,-2/15] ); R=-8/15;\n",
"añadeRes(v,R)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO 1.0714 NO NO NO 2.5000 4.3229 -2.0833 NO]"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(7)"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {\\color{green}{s}}_{7} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{1076}{411} & 0 & 0 & 0 & 0 & -\\frac{10849}{17810} & \\frac{1127}{822} & -\\frac{1105}{274} & \\frac{8693}{411} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{4636}{411} & 1 & 0 & 0 & 0 & -\\frac{111823}{53430} & \\frac{10117}{2466} & -\\frac{3065}{274} & \\frac{51361}{411} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & -\\frac{5}{39} & \\frac{1}{9} & 0 & 29 \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{1174}{411} & 0 & 0 & 0 & 0 & \\frac{11571}{17810} & -\\frac{357}{274} & \\frac{965}{274} & \\frac{1283}{411} \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & -\\frac{84}{137} & 0 & 1 & 0 & 0 & -\\frac{8961}{17810} & \\frac{1363}{822} & -\\frac{1695}{274} & -\\frac{96}{137} \\\\\n",
"{\\color{green}{s}}_{3} & 0 & 0 & 0 & \\frac{6}{137} & 0 & 0 & 1 & 0 & \\frac{2194}{8905} & -\\frac{331}{137} & \\frac{1455}{137} & \\frac{46}{137} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{30149}{411} & 0 & 0 & 0 & 1 & \\frac{68912}{8905} & \\frac{8614}{411} & -\\frac{21535}{137} & \\frac{34456}{411} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{2}{3} & 0 & 0 & 0 & 0 & -\\frac{7}{78} & \\frac{11}{18} & -\\frac{5}{2} & \\frac{590}{3}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7| rhs]\n",
"[-------------+---------------------------------------------------------------------------------------------------------------------------------------------------------+-------------]\n",
"[ x1| 1 0 0 -1076/411 0 0 0 0 -10849/17810 1127/822 -1105/274| 8693/411]\n",
"[ s1| 0 0 0 -4636/411 1 0 0 0 -111823/53430 10117/2466 -3065/274| 51361/411]\n",
"[ x3| 0 0 1 0 0 0 0 0 -5/39 1/9 0| 29]\n",
"[ x2| 0 1 0 1174/411 0 0 0 0 11571/17810 -357/274 965/274| 1283/411]\n",
"[ s2| 0 0 0 -84/137 0 1 0 0 -8961/17810 1363/822 -1695/274| -96/137]\n",
"[ s3| 0 0 0 6/137 0 0 1 0 2194/8905 -331/137 1455/137| 46/137]\n",
"[ s4| 0 0 0 30149/411 0 0 0 1 68912/8905 8614/411 -21535/137| 34456/411]\n",
"[-------------+---------------------------------------------------------------------------------------------------------------------------------------------------------+-------------]\n",
"[ Z| 0 0 0 2/3 0 0 0 0 -7/78 11/18 -5/2| 590/3]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(7,8)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -1.0873 NO NO NO NO 0.17837 NO 0.40413]"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(5)"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\begin{array}{r|rrrrrrrrrrr|r}\n",
"{Basic} & {x}_{1} & {x}_{2} & {x}_{3} & {x}_{4} & {\\color{green}{s}}_{1} & {\\color{green}{s}}_{2} & {\\color{green}{s}}_{3} & {\\color{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {\\color{green}{s}}_{7} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{16808}{8961} & 0 & -\\frac{10849}{8961} & 0 & 0 & 0 & -\\frac{590}{927} & \\frac{10325}{2987} & \\frac{197135}{8961} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{78224}{8961} & 1 & -\\frac{111823}{26883} & 0 & 0 & 0 & -\\frac{7772}{2781} & \\frac{130345}{8961} & \\frac{1145939}{8961} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{1400}{8961} & 0 & -\\frac{6850}{26883} & 0 & 0 & 0 & -\\frac{866}{2781} & \\frac{14125}{8961} & \\frac{261469}{8961} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{638}{309} & 0 & \\frac{133}{103} & 0 & 0 & 0 & \\frac{259}{309} & -\\frac{460}{103} & \\frac{685}{309} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & \\frac{3640}{2987} & 0 & -\\frac{17810}{8961} & 0 & 0 & 1 & -\\frac{3055}{927} & \\frac{36725}{2987} & \\frac{4160}{2987} \\\\\n",
"{\\color{green}{s}}_{3} & 0 & 0 & 0 & -\\frac{766}{2987} & 0 & \\frac{4388}{8961} & 1 & 0 & 0 & -\\frac{1487}{927} & \\frac{22675}{2987} & -\\frac{22}{2987} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{572831}{8961} & 0 & \\frac{137824}{8961} & 0 & 1 & 0 & \\frac{43070}{927} & -\\frac{753725}{2987} & \\frac{654664}{8961} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{2318}{2987} & 0 & -\\frac{4795}{26883} & 0 & 0 & 0 & \\frac{877}{2781} & -\\frac{12515}{8961} & \\frac{587070}{2987}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7| rhs]\n",
"[-------------+---------------------------------------------------------------------------------------------------------------------------------------------------------+-------------]\n",
"[ x1| 1 0 0 -16808/8961 0 -10849/8961 0 0 0 -590/927 10325/2987| 197135/8961]\n",
"[ s1| 0 0 0 -78224/8961 1 -111823/26883 0 0 0 -7772/2781 130345/8961| 1145939/8961]\n",
"[ x3| 0 0 1 1400/8961 0 -6850/26883 0 0 0 -866/2781 14125/8961| 261469/8961]\n",
"[ x2| 0 1 0 638/309 0 133/103 0 0 0 259/309 -460/103| 685/309]\n",
"[ s5| 0 0 0 3640/2987 0 -17810/8961 0 0 1 -3055/927 36725/2987| 4160/2987]\n",
"[ s3| 0 0 0 -766/2987 0 4388/8961 1 0 0 -1487/927 22675/2987| -22/2987]\n",
"[ s4| 0 0 0 572831/8961 0 137824/8961 0 1 0 43070/927 -753725/2987| 654664/8961]\n",
"[-------------+---------------------------------------------------------------------------------------------------------------------------------------------------------+-------------]\n",
"[ Z| 0 0 0 2318/2987 0 -4795/26883 0 0 0 877/2781 -12515/8961| 587070/2987]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(5,9)"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"[ 197135/8961]\n",
"[1145939/8961]\n",
"[ 261469/8961]\n",
"[ 685/309]\n",
"[ 4160/2987]\n",
"[ -22/2987]\n",
"[ 654664/8961]"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"b # Vector de la derecha del cuadro actual"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Para obtener la parte decimal de $b_1$"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"8954/8961"
]
},
"execution_count": 22,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"b[0][0]- floor(b[0][0])"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"8954/8961"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"b[0][0]- floor(b[0][0])"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [
{
"data": {
"text/plain": [
"1114/8961"
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"A[0][3]-floor(A[0][3])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "SageMath 9.5",
"language": "sage",
"name": "sagemath"
},
"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.11.2"
}
},
"nbformat": 4,
"nbformat_minor": 4
}