{
"cells": [
{
"cell_type": "code",
"execution_count": 102,
"metadata": {
"jupyter": {
"outputs_hidden": true
},
"scrolled": false,
"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": [
"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"
]
},
{
"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} & \\frac{101}{3} & \\frac{82}{3} & 0 & \\frac{4}{3} & -1 & 0 & \\frac{8}{3} & 1 & 0 & -\\frac{8}{3} & \\frac{2020}{3} \\\\\n",
"{\\color{red}{t}}_{2} & \\frac{73}{2} & \\frac{65}{2} & 0 & -2 & 0 & -1 & \\frac{5}{2} & 0 & 1 & -\\frac{5}{2} & 875 \\\\\n",
"{x}_{3} & -\\frac{7}{6} & -\\frac{5}{6} & 1 & \\frac{2}{3} & 0 & 0 & -\\frac{1}{6} & 0 & 0 & \\frac{1}{6} & \\frac{5}{3} \\\\\n",
"\\hline\n",
" {W_{ind}} & -\\frac{421}{6} & -\\frac{359}{6} & 0 & \\frac{2}{3} & 1 & 1 & -\\frac{31}{6} & 0 & 0 & \\frac{37}{6} & \\frac{4645}{3}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 t1 t2 t3| rhs]\n",
"[------+---------------------------------------------------------------------+------]\n",
"[ t1| 101/3 82/3 0 4/3 -1 0 8/3 1 0 -8/3|2020/3]\n",
"[ t2| 73/2 65/2 0 -2 0 -1 5/2 0 1 -5/2| 875]\n",
"[ x3| -7/6 -5/6 1 2/3 0 0 -1/6 0 0 1/6| 5/3]\n",
"[------+---------------------------------------------------------------------+------]\n",
"[ W|-421/6 -359/6 0 2/3 1 1 -31/6 0 0 37/6|4645/3]"
]
},
"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|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}_{1} & 1 & \\frac{82}{101} & 0 & \\frac{4}{101} & -\\frac{3}{101} & 0 & \\frac{8}{101} & \\frac{3}{101} & 0 & -\\frac{8}{101} & 20 \\\\\n",
"{\\color{red}{t}}_{2} & 0 & \\frac{579}{202} & 0 & -\\frac{348}{101} & \\frac{219}{202} & -1 & -\\frac{79}{202} & -\\frac{219}{202} & 1 & \\frac{79}{202} & 145 \\\\\n",
"{x}_{3} & 0 & \\frac{23}{202} & 1 & \\frac{72}{101} & -\\frac{7}{202} & 0 & -\\frac{15}{202} & \\frac{7}{202} & 0 & \\frac{15}{202} & 25 \\\\\n",
"\\hline\n",
" {W_{ind}} & 0 & -\\frac{579}{202} & 0 & \\frac{348}{101} & -\\frac{219}{202} & 1 & \\frac{79}{202} & \\frac{421}{202} & 0 & \\frac{123}{202} & 145\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 t1 t2 t3| rhs]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ x1| 1 82/101 0 4/101 -3/101 0 8/101 3/101 0 -8/101| 20]\n",
"[ t2| 0 579/202 0 -348/101 219/202 -1 -79/202 -219/202 1 79/202| 145]\n",
"[ x3| 0 23/202 1 72/101 -7/202 0 -15/202 7/202 0 15/202| 25]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ W| 0 -579/202 0 348/101 -219/202 1 79/202 421/202 0 123/202| 145]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"cambio( 1 , 2 )\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{101}{82} & 1 & 0 & \\frac{2}{41} & -\\frac{3}{82} & 0 & \\frac{4}{41} & \\frac{3}{82} & 0 & -\\frac{4}{41} & \\frac{1010}{41} \\\\\n",
"{\\color{red}{t}}_{2} & -\\frac{579}{164} & 0 & 0 & -\\frac{147}{41} & \\frac{195}{164} & -1 & -\\frac{55}{82} & -\\frac{195}{164} & 1 & \\frac{55}{82} & \\frac{3050}{41} \\\\\n",
"{x}_{3} & -\\frac{23}{164} & 0 & 1 & \\frac{29}{41} & -\\frac{5}{164} & 0 & -\\frac{7}{82} & \\frac{5}{164} & 0 & \\frac{7}{82} & \\frac{910}{41} \\\\\n",
"\\hline\n",
" {W_{ind}} & \\frac{579}{164} & 0 & 0 & \\frac{147}{41} & -\\frac{195}{164} & 1 & \\frac{55}{82} & \\frac{359}{164} & 0 & \\frac{27}{82} & \\frac{3050}{41}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 t1 t2 t3| rhs]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ x2| 101/82 1 0 2/41 -3/82 0 4/41 3/82 0 -4/41| 1010/41]\n",
"[ t2|-579/164 0 0 -147/41 195/164 -1 -55/82 -195/164 1 55/82| 3050/41]\n",
"[ x3| -23/164 0 1 29/41 -5/164 0 -7/82 5/164 0 7/82| 910/41]\n",
"[--------+-----------------------------------------------------------------------------------------+--------]\n",
"[ W| 579/164 0 0 147/41 -195/164 1 55/82 359/164 0 27/82| 3050/41]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"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"
},
{
"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"
},
{
"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",
"\n",
"\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": 101,
"metadata": {
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [],
"source": [
"#v=matrix(QQ,1,nv,[0,-193/73,0,-(232/73)*(163/56),0,-202*163/(219*56),-(79/219)*(163/56)] ); R=-163/219;\n",
"#añadeRes(v,R) # hecho sin mejoras"
]
},