Trabalho2Exercício1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "280d1700",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pyomo.environ as pyEnv"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "ffd04fff",
   "metadata": {},
   "outputs": [],
   "source": [
    "p_rt = [[260000, 280000, 380000, 350000, 230000, 150000, 0, 0, 0, 0, 140000, 180000], \n",
    "     [380000, 370000, 350000, 230000, 120000, 0, 0, 0, 0, 110000, 250000, 360000], \n",
    "     [370000, 400000, 420000, 380000, 210000, 110000, 0, 0, 0, 0, 130000, 240000], \n",
    "     [260000, 250000, 200000, 120000, 0, 0, 0, 0, 100000, 110000, 150000, 250000]]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "0b17fd8d",
   "metadata": {},
   "outputs": [],
   "source": [
    "ep_f = [240000,180000,160000,170000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "616d8d77",
   "metadata": {},
   "outputs": [],
   "source": [
    "cp_f = [10.3,12.4,9.7,10.5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "51e36034",
   "metadata": {},
   "outputs": [],
   "source": [
    "eb_f = [1500000,750000,1250000,1000000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "bc9e3657",
   "metadata": {},
   "outputs": [],
   "source": [
    "cb_f =[8.10,7.20,7.60,5.40]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "99907bca",
   "metadata": {},
   "outputs": [],
   "source": [
    "fp_f = [1000000,800000,1000000,900000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "1ab58881",
   "metadata": {},
   "outputs": [],
   "source": [
    "cfp_f = [5.3,4.1,4.6,3.8]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "5ef37199",
   "metadata": {},
   "outputs": [],
   "source": [
    "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]]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "ef8593af",
   "metadata": {},
   "outputs": [],
   "source": [
    "ct_ffLinha =[[0, 7.2, 8.4, 18.4], \n",
    "       [7.2, 0, 10, 17.2], \n",
    "       [8.4, 10, 0, 22.8], \n",
    "       [18.4, 17.2, 22.8, 0]]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "2eababa0",
   "metadata": {},
   "outputs": [],
   "source": [
    "d_t = [600000, 500000, 500000, 400000, 300000, 500000, 700000, 900000, 800000, 800000, 700000, 600000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "465cf79d",
   "metadata": {},
   "outputs": [],
   "source": [
    "#indices\n",
    "r=[1,2,3,4]\n",
    "f=[1,2,3,4]\n",
    "fLinha=[1,2,3,4]\n",
    "t=[1,2,3,4,5,6,7,8,9,10,11,12]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "310d64e1",
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "modelo = pyEnv.ConcreteModel()\n",
    "\n",
    "modelo.R = pyEnv.RangeSet(len(r))\n",
    "modelo.F = pyEnv.RangeSet(len(f))\n",
    "modelo.FLinha = pyEnv.RangeSet(len(fLinha))\n",
    "modelo.T = pyEnv.RangeSet(len(t))\n",
    "\n",
    "\n",
    "\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.z_fLinhaft = pyEnv.Var(modelo.FLinha,modelo.F, modelo.T , within = pyEnv.NonNegativeReals)\n",
    "modelo.h_fLinhaft = pyEnv.Var(modelo.FLinha,modelo.F, modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.s_ft = pyEnv.Var(modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.he = pyEnv.Var(modelo.FLinha,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "modelo.ze = pyEnv.Var(modelo.F,modelo.T,within = pyEnv.NonNegativeReals)\n",
    "\n",
    "modelo.P_rt = pyEnv.Param(modelo.R,modelo.T,initialize = lambda modelo,r,t:p_rt[r-1][t-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.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.Fp_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:fp_f[f-1])\n",
    "modelo.CFP_f = pyEnv.Param(modelo.F,initialize = lambda modelo,f:cfp_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_FLinhaF =  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(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.eb_ft[f,t] * modelo.CB_f[f] for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.h_fLinhaft[fLinha,f,t] * modelo.CT_FLinhaF[fLinha,f] for fLinha in modelo.FLinha for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.y_ft[f,t] * modelo.CFP_f[f] for f in modelo.F for t in modelo.T) + sum(\n",
    "        modelo.z_fLinhaft[fLinha,f,t] * modelo.CT_FLinhaF[fLinha,f] for fLinha in modelo.FLinha 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",
    "\n",
    "modelo.objetivo = pyEnv.Objective(rule = func_objetivo ,sense = pyEnv.minimize) \n",
    "\n",
    "#Capacidade de distribuição de fumo bruto\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[fLinha,f,t] for fLinha in modelo.FLinha if f != fLinha) - modelo.he[f,t] - modelo.y_ft[f,t]\n",
    "    else:\n",
    "         return  modelo.eb_ft[f,t] ==  modelo.eb_ft[f,t] + sum(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) - modelo.he[f,t] - modelo.y_ft[f,t]\n",
    "\n",
    "modelo.estoqueBrutoDisponivelParaTransporte = pyEnv.Constraint(modelo.F,modelo.T, rule = EstoqueBrutoDisponivelParaTransporte)\n",
    "\n",
    "def f_restr4(modelo, fLinha, f, t):\n",
    "    if fLinha == f:\n",
    "        return modelo.he[f,t] == sum(modelo.h_fLinhaft[fLinha,f,t] for f in modelo.F)\n",
    "    else:\n",
    "        return pyEnv.Constraint.Skip\n",
    "    \n",
    "modelo.restr_4 = pyEnv.Constraint(modelo.FLinha, modelo.F, modelo.T, rule=f_restr4)\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.y_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)\n",
    "    else:\n",
    "        return modelo.ep_ft[f,t] == modelo.ep_ft[f,t] + modelo.y_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)\n",
    "\n",
    "modelo.estoqueProcessadoDisponivelTransporte = pyEnv.Constraint(modelo.F, modelo.T, rule = EstoqueProcessadoDisponivelTransporte)\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",
    "\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.Fp_f[f]\n",
    "\n",
    "modelo.capacidadeProcessamento = pyEnv.Constraint(modelo.F, modelo.T,rule = CapacidadeProcessamento)\n",
    "\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"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "046cf7b0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "GLPSOL--GLPK LP/MIP Solver 5.0\n",
      "Parameter(s) specified in the command line:\n",
      " --write C:\\Users\\User\\AppData\\Local\\Temp\\tmpfz4ln9h5.glpk.raw --wglp C:\\Users\\User\\AppData\\Local\\Temp\\tmp9r261l7f.glpk.glp\n",
      " --cpxlp C:\\Users\\User\\AppData\\Local\\Temp\\tmpzmcrbwjk.pyomo.lp\n",
      "Reading problem data from 'C:\\Users\\User\\AppData\\Local\\Temp\\tmpzmcrbwjk.pyomo.lp'...\n",
      "349 rows, 769 columns, 1437 non-zeros\n",
      "3884 lines were read\n",
      "Writing problem data to 'C:\\Users\\User\\AppData\\Local\\Temp\\tmp9r261l7f.glpk.glp'...\n",
      "3387 lines were written\n",
      "GLPK Simplex Optimizer 5.0\n",
      "349 rows, 769 columns, 1437 non-zeros\n",
      "Preprocessing...\n",
      "188 rows, 628 columns, 1088 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 180\n",
      "*     0: obj =   0.000000000e+00 inf =   0.000e+00 (98)\n",
      "*    63: obj =   0.000000000e+00 inf =   0.000e+00 (0)\n",
      "OPTIMAL LP SOLUTION FOUND\n",
      "Time used:   0.0 secs\n",
      "Memory used: 0.6 Mb (652416 bytes)\n",
      "Writing basic solution to 'C:\\Users\\User\\AppData\\Local\\Temp\\tmpfz4ln9h5.glpk.raw'...\n",
      "1127 lines were written\n"
     ]
    }
   ],
   "source": [
    "solver = pyEnv.SolverFactory('glpk')\n",
    "resultado_objetivo = solver.solve(modelo,tee = True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "d75f70d8",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Problem: \n",
      "- Name: unknown\n",
      "  Lower bound: 0.0\n",
      "  Upper bound: 0.0\n",
      "  Number of objectives: 1\n",
      "  Number of constraints: 349\n",
      "  Number of variables: 769\n",
      "  Number of nonzeros: 1437\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.06034135818481445\n",
      "Solution: \n",
      "- number of solutions: 0\n",
      "  number of solutions displayed: 0\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(resultado_objetivo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "857fa4c9",
   "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 -----> 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 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 -----> 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 6\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 7\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 8\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 9\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 10\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 11\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"
     ]
    }
   ],
   "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": "7fda0fc6",
   "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
}