{ "cells": [ { "cell_type": "code", "execution_count": 4, "metadata": {}, "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", " \\max \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 17 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 12 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 7 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 20 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 6 x_{2} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 7 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 20 \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 6 x_{1} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 5 x_{3} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 16 \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 4 x_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 x_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 2 x_{3} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 2 x_{4} \\mspace{-6mu}&\\mspace{-6mu} \\leq \\mspace{-6mu}&\\mspace{-6mu} 14 \\\\\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{8821}{40}$$" ], "text/plain": [ "8821/40" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "" ], "text/latex": [ "$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(0,\\,\\frac{69}{20},\\,\\frac{161}{24},\\,\\frac{793}{120}\\right)$$" ], "text/plain": [ "(0, 69/20, 161/24, 793/120)" ] }, "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", " \\min \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 20 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 16 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 14 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 6 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 17 \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 6 y_{1} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 4 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 12 \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 7 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 5 y_{2} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} 2 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 7 \\\\\n", " \\mspace{-6mu}&\\mspace{-6mu} \\mspace{-6mu}&\\mspace{-6mu} 7 y_{1} \\mspace{-6mu}&\\mspace{-6mu} + \\mspace{-6mu}&\\mspace{-6mu} y_{2} \\mspace{-6mu}&\\mspace{-6mu} - \\mspace{-6mu}&\\mspace{-6mu} 2 y_{3} \\mspace{-6mu}&\\mspace{-6mu} \\geq \\mspace{-6mu}&\\mspace{-6mu} 20 \\\\\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{8821}{40}$$" ], "text/plain": [ "8821/40" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "" ], "text/latex": [ "$$\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\frac{149}{40},\\,\\frac{9}{2},\\,\\frac{423}{80}\\right)$$" ], "text/plain": [ "(149/40, 9/2, 423/80)" ] }, "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}}\\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", " {\\color{green}{s}}_{1} & 1 & 6 & -7 & 7 & 1 & 0 & 0 & 20 \\\\\n", "{\\color{green}{s}}_{2} & 6 & -7 & 5 & 1 & 0 & 1 & 0 & 16 \\\\\n", "{\\color{green}{s}}_{3} & 4 & 4 & 2 & -2 & 0 & 0 & 1 & 14 \\\\\n", "\\hline\n", " {Z_{ind}} & 17 & 12 & 7 & 20 & 0 & 0 & 0 & 0\n", "\\end{array}\\right)$$" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3|rhs]\n", "[---+---------------------------+---]\n", "[ s1| 1 6 -7 7 1 0 0| 20]\n", "[ s2| 6 -7 5 1 0 1 0| 16]\n", "[ s3| 4 4 2 -2 0 0 1| 14]\n", "[---+---------------------------+---]\n", "[ Z| 17 12 7 20 0 0 0| 0]" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ " \n", "(Es un problema de tipo maximizar)\n", " \n", "cambio( 1 , 4 )\n", "cambio( 2 , 3 )\n", "cambio( 3 , 2 )\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}_{4} & \\frac{35}{24} & 0 & 0 & 1 & \\frac{17}{120} & \\frac{1}{6} & \\frac{19}{240} & \\frac{793}{120} \\\\\n", "{x}_{3} & \\frac{47}{24} & 0 & 1 & 0 & \\frac{1}{24} & \\frac{1}{6} & \\frac{11}{48} & \\frac{161}{24} \\\\\n", "{x}_{2} & \\frac{3}{4} & 1 & 0 & 0 & \\frac{1}{20} & 0 & \\frac{7}{40} & \\frac{69}{20} \\\\\n", "\\hline\n", " {Z_{ind}} & -\\frac{279}{8} & 0 & 0 & 0 & -\\frac{149}{40} & -\\frac{9}{2} & -\\frac{423}{80} & \\frac{8821}{40}\n", "\\end{array}\\right)$$" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3| rhs]\n", "[-------+-------------------------------------------------------+-------]\n", "[ x4| 35/24 0 0 1 17/120 1/6 19/240|793/120]\n", "[ x3| 47/24 0 1 0 1/24 1/6 11/48| 161/24]\n", "[ x2| 3/4 1 0 0 1/20 0 7/40| 69/20]\n", "[-------+-------------------------------------------------------+-------]\n", "[ Z| -279/8 0 0 0 -149/40 -9/2 -423/80|8821/40]" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "El problema es óptimo\n" ] } ], "source": [ "tipo='max' # 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,[17,12,7,20] ) # Vector de costos\n", "\n", "A=matrix(QQ,m,n,[ # coeficientes de las variables x's\n", "[1,6,-7,7],\n", "[6,-7,5,1],\n", "[4,4,2,-2]\n", "])\n", "\n", "b=matrix(QQ,m,1,[20,16,14]) # En miles de unidades. 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", "############ HASTA AQUÍ LOS DATOS ################################\n", "####################################################################\n", "\n", "load('https://sage.unex.es/501708/SAGE/simplex2023Auto_MC.sage')\n", "#load('https://sage.unex.es/501708/SAGE/simplex2023AutoREVISED.sage')\n", "#load('https://sage.unex.es/501708/SAGE/simplex2023Auto.sage')\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1171.42857142857" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "8200/7." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "########################################\n", "\n", "Añadimos $x_2 \\geq 1172$" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8| rhs]\n", "[----------+-----------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ s1| 0 0 0 0 1 -6/7 0 -20/7 0 19/7 -1/7 0| 29900/7]\n", "[ s5| 0 0 0 0 0 1/7 0 -6/7 1 5/7 -8/7 0| 6800/7]\n", "[ s3| 0 0 0 0 0 -4/7 1 -4/7 0 -6/7 -24/7 0| 23200/7]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 1/7 0 -6/7 0 5/7 -8/7 0| 8200/7]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0| 400]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 1 0| 500]\n", "[ s8| 0 0 0 0 0 1/7 0 -6/7 0 5/7 -8/7 1| -4/7]\n", "[----------+-----------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ Z| 0 0 0 0 0 -1200/7 0 -5050/7 0 -1100/7 -4400/7 0|22475000/7]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "v=matrix(QQ,1,nv,[0,-1,0,0,0,0,0,0,0,0,0]) ; R=-1172\n", "añadeRes(v,R)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[ NO NO NO NO NO NO NO 841.67 NO NO 550.00 NO]" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cocientesD(8)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8| rhs]\n", "[-------+-----------------------------------------------------------------------------------------------+-------]\n", "[ s1| 0 0 0 0 1 -7/8 0 -11/4 0 21/8 0 -1/8| 8543/2]\n", "[ s5| 0 0 0 0 0 0 0 0 1 0 0 -1| 972]\n", "[ s3| 0 0 0 0 0 -1 1 2 0 -3 0 -3| 3316]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 0 0 0 0 0 0 -1| 1172]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0| 400]\n", "[ x4| 0 0 0 1 0 1/8 0 -3/4 0 5/8 0 7/8| 999/2]\n", "[ s7| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 1 -7/8| 1/2]\n", "[-------+-----------------------------------------------------------------------------------------------+-------]\n", "[ Z| 0 0 0 0 0 -250 0 -250 0 -550 0 -550|3210400]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "cambio(8,11)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "########################################\n", "\n", "Añadimos $x_2 \\leq 499$" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9| rhs]\n", "[-------+-------------------------------------------------------------------------------------------------------+-------]\n", "[ s1| 0 0 0 0 1 -7/8 0 -11/4 0 21/8 0 -1/8 0| 8543/2]\n", "[ s5| 0 0 0 0 0 0 0 0 1 0 0 -1 0| 972]\n", "[ s3| 0 0 0 0 0 -1 1 2 0 -3 0 -3 0| 3316]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 0 0 0 0 0 0 -1 0| 1172]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 1/8 0 -3/4 0 5/8 0 7/8 0| 999/2]\n", "[ s7| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 1 -7/8 0| 1/2]\n", "[ s9| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 0 -7/8 1| -1/2]\n", "[-------+-------------------------------------------------------------------------------------------------------+-------]\n", "[ Z| 0 0 0 0 0 -250 0 -250 0 -550 0 -550 0|3210400]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "v=matrix(QQ,1,nv,[0,0,0,1,0,0,0,0,0,0,0,0]) ; R=499\n", "añadeRes(v,R)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[ NO NO NO NO NO 2000.0 NO NO NO 880.00 NO 628.57 NO]" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cocientesD(9)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9| rhs]\n", "[----------+----------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ s1| 0 0 0 0 1 -6/7 0 -20/7 0 19/7 0 0 -1/7| 29901/7]\n", "[ s5| 0 0 0 0 0 1/7 0 -6/7 1 5/7 0 0 -8/7| 6808/7]\n", "[ s3| 0 0 0 0 0 -4/7 1 -4/7 0 -6/7 0 0 -24/7| 23224/7]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 1/7 0 -6/7 0 5/7 0 0 -8/7| 8208/7]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 0 0 1| 499]\n", "[ s7| 0 0 0 0 0 0 0 0 0 0 1 0 -1| 1]\n", "[ s8| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 1 -8/7| 4/7]\n", "[----------+----------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ Z| 0 0 0 0 0 -1200/7 0 -5050/7 0 -1100/7 0 0 -4400/7|22470600/7]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "cambio(9,12)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "########################################\n", "\n", "Añadimos $x_2 \\geq 1173$" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10| rhs]\n", "[----------+---------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ s1| 0 0 0 0 1 -6/7 0 -20/7 0 19/7 0 0 -1/7 0| 29901/7]\n", "[ s5| 0 0 0 0 0 1/7 0 -6/7 1 5/7 0 0 -8/7 0| 6808/7]\n", "[ s3| 0 0 0 0 0 -4/7 1 -4/7 0 -6/7 0 0 -24/7 0| 23224/7]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 1/7 0 -6/7 0 5/7 0 0 -8/7 0| 8208/7]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 0 0 1 0| 499]\n", "[ s7| 0 0 0 0 0 0 0 0 0 0 1 0 -1 0| 1]\n", "[ s8| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 1 -8/7 0| 4/7]\n", "[ s10| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 0 -8/7 1| -3/7]\n", "[----------+---------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ Z| 0 0 0 0 0 -1200/7 0 -5050/7 0 -1100/7 0 0 -4400/7 0|22470600/7]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "v=matrix(QQ,1,nv,[0,-1,0,0,0,0,0,0,0,0,0,0,0]) ; R=-1173\n", "añadeRes(v,R)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10| rhs]\n", "[-------+---------------------------------------------------------------------------------------------------------------+-------]\n", "[ s1| 0 0 0 0 1 -7/8 0 -11/4 0 21/8 0 0 0 -1/8|34173/8]\n", "[ s5| 0 0 0 0 0 0 0 0 1 0 0 0 0 -1| 973]\n", "[ s3| 0 0 0 0 0 -1 1 2 0 -3 0 0 0 -3| 3319]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 0 0 0 0 0 0 0 0 -1| 1173]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 1/8 0 -3/4 0 5/8 0 0 0 7/8| 3989/8]\n", "[ s7| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 1 0 0 -7/8| 11/8]\n", "[ s8| 0 0 0 0 0 0 0 0 0 0 0 1 0 -1| 1]\n", "[ s9| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 0 0 1 -7/8| 3/8]\n", "[-------+---------------------------------------------------------------------------------------------------------------+-------]\n", "[ Z| 0 0 0 0 0 -250 0 -250 0 -550 0 0 0 -550|3209850]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "cambio(10,13)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "########################################\n", "\n", "Añadimos $x_4 \\leq 498$ (P7)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11| rhs]\n", "[-------+-----------------------------------------------------------------------------------------------------------------------+-------]\n", "[ s1| 0 0 0 0 1 -7/8 0 -11/4 0 21/8 0 0 0 -1/8 0|34173/8]\n", "[ s5| 0 0 0 0 0 0 0 0 1 0 0 0 0 -1 0| 973]\n", "[ s3| 0 0 0 0 0 -1 1 2 0 -3 0 0 0 -3 0| 3319]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0| 1173]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 1/8 0 -3/4 0 5/8 0 0 0 7/8 0| 3989/8]\n", "[ s7| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 1 0 0 -7/8 0| 11/8]\n", "[ s8| 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 0| 1]\n", "[ s9| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 0 0 1 -7/8 0| 3/8]\n", "[ s11| 0 0 0 0 0 -1/8 0 3/4 0 -5/8 0 0 0 -7/8 1| -5/8]\n", "[-------+-----------------------------------------------------------------------------------------------------------------------+-------]\n", "[ Z| 0 0 0 0 0 -250 0 -250 0 -550 0 0 0 -550 0|3209850]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "v=matrix(QQ,1,nv,[0,0,0,1,0,0,0,0,0,0,0,0,0,0]) ; R=498\n", "añadeRes(v,R)" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11| rhs]\n", "[----------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ s1| 0 0 0 0 1 -6/7 0 -20/7 0 19/7 0 0 0 0 -1/7| 29902/7]\n", "[ s5| 0 0 0 0 0 1/7 0 -6/7 1 5/7 0 0 0 0 -8/7| 6816/7]\n", "[ s3| 0 0 0 0 0 -4/7 1 -4/7 0 -6/7 0 0 0 0 -24/7| 23248/7]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 1/7 0 -6/7 0 5/7 0 0 0 0 -8/7| 8216/7]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1| 498]\n", "[ s7| 0 0 0 0 0 0 0 0 0 0 1 0 0 0 -1| 2]\n", "[ s8| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 1 0 0 -8/7| 12/7]\n", "[ s9| 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1| 1]\n", "[ s10| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 0 0 1 -8/7| 5/7]\n", "[----------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ Z| 0 0 0 0 0 -1200/7 0 -5050/7 0 -1100/7 0 0 0 0 -4400/7|22466200/7]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "cambio(11,14)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "########################################\n", "\n", "Añadimos $x_2 \\leq 1173$ (P9)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12| rhs]\n", "[----------+-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ s1| 0 0 0 0 1 -6/7 0 -20/7 0 19/7 0 0 0 0 -1/7 0| 29902/7]\n", "[ s5| 0 0 0 0 0 1/7 0 -6/7 1 5/7 0 0 0 0 -8/7 0| 6816/7]\n", "[ s3| 0 0 0 0 0 -4/7 1 -4/7 0 -6/7 0 0 0 0 -24/7 0| 23248/7]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 1/7 0 -6/7 0 5/7 0 0 0 0 -8/7 0| 8216/7]\n", "[ x3| 0 0 1 0 0 0 0 0 0 -1 0 0 0 0 0 0| 400]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0| 498]\n", "[ s7| 0 0 0 0 0 0 0 0 0 0 1 0 0 0 -1 0| 2]\n", "[ s8| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 1 0 0 -8/7 0| 12/7]\n", "[ s9| 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 0| 1]\n", "[ s10| 0 0 0 0 0 1/7 0 -6/7 0 5/7 0 0 0 1 -8/7 0| 5/7]\n", "[ s12| 0 0 0 0 0 -1/7 0 6/7 0 -5/7 0 0 0 0 8/7 1| -5/7]\n", "[----------+-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+----------]\n", "[ Z| 0 0 0 0 0 -1200/7 0 -5050/7 0 -1100/7 0 0 0 0 -4400/7 0|22466200/7]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "v=matrix(QQ,1,nv,[0,1,0,0,0,0,0,0,0,0,0,0,0,0,0]) ; R=1173\n", "añadeRes(v,R)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "[ B| x1 x2 x3 x4 s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12| rhs]\n", "[-------+-------------------------------------------------------------------------------------------------------------------------------+-------]\n", "[ s1| 0 0 0 0 1 -7/5 0 2/5 0 0 0 0 0 0 21/5 19/5| 4269]\n", "[ s5| 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1| 973]\n", "[ s3| 0 0 0 0 0 -2/5 1 -8/5 0 0 0 0 0 0 -24/5 -6/5| 3322]\n", "[ x1| 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0| 300]\n", "[ x2| 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1| 1173]\n", "[ x3| 0 0 1 0 0 1/5 0 -6/5 0 0 0 0 0 0 -8/5 -7/5| 401]\n", "[ x4| 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0| 498]\n", "[ s7| 0 0 0 0 0 0 0 0 0 0 1 0 0 0 -1 0| 2]\n", "[ s8| 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1| 1]\n", "[ s9| 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 0| 1]\n", "[ s10| 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1| 0]\n", "[ s6| 0 0 0 0 0 1/5 0 -6/5 0 1 0 0 0 0 -8/5 -7/5| 1]\n", "[-------+-------------------------------------------------------------------------------------------------------------------------------+-------]\n", "[ Z| 0 0 0 0 0 -140 0 -910 0 0 0 0 0 0 -880 -220|3209300]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "cambio(12,10)" ] } ], "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 }