fomula (1).R

SolverLP<-function(model,method="CPLEX_LP",decimal=0)
{
  library(Rglpk)
  model1.lp<-Rglpk_read_file(model,type=method,
                             verbose=F)
  model1.lp.sol<-Rglpk_solve_LP(model1.lp$objective,
                                model1.lp$constraints[[1]],
                                model1.lp$constraints[[2]],
                                model1.lp$constraints[[3]],
                                model1.lp$bounds,
                                model1.lp$types,
                                model1.lp$maximum)
  library(xtable)
  model1.lp.sol.df<-as.data.frame(model1.lp.sol$solution)
  model1.lp.sol.df<-rbind(model1.lp.sol.df,c(model1.lp.sol$optimum))
  rownames(model1.lp.sol.df)<-c(attr(model1.lp,"objective_vars_names"),"obj")
  colnames(model1.lp.sol.df)<-"Solution"
  table.sol<-xtable(model1.lp.sol.df,digits=decimal)
  results<-list(sol=model1.lp.sol,df=model1.lp.sol.df)
  return(results)
  printRangesGLPK(model1.lp.sol, fname = "report.txt")
  file.show("report.txt")
}