rm(list = setdiff(ls(), lsf.str())) load(paste(dirname(sys.frame(1)$ofile),"/problem_3_data.RData",sep = "")) input.model <- list() input.model$risk <-c(-1,"log_eut") input.model$H <- H input.model$Aeq <- Aeq input.model$beq <- 1 input.model$lb <- lb input.model$ub<- ub results <- rpsg_riskprog(input.model) print(results)