library(lpSolve)

#Defino el modelo

#Funcional
coef_funcional<-c(2,2.5,1.5,2)

#Restricciones
matriz_restr<-matrix(c(6500,-6500,0.002,-0.002,0.0005,-0.0005,5800,-5800,0.0035,-0.0035,0.00015,-0.00015,6200,-6200,0.0015,-0.0015,0.00065,-0.00065,5900,-5900,0.0011,-0.0011,0.001,-0.001),ncol=4)
restr<-c(5950,-6050,0.001,-0.003,0.00045,-0.00055)

#Minimización 
desigualdad<-rep('>=',6)

#Solucion: Costo mínimo 

costo_min<-lp('min',coef_funcional,matriz_restr,desigualdad,restr)
costo_min

#Kg de A,B,C,D para el costo mín

costo_min$solution