#Defino las variables del problema #Pij:Peso de la carga i (con i=1,2,3,4) en el compartimento j (con j=D,C,T)[tn] #Vij:Volumen de la carga i (con i=1,2,3,4) en el compartimento j (con j=D,C,T) [pies3/tn]. #Siendo Vij = Ci*Pij con Ci = 500, 700, 600, 400 #Defino las variables totalizadoras #-P1+P1D+P1C+P1T=0 #-P2+P2D+P2C+P2T=0 #-P3+P3D+P3C+P3T=0 #-P4+P4D+P4C+P4T=0 #-PD+P1D+P2D+P3D+P4D=0 #-PC+P1C+P2C+P3C+P4C=0 #-PT+P1T+P2T+P3T+P4T=0 #Defino el funcional o función objetivo que estoy buscando maximizar #[MAX] [MAX] Z= 320.P1+400.P2+360.P3+290.P4 z<-c(320,400,360,290,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) #Sujeto a: #Armo la matriz de coeficientes tecnológicos A A<-read_csv("~/Documents/FIUBA/IO/Tabla_2.9.csv") #Defino el vector de términos independientes b<-c(0,0,0,0,0,0,0,0,0,0,12,16,25,13,12,18,10,7000,9000,5000) dir<-c("=","=","=","=","=","=","=","=","=","=","=","=","<=","<=","<=","<=","<=","<=","<=","<=") library(linprog) library(lpSolve) Sol<-lp("max",z,A,dir,b) Detalle<-lp("max",z,A,dir,b)$solution print(Sol) print(Detalle)