Objetivo<-c(2,2.5,1.5,2) c1<-c(0.002,0.0035,0.0015,0.0011,-1,0,0) c2<-c(0,0,0,0,1,0,0) c3<-c(0,0,0,0,1,0,0) c4<-c(0.0005,0.00015,0.00065,0.001,0,-1,0) c5<-c(0,0,0,0,0,1,0) c6<-c(0,0,0,0,0,1,0) c7<-c(1/6500,1/5800,1/6200,1/5900,0,0,-1) c8<-c(0,0,0,0,0,0,6050) c9<-c(0,0,0,0,0,0,5950) c10<-c(1,1,1,1,0,0,0) Coef_restricciones<-rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10) Sol_restricciones<-c(0,0.001,0.003,0,0.00045,0.00055,0,1,1,1) Dir_restricciones<-c("=",">=","<=","=",">=","<=","=",">=","<=","=") library(lpSolve) lp(direction = "min",objective.in = Objetivo,const.mat = Coef_restricciones,const.dir = Dir_restricciones,const.rhs = Sol_restricciones) Solucion<-lp(direction = "min",objective.in = Objetivo,const.mat = Coef_restricciones,const.dir = Dir_restricciones,const.rhs = Sol_restricciones)