
{
"cell_type": "code",
"execution_count": 103,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"[71/73, 163/219, 139/219]"
]
},
"execution_count": 103,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pardec()"
]
},
{
"cell_type": "code",
"execution_count": 104,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"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 & -\\frac{65}{73} & 0 & -\\frac{69}{73} & 0 & -\\frac{71}{73} & -\\frac{5}{73} & 1 & -\\frac{71}{73} \\\\\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/73 0 -69/73 0 -71/73 -5/73 1| -71/73]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 14/73 0 80/73 0 47/219 65/219 0|42950/219]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"corte(1)"
]
},
{
"cell_type": "code",
"execution_count": 105,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[ NO -0.21538 NO -1.1594 NO -0.22066 -4.3333 NO]"
]
},
"execution_count": 105,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(4)"
]
},
{
"cell_type": "code",
"execution_count": 106,
"metadata": {},
"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 & 1 & 23 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{389}{65} & 1 & -\\frac{743}{195} & -\\frac{22}{39} & \\frac{193}{65} & \\frac{25517}{195} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{13} & 0 & -\\frac{10}{39} & -\\frac{4}{39} & \\frac{3}{13} & \\frac{1147}{39} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{69}{65} & 0 & \\frac{71}{65} & \\frac{1}{13} & -\\frac{73}{65} & \\frac{71}{65} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{58}{65} & 0 & \\frac{1}{195} & \\frac{11}{39} & \\frac{14}{65} & \\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| 23]\n",
"[ s1| 0 0 0 -389/65 1 -743/195 -22/39 193/65|25517/195]\n",
"[ x3| 0 0 1 5/13 0 -10/39 -4/39 3/13| 1147/39]\n",
"[ x2| 0 1 0 69/65 0 71/65 1/13 -73/65| 71/65]\n",
"[---------+-------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 0 0 58/65 0 1/195 11/39 14/65|38284/195]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(4,2)"
]
},
{
"cell_type": "code",
"execution_count": 107,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"[0, 167/195, 16/39, 6/65]"
]
},
"execution_count": 107,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pardec()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"193"
]
},
{
"cell_type": "code",
"execution_count": 83,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"292/615"
]
},
"execution_count": 83,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"292/(41*15)"
]
},