rm(list = setdiff(ls(), lsf.str()))
load(paste(dirname(sys.frame(1)$ofile),"/problem_2_data.RData",sep = ""))

input.model <- list()
input.model$w <- c(-10,-3)
input.model$d <- d
input.model$risk <-c(1.3571,"pm2_pen",-1.3571,"pm2_pen")
input.model$H[[1]] <- H1
input.model$H[[2]] <- H2
input.model$Aeq <- Aeq
input.model$beq <- 1
input.model$lb <- lb
input.model$ub<- ub

results <- rpsg_riskprog(input.model)

print(results)