HesauHugo_Trabalho2Excercicio1Jupyter.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "7acf1c67",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pyomo.environ as pyEnv"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "e23fd4a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "#Produção de fumo bruto em toneladas da região r no mês t\n",
    "p_rt =[[230000,150000,0,0,0,0,140000,180000,260000,280000,380000,350000]\n",
    "       ,[120000,0,0,0,0,110000,250000,360000,380000,370000,350000,230000]\n",
    "       ,[210000,110000,0,0,0,0,130000,240000,370000,400000,420000,380000]\n",
    "       ,[0,0,0,0,100000,110000,150000,250000,260000,250000,200000,120000]\n",
    "      ]\n",
    "#Capacidade de armazenamento de fumo bruto na fábrica fem toneladas\n",
    "eb_f =[1500000,750000,1250000,1000000]\n",
    "#Custo do estoque de fumo bruto por tonelada na fábrica f\n",
    "cb_f =[8.1,7.2,7.6,5.4]\n",
    "#Capacidade de processamento de fumo na fábrica f em toneladas\n",
    "f_f = [240000,180000,160000,170000]\n",
    "#Custo do processamento de fumo por tonelada na fábrica f\n",
    "cf_f = [10.3,12.4,9.7,10.5]\n",
    "#Capacidade de estoque de fumo processado\n",
    "ep_f = [1000000,800000,1000000,900000]\n",
    "#Custo do estoque de fumo processado por tonelada na fábrica f\n",
    "cp_f = [5.3,4.1,4.6,3.8]\n",
    "#Custo do transporte de fumo por tonelada da região r  para fábrica f \n",
    "#(RS,SC,PR,MS) (blumenau,rio negro,cruz alta,uberaba)\n",
    "ct_rf = [[10.8,12,7.2,26]\n",
    "         ,[3.6,4,6.8,25.2]\n",
    "         ,[10.8,9.2,14.4,21.6]\n",
    "         ,[16.4,14.8,23.4,10.8]]\n",
    "\n",
    "#Custo do transporte de fumo por tonelada da fábrica f  para fábrica f'\n",
    "#(blumenau,rio negro,cruz alta,uberaba)(blumenau,rio negro,cruz alta,uberaba)\n",
    "ct_ffLinha = [[0,7.2,8.4,18.4]\n",
    "             ,[7.2,0,10,17.20]\n",
    "             ,[8.4,10,0,22.8]\n",
    "              ,[18.4,17.2,22.8,0]\n",
    "             ]\n",
    "#Demanda total de fumo no mês t em toneladas\n",
    "d_t = [300000,500000,700000,900000,800000,800000,700000,600000,600000,500000,500000,400000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "f969da1f",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "GLPSOL: GLPK LP/MIP Solver, v4.65\n",
      "Parameter(s) specified in the command line:\n",
      " --write C:\\Users\\User\\AppData\\Local\\Temp\\tmpt0jv60aq.glpk.raw --wglp C:\\Users\\User\\AppData\\Local\\Temp\\tmpywfr5w1a.glpk.glp\n",
      " --cpxlp C:\\Users\\User\\AppData\\Local\\Temp\\tmpz12ydvr0.pyomo.lp\n",
      "Reading problem data from 'C:\\Users\\User\\AppData\\Local\\Temp\\tmpz12ydvr0.pyomo.lp'...\n",
      "301 rows, 673 columns, 1473 non-zeros\n",
      "3680 lines were read\n",
      "Writing problem data to 'C:\\Users\\User\\AppData\\Local\\Temp\\tmpywfr5w1a.glpk.glp'...\n",
      "3279 lines were written\n",
      "GLPK Simplex Optimizer, v4.65\n",
      "301 rows, 673 columns, 1473 non-zeros\n",
      "Preprocessing...\n",
      "140 rows, 592 columns, 1184 non-zeros\n",
      "Scaling...\n",
      " A: min|aij| =  1.000e+00  max|aij| =  1.000e+00  ratio =  1.000e+00\n",
      "Problem data seem to be well scaled\n",
      "Constructing initial basis...\n",
      "Size of triangular part is 131\n",
      "      0: obj =   1.900000000e+06 inf =   6.900e+06 (11)\n",
      "     53: obj =   1.465490000e+08 inf =   0.000e+00 (0)\n",
      "*   124: obj =   1.072970000e+08 inf =   0.000e+00 (0)\n",
      "OPTIMAL LP SOLUTION FOUND\n",
      "Time used:   0.0 secs\n",
      "Memory used: 0.6 Mb (603888 bytes)\n",
      "Writing basic solution to 'C:\\Users\\User\\AppData\\Local\\Temp\\tmpt0jv60aq.glpk.raw'...\n",
      "983 lines were written\n"
     ]
    }
   ],
   "source": [
    "#indices\n",
    "modelo = pyEnv.ConcreteModel()\n",
    "\n",
    "modelo.R = pyEnv.RangeSet(len([1,2,3,4]))\n",
    "modelo.F = pyEnv.RangeSet(len([1,2,3,4]))\n",
    "modelo.FLinha = pyEnv.RangeSet(len([1,2,3,4]))\n",
    "modelo.T = pyEnv.RangeSet(len([1,2,3,4,5,6,7,8,9,10,11,12]))\n",
    "\n",
    "#Variáveis\n",
    "modelo.x_rft = pyEnv.Var(modelo.R,modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.eb_ft = pyEnv.Var(modelo.F,modelo.T, within = pyEnv.NonNegativeReals)\n",
    "modelo.y_ft = pyEnv.Var(modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.ep_ft = pyEnv.Var(modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.s_ft = pyEnv.Var(modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.h_fLinhaft = pyEnv.Var(modelo.FLinha,modelo.F, modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.z_fLinhaft = pyEnv.Var(modelo.FLinha,modelo.F, modelo.T , within = pyEnv.NonNegativeReals)\n",
    "\n",
    "\n",
    "#parâmetros\n",
    "modelo.P_rt = pyEnv.Param(modelo.R,modelo.T,initialize = lambda modelo,r,t:p_rt[r-1][t-1])\n",
    "modelo.EB_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:eb_f[f-1])\n",
    "modelo.CB_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:cb_f[f-1])\n",
    "modelo.F_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:f_f[f-1])\n",
    "modelo.CF_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:cf_f[f-1])\n",
    "modelo.EP_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:ep_f[f-1])\n",
    "modelo.CP_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:cp_f[f-1])\n",
    "modelo.CT_rf = pyEnv.Param(modelo.R,modelo.F,initialize = lambda modelo,r,f:ct_rf[r-1][f-1])\n",
    "modelo.CT_FFLinha =  pyEnv.Param(modelo.FLinha,modelo.F,initialize = lambda modelo,f,fLinha:ct_ffLinha[f-1][fLinha-1])\n",
    "modelo.D_t =  pyEnv.Param(modelo.T,initialize = lambda modelo,t:d_t[t-1])\n",
    "\n",
    "#Função objetivo\n",
    "def func_objetivo(modelo):\n",
    "    return sum(\n",
    "        modelo.x_rft[r,f,t] * modelo.CT_rf[r,f] for r in modelo.R for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.h_fLinhaft[f,fLinha,t] * modelo.CT_FFLinha[f,fLinha] for f in modelo.F for fLinha in modelo.FLinha for t in modelo.T) + sum(\n",
    "        modelo.z_fLinhaft[f,fLinha,t] * modelo.CT_FFLinha[f,fLinha] for f in modelo.F for fLinha in modelo.FLinha for t in modelo.T)  + sum(\n",
    "        modelo.y_ft[f,t] * modelo.CF_f[f] for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.eb_ft[f,t] * modelo.CB_f[f] for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.ep_ft[f,t] * modelo.CP_f[f] for f in modelo.F for t in modelo.T) \n",
    "\n",
    "modelo.objetivo = pyEnv.Objective(rule = func_objetivo ,sense = pyEnv.minimize) \n",
    "\n",
    "#Capacidade de distribuição de fumo bruto\n",
    "\n",
    "def CapacidadeDistribuicaoDeFumoBruto(modelo,r,t):\n",
    "    return modelo.P_rt[r,t] >= sum(modelo.x_rft[r,f,t] for f in modelo.F)\n",
    "\n",
    "modelo.capacidadeDistribuicaoDeFumoBruto = pyEnv.Constraint(modelo.R,modelo.T,rule = CapacidadeDistribuicaoDeFumoBruto)\n",
    "\n",
    "#Estoque Bruto e disponível para transporte\n",
    "def EstoqueBrutoDisponivelParaTransporte(modelo,f,t): \n",
    "    if t>1 :\n",
    "        return   modelo.eb_ft[f,t]  == modelo.eb_ft[f,t-1] + sum(\n",
    "            modelo.x_rft[r,f,t] for r in modelo.R) + sum(\n",
    "            modelo.h_fLinhaft[fLinha,f,t] for fLinha in modelo.FLinha if f != fLinha )  - sum(\n",
    "            modelo.h_fLinhaft[f,fLinha,t] for fLinha in modelo.FLinha if f != fLinha) - modelo.y_ft[f,t]\n",
    "    \n",
    "    else:\n",
    "        return   modelo.eb_ft[f,t]  == modelo.eb_ft[f,t] + sum(\n",
    "            modelo.x_rft[r,f,t] for r in modelo.R) + sum(\n",
    "            modelo.h_fLinhaft[fLinha,f,t] for fLinha in modelo.FLinha if f != fLinha )  - sum(\n",
    "            modelo.h_fLinhaft[f,fLinha,t] for fLinha in modelo.FLinha if f != fLinha) - modelo.y_ft[f,t]\n",
    "    \n",
    "modelo.estoqueBrutoDisponivelParaTransporte = pyEnv.Constraint(modelo.F,modelo.T, rule = EstoqueBrutoDisponivelParaTransporte)    \n",
    "\n",
    "#Demanda\n",
    "def Demanda(modelo,t):\n",
    "    return modelo.D_t[t]  == sum(modelo.s_ft[f,t] for f in modelo.F)\n",
    "\n",
    "modelo.demanda = pyEnv.Constraint(modelo.T, rule = Demanda)\n",
    "\n",
    "#Estoque processado e disponível para transporte\n",
    "def EstoqueProcessadoDisponivelTransporte(modelo,f,t):\n",
    "    if t > 1 :\n",
    "        return modelo.ep_ft[f,t] == modelo.ep_ft[f,t-1] +  sum(\n",
    "            modelo.z_fLinhaft[fLinha,f,t] for fLinha in modelo.FLinha if f != fLinha) - sum(\n",
    "            modelo.z_fLinhaft[f,fLinha,t] for fLinha in modelo.FLinha if f != fLinha) + modelo.y_ft[f,t] -  modelo.s_ft[f,t]\n",
    "    else:\n",
    "        return modelo.ep_ft[f,t] == modelo.ep_ft[f,t] +  sum(\n",
    "            modelo.z_fLinhaft[fLinha,f,t] for fLinha in modelo.FLinha if f != fLinha) - sum(\n",
    "            modelo.z_fLinhaft[f,fLinha,t] for fLinha in modelo.FLinha if f != fLinha) + modelo.y_ft[f,t] -  modelo.s_ft[f,t]\n",
    "\n",
    "modelo.estoqueProcessadoDisponivelTransporte = pyEnv.Constraint(modelo.F,modelo.T, rule = EstoqueProcessadoDisponivelTransporte)    \n",
    "\n",
    "#Capacidade Estoque bruto\n",
    "def CapacidadeEstoqueBruto(modelo,f,t):\n",
    "    return modelo.eb_ft[f,t] <= modelo.EB_f[f]\n",
    "\n",
    "modelo.capacidadeEstoqueBruto = pyEnv.Constraint(modelo.F, modelo.T,rule = CapacidadeEstoqueBruto)\n",
    "\n",
    "#Capacidade de processamento\n",
    "def CapacidadeProcessamento(modelo,f,t):\n",
    "    return modelo.y_ft[f,t] <= modelo.F_f[f]\n",
    "\n",
    "modelo.capacidadeProcessamento = pyEnv.Constraint(modelo.F, modelo.T,rule = CapacidadeProcessamento)\n",
    "\n",
    "#Capacidade do Estoque processamento\n",
    "def CapacidadeEstoqueProcessado(modelo,f,t):\n",
    "    return modelo.ep_ft[f,t] <= modelo.EP_f[f]\n",
    "\n",
    "modelo.capacidadeEstoqueProcessado = pyEnv.Constraint(modelo.F, modelo.T,rule = CapacidadeEstoqueProcessado)\n",
    "\n",
    "solver = pyEnv.SolverFactory('glpk')\n",
    "resultado_objetivo = solver.solve(modelo,tee = True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "461af9b5",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Problem: \n",
      "- Name: unknown\n",
      "  Lower bound: 107297000.0\n",
      "  Upper bound: 107297000.0\n",
      "  Number of objectives: 1\n",
      "  Number of constraints: 301\n",
      "  Number of variables: 673\n",
      "  Number of nonzeros: 1473\n",
      "  Sense: minimize\n",
      "Solver: \n",
      "- Status: ok\n",
      "  Termination condition: optimal\n",
      "  Statistics: \n",
      "    Branch and bound: \n",
      "      Number of bounded subproblems: 0\n",
      "      Number of created subproblems: 0\n",
      "  Error rc: 0\n",
      "  Time: 0.050728797912597656\n",
      "Solution: \n",
      "- number of solutions: 0\n",
      "  number of solutions displayed: 0\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(resultado_objetivo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b476ee63",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mes 0\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 160000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 120000.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 20000.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 1\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 0.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 0.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 2\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 0.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 0.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 3\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 0.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 0.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 4\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 0.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 0.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 5\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 0.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 110000.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 110000.0\n",
      "\n",
      "mes 6\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 140000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 240000.0\n",
      "De 1 para 1 -----> 10000.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 130000.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 150000.0\n",
      "\n",
      "mes 7\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 160000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 240000.0\n",
      "De 1 para 1 -----> 120000.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 80000.0\n",
      "\n",
      "mes 8\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 160000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 240000.0\n",
      "De 1 para 1 -----> 140000.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 60000.0\n",
      "\n",
      "mes 9\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 130000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 240000.0\n",
      "De 1 para 1 -----> 130000.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 10\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 160000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 240000.0\n",
      "De 1 para 1 -----> 110000.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n",
      "mes 11\n",
      "De 0 para 0 -----> 0.0\n",
      "De 0 para 1 -----> 0.0\n",
      "De 0 para 2 -----> 160000.0\n",
      "De 0 para 3 -----> 0.0\n",
      "\n",
      "De 1 para 0 -----> 230000.0\n",
      "De 1 para 1 -----> 0.0\n",
      "De 1 para 2 -----> 0.0\n",
      "De 1 para 3 -----> 0.0\n",
      "\n",
      "De 2 para 0 -----> 0.0\n",
      "De 2 para 1 -----> 0.0\n",
      "De 2 para 2 -----> 0.0\n",
      "De 2 para 3 -----> 0.0\n",
      "\n",
      "De 3 para 0 -----> 0.0\n",
      "De 3 para 1 -----> 0.0\n",
      "De 3 para 2 -----> 0.0\n",
      "De 3 para 3 -----> 0.0\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for t in modelo.T:\n",
    "    print(\"mes\",t-1)\n",
    "    for r in modelo.R:\n",
    "        for f in modelo.F:\n",
    "            print(\"De\", r-1, \"para\", f-1, \"----->\", modelo.x_rft[r,f,t]())\n",
    "        print(\"\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2704fd86",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}