
{
"cell_type": "code",
"execution_count": 108,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"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 & 1 & 0 & 23 \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{389}{65} & 1 & -\\frac{743}{195} & -\\frac{22}{39} & \\frac{193}{65} & 0 & \\frac{25517}{195} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5}{13} & 0 & -\\frac{10}{39} & -\\frac{4}{39} & \\frac{3}{13} & 0 & \\frac{1147}{39} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{69}{65} & 0 & \\frac{71}{65} & \\frac{1}{13} & -\\frac{73}{65} & 0 & \\frac{71}{65} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & -\\frac{1}{65} & 0 & -\\frac{37}{195} & -\\frac{17}{39} & -\\frac{193}{65} & 1 & -\\frac{167}{195} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{58}{65} & 0 & \\frac{1}{195} & \\frac{11}{39} & \\frac{14}{65} & 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 0| 23]\n",
"[ s1| 0 0 0 -389/65 1 -743/195 -22/39 193/65 0|25517/195]\n",
"[ x3| 0 0 1 5/13 0 -10/39 -4/39 3/13 0| 1147/39]\n",
"[ x2| 0 1 0 69/65 0 71/65 1/13 -73/65 0| 71/65]\n",
"[ s5| 0 0 0 -1/65 0 -37/195 -17/39 -193/65 1| -167/195]\n",
"[---------+-----------------------------------------------------------------------------------------+---------]\n",
"[ Z| 0 0 0 58/65 0 1/195 11/39 14/65 0|38284/195]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"corte(2)"
]
},
{
"cell_type": "code",
"execution_count": 85,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -82.000 NO -2.2864 -6.0068 -1.9658 NO]"
]
},
"execution_count": 85,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(5)"
]
},
{
"cell_type": "code",
"execution_count": 86,
"metadata": {
"scrolled": false
},
"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{199}{73} & 0 & -\\frac{67}{73} & -\\frac{60}{73} & 0 & \\frac{2665}{292} & \\frac{1490}{73} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{811}{73} & 1 & -\\frac{781}{219} & -\\frac{658}{219} & 0 & \\frac{7913}{292} & \\frac{26974}{219} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & -\\frac{1}{73} & 0 & -\\frac{52}{219} & -\\frac{64}{219} & 0 & \\frac{615}{292} & \\frac{6310}{219} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & 3 & 0 & 1 & 1 & 0 & -\\frac{41}{4} & 4 \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{1}{73} & 0 & \\frac{2627}{17520} & \\frac{959}{17520} & 1 & -\\frac{615}{292} & \\frac{41}{219} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{38}{73} & 0 & \\frac{5}{219} & \\frac{23}{219} & 0 & \\frac{287}{146} & \\frac{43118}{219}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5| rhs]\n",
"[----------+--------------------------------------------------------------------------------------------------+----------]\n",
"[ x1| 1 0 0 -199/73 0 -67/73 -60/73 0 2665/292| 1490/73]\n",
"[ s1| 0 0 0 -811/73 1 -781/219 -658/219 0 7913/292| 26974/219]\n",
"[ x3| 0 0 1 -1/73 0 -52/219 -64/219 0 615/292| 6310/219]\n",
"[ x2| 0 1 0 3 0 1 1 0 -41/4| 4]\n",
"[ s4| 0 0 0 1/73 0 2627/17520 959/17520 1 -615/292| 41/219]\n",
"[----------+--------------------------------------------------------------------------------------------------+----------]\n",
"[ Z| 0 0 0 38/73 0 5/219 23/219 0 287/146| 43118/219]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(5,8)"
]
},
{
"cell_type": "code",
"execution_count": 87,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"[30/73, 37/219, 178/219, 0, 41/219]"
]
},
"execution_count": 87,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pardec()"
]
},