rm(list = setdiff(ls(), lsf.str())) load(paste(dirname(sys.frame(1)$ofile),"/problem_partial_moment_penalty_data.RData",sep = "")) input.model <- list() input.model$risk <-c(2.875,"pm_pen", 0.25,"pm_pen_g", 3.9167,"pm_pen", 0.25,"pm_pen_g", 6,"pm_pen", 0.25,"pm_pen_g", 12.25,"pm_pen", 0.25,"pm_pen_g") input.model$w <- c(0,0,0,0,0,0,0,0) input.model$H <- H input.model$Aeq <- Aeq input.model$beq<- 0 results <- rpsg_riskprog(input.model) print(results)