
{
"cell_type": "code",
"execution_count": 88,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"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{green}{s}}_{4} & {\\color{green}{s}}_{5} & {\\color{green}{s}}_{6} & {b} \\\\\n",
"\\hline\n",
" {x}_{1} & 1 & 0 & 0 & -\\frac{199}{73} & 0 & -\\frac{67}{73} & -\\frac{60}{73} & 0 & \\frac{2665}{292} & 0 & \\frac{1490}{73} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{811}{73} & 1 & -\\frac{781}{219} & -\\frac{658}{219} & 0 & \\frac{7913}{292} & 0 & \\frac{26974}{219} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & -\\frac{1}{73} & 0 & -\\frac{52}{219} & -\\frac{64}{219} & 0 & \\frac{615}{292} & 0 & \\frac{6310}{219} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & 3 & 0 & 1 & 1 & 0 & -\\frac{41}{4} & 0 & 4 \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{1}{73} & 0 & \\frac{2627}{17520} & \\frac{959}{17520} & 1 & -\\frac{615}{292} & 0 & \\frac{41}{219} \\\\\n",
"{\\color{green}{s}}_{6} & 0 & 0 & 0 & -\\frac{148}{6643} & 0 & -\\frac{2294}{19929} & -\\frac{37}{39858} & 0 & -\\frac{7913}{292} & 1 & -\\frac{37}{219} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{38}{73} & 0 & \\frac{5}{219} & \\frac{23}{219} & 0 & \\frac{287}{146} & 0 & \\frac{43118}{219}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6| rhs]\n",
"[-----------+-----------------------------------------------------------------------------------------------------------------------+-----------]\n",
"[ x1| 1 0 0 -199/73 0 -67/73 -60/73 0 2665/292 0| 1490/73]\n",
"[ s1| 0 0 0 -811/73 1 -781/219 -658/219 0 7913/292 0| 26974/219]\n",
"[ x3| 0 0 1 -1/73 0 -52/219 -64/219 0 615/292 0| 6310/219]\n",
"[ x2| 0 1 0 3 0 1 1 0 -41/4 0| 4]\n",
"[ s4| 0 0 0 1/73 0 2627/17520 959/17520 1 -615/292 0| 41/219]\n",
"[ s6| 0 0 0 -148/6643 0 -2294/19929 -37/39858 0 -7913/292 1| -37/219]\n",
"[-----------+-----------------------------------------------------------------------------------------------------------------------+-----------]\n",
"[ Z| 0 0 0 38/73 0 5/219 23/219 0 287/146 0| 43118/219]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"corte(2)"
]
},
{
"cell_type": "code",
"execution_count": 89,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -23.365 NO -0.19834 -113.14 NO -0.072539 NO]"
]
},
"execution_count": 89,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(6)"
]
},
{
"cell_type": "code",
"execution_count": 90,
"metadata": {},
"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{3693}{1351} & 0 & -\\frac{3877}{4053} & -\\frac{6665}{8106} & 0 & 0 & \\frac{65}{193} & \\frac{11785}{579} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{1013}{91} & 1 & -\\frac{335}{91} & -\\frac{547}{182} & 0 & 0 & 1 & 123 \\\\\n",
"{x}_{3} & 0 & 0 & 1 & -\\frac{271}{17563} & 0 & -\\frac{12982}{52689} & -\\frac{30803}{105378} & 0 & 0 & \\frac{15}{193} & \\frac{16675}{579} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{52837}{17563} & 0 & \\frac{54983}{52689} & \\frac{105415}{105378} & 0 & 0 & -\\frac{73}{193} & \\frac{2353}{579} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{271}{17563} & 0 & \\frac{669737}{4215120} & \\frac{231029}{4215120} & 1 & 0 & -\\frac{15}{193} & \\frac{116}{579} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & \\frac{592}{720083} & 0 & \\frac{9176}{2160249} & \\frac{74}{2160249} & 0 & 1 & -\\frac{292}{7913} & \\frac{148}{23739} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{1302}{2509} & 0 & \\frac{109}{7527} & \\frac{790}{7527} & 0 & 0 & \\frac{14}{193} & \\frac{114004}{579}\n",
"\\end{array}\\right)$$"
],
"text/plain": [
"[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6| rhs]\n",
"[--------------+-----------------------------------------------------------------------------------------------------------------------------------------------------+--------------]\n",
"[ x1| 1 0 0 -3693/1351 0 -3877/4053 -6665/8106 0 0 65/193| 11785/579]\n",
"[ s1| 0 0 0 -1013/91 1 -335/91 -547/182 0 0 1| 123]\n",
"[ x3| 0 0 1 -271/17563 0 -12982/52689 -30803/105378 0 0 15/193| 16675/579]\n",
"[ x2| 0 1 0 52837/17563 0 54983/52689 105415/105378 0 0 -73/193| 2353/579]\n",
"[ s4| 0 0 0 271/17563 0 669737/4215120 231029/4215120 1 0 -15/193| 116/579]\n",
"[ s5| 0 0 0 592/720083 0 9176/2160249 74/2160249 0 1 -292/7913| 148/23739]\n",
"[--------------+-----------------------------------------------------------------------------------------------------------------------------------------------------+--------------]\n",
"[ Z| 0 0 0 1302/2509 0 109/7527 790/7527 0 0 14/193| 114004/579]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(6,9)"
]
},
{
"cell_type": "code",
"execution_count": 91,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"[205/579, 0, 463/579, 37/579, 116/579, 148/23739]"
]
},
"execution_count": 91,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pardec()"
]
},
{
"cell_type": "code",
"execution_count": 92,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/latex": [
"$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left[1, 2, 3, 4, 5, 6, 7\\right]$$"
],
"text/plain": [
"[1, 2, 3, 4, 5, 6, 7]"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"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{3693}{1351} & 0 & -\\frac{3877}{4053} & -\\frac{6665}{8106} & 0 & 0 & \\frac{65}{193} & 0 & \\frac{11785}{579} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{1013}{91} & 1 & -\\frac{335}{91} & -\\frac{547}{182} & 0 & 0 & 1 & 0 & 123 \\\\\n",
"{x}_{3} & 0 & 0 & 1 & -\\frac{271}{17563} & 0 & -\\frac{12982}{52689} & -\\frac{30803}{105378} & 0 & 0 & \\frac{15}{193} & 0 & \\frac{16675}{579} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{52837}{17563} & 0 & \\frac{54983}{52689} & \\frac{105415}{105378} & 0 & 0 & -\\frac{73}{193} & 0 & \\frac{2353}{579} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & \\frac{271}{17563} & 0 & \\frac{669737}{4215120} & \\frac{231029}{4215120} & 1 & 0 & -\\frac{15}{193} & 0 & \\frac{116}{579} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & \\frac{592}{720083} & 0 & \\frac{9176}{2160249} & \\frac{74}{2160249} & 0 & 1 & -\\frac{292}{7913} & 0 & \\frac{148}{23739} \\\\\n",
"{\\color{green}{s}}_{7} & 0 & 0 & 0 & -\\frac{148}{17563} & 0 & -\\frac{2294}{52689} & -\\frac{37}{105378} & 0 & 0 & -\\frac{2701}{104606} & 1 & -\\frac{37}{579} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{1302}{2509} & 0 & \\frac{109}{7527} & \\frac{790}{7527} & 0 & 0 & \\frac{14}{193} & 0 & \\frac{114004}{579}\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 -3693/1351 0 -3877/4053 -6665/8106 0 0 65/193 0| 11785/579]\n",
"[ s1| 0 0 0 -1013/91 1 -335/91 -547/182 0 0 1 0| 123]\n",
"[ x3| 0 0 1 -271/17563 0 -12982/52689 -30803/105378 0 0 15/193 0| 16675/579]\n",
"[ x2| 0 1 0 52837/17563 0 54983/52689 105415/105378 0 0 -73/193 0| 2353/579]\n",
"[ s4| 0 0 0 271/17563 0 669737/4215120 231029/4215120 1 0 -15/193 0| 116/579]\n",
"[ s5| 0 0 0 592/720083 0 9176/2160249 74/2160249 0 1 -292/7913 0| 148/23739]\n",
"[ s7| 0 0 0 -148/17563 0 -2294/52689 -37/105378 0 0 -2701/104606 1| -37/579]\n",
"[--------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------+--------------]\n",
"[ Z| 0 0 0 1302/2509 0 109/7527 790/7527 0 0 14/193 0| 114004/579]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"corte(4)"
]
},
{
"cell_type": "code",
"execution_count": 93,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -61.581 NO -0.33261 -298.92 NO NO -2.8093 NO]"
]
},
"execution_count": 93,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(7)"
]
},
{
"cell_type": "code",
"execution_count": 94,
"metadata": {},
"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{79}{31} & 0 & 0 & -\\frac{101}{124} & 0 & 0 & \\frac{30381}{33604} & -\\frac{50401}{2294} & \\frac{1349}{62} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{323}{31} & 1 & 0 & -\\frac{369}{124} & 0 & 0 & \\frac{106969}{33604} & -\\frac{193965}{2294} & \\frac{7961}{62} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{1}{31} & 0 & 0 & -\\frac{9}{31} & 0 & 0 & \\frac{3761}{16802} & -\\frac{6491}{1147} & \\frac{904}{31} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{87}{31} & 0 & 0 & \\frac{123}{124} & 0 & 0 & -\\frac{33507}{33604} & \\frac{54983}{2294} & \\frac{157}{62} \\\\\n",
"{\\color{green}{s}}_{4} & 0 & 0 & 0 & -\\frac{19}{1240} & 0 & 0 & \\frac{531}{9920} & 1 & 0 & -\\frac{462257}{2688320} & \\frac{18101}{4960} & -\\frac{163}{4960} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 1 & -\\frac{438}{11111} & \\frac{4}{41} & 0 \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & \\frac{6}{31} & 0 & 1 & \\frac{1}{124} & 0 & 0 & \\frac{19929}{33604} & -\\frac{52689}{2294} & \\frac{91}{62} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{16}{31} & 0 & 0 & \\frac{13}{124} & 0 & 0 & \\frac{2149}{33604} & \\frac{763}{2294} & \\frac{12209}{62}\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 -79/31 0 0 -101/124 0 0 30381/33604 -50401/2294| 1349/62]\n",
"[ s1| 0 0 0 -323/31 1 0 -369/124 0 0 106969/33604 -193965/2294| 7961/62]\n",
"[ x3| 0 0 1 1/31 0 0 -9/31 0 0 3761/16802 -6491/1147| 904/31]\n",
"[ x2| 0 1 0 87/31 0 0 123/124 0 0 -33507/33604 54983/2294| 157/62]\n",
"[ s4| 0 0 0 -19/1240 0 0 531/9920 1 0 -462257/2688320 18101/4960| -163/4960]\n",
"[ s5| 0 0 0 0 0 0 0 0 1 -438/11111 4/41| 0]\n",
"[ s2| 0 0 0 6/31 0 1 1/124 0 0 19929/33604 -52689/2294| 91/62]\n",
"[---------------+-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+---------------]\n",
"[ Z| 0 0 0 16/31 0 0 13/124 0 0 2149/33604 763/2294| 12209/62]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(7,6)"
]
},
{
"cell_type": "code",
"execution_count": 95,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[ NO NO NO -33.684 NO NO NO NO NO -0.37191 NO]"
]
},
"execution_count": 95,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cocientesD(5)"
]
},
{
"cell_type": "code",
"execution_count": 96,
"metadata": {},
"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{1215251}{462257} & 0 & 0 & -\\frac{492833}{924514} & \\frac{2430480}{462257} & 0 & 0 & -\\frac{47595730}{17103509} & \\frac{9977945}{462257} \\\\\n",
"{\\color{green}{s}}_{1} & 0 & 0 & 0 & -\\frac{4947543}{462257} & 1 & 0 & -\\frac{1835037}{924514} & \\frac{8557520}{462257} & 0 & 0 & -\\frac{290654546}{17103509} & \\frac{59074065}{462257} \\\\\n",
"{x}_{3} & 0 & 0 & 1 & \\frac{5691}{462257} & 0 & 0 & -\\frac{203985}{924514} & \\frac{601760}{462257} & 0 & 0 & -\\frac{15536430}{17103509} & \\frac{13460235}{462257} \\\\\n",
"{x}_{2} & 0 & 1 & 0 & \\frac{1338375}{462257} & 0 & 0 & \\frac{630087}{924514} & -\\frac{2680560}{462257} & 0 & 0 & \\frac{47990306}{17103509} & \\frac{1258645}{462257} \\\\\n",
"{\\color{green}{s}}_{6} & 0 & 0 & 0 & \\frac{41192}{462257} & 0 & 0 & -\\frac{143901}{462257} & -\\frac{2688320}{462257} & 0 & 1 & -\\frac{9810742}{462257} & \\frac{88346}{462257} \\\\\n",
"{\\color{green}{s}}_{5} & 0 & 0 & 0 & \\frac{66576}{18952537} & 0 & 0 & -\\frac{232578}{18952537} & -\\frac{4344960}{18952537} & 1 & 0 & -\\frac{14007448}{18952537} & \\frac{142788}{18952537} \\\\\n",
"{\\color{green}{s}}_{2} & 0 & 0 & 0 & \\frac{65040}{462257} & 0 & 1 & \\frac{89069}{462257} & \\frac{1594320}{462257} & 0 & 0 & -\\frac{177559200}{17103509} & \\frac{626080}{462257} \\\\\n",
"\\hline\n",
" {Z_{ind}} & 0 & 0 & 0 & \\frac{235950}{462257} & 0 & 0 & \\frac{57665}{462257} & \\frac{171920}{462257} & 0 & 0 & \\frac{28902692}{17103509} & \\frac{91033000}{462257}\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 -1215251/462257 0 0 -492833/924514 2430480/462257 0 0 -47595730/17103509| 9977945/462257]\n",
"[ s1| 0 0 0 -4947543/462257 1 0 -1835037/924514 8557520/462257 0 0 -290654546/17103509| 59074065/462257]\n",
"[ x3| 0 0 1 5691/462257 0 0 -203985/924514 601760/462257 0 0 -15536430/17103509| 13460235/462257]\n",
"[ x2| 0 1 0 1338375/462257 0 0 630087/924514 -2680560/462257 0 0 47990306/17103509| 1258645/462257]\n",
"[ s6| 0 0 0 41192/462257 0 0 -143901/462257 -2688320/462257 0 1 -9810742/462257| 88346/462257]\n",
"[ s5| 0 0 0 66576/18952537 0 0 -232578/18952537 -4344960/18952537 1 0 -14007448/18952537| 142788/18952537]\n",
"[ s2| 0 0 0 65040/462257 0 1 89069/462257 1594320/462257 0 0 -177559200/17103509| 626080/462257]\n",
"[-------------------+---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+-------------------]\n",
"[ Z| 0 0 0 235950/462257 0 0 57665/462257 171920/462257 0 0 28902692/17103509| 91033000/462257]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cambio(5,10)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"v=matrix(QQ,1,nv,[0,-47/73,0,-(60/73)*(163/56),0,-202*163/(219*56),-(79/219)*(163/56)] ); R=-163/219;\n",
"añadeRes(v,R) # hecho con la segunda mejora"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"cocientesD(4)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"cambio(4,6)"
]
}
],
"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
}