library(lpSolve)
## Warning: package 'lpSolve' was built under R version 4.3.3
z=c(12,8)
x=matrix(c(3,7,
         2,4,
         2,4), nrow = 3, byrow = TRUE)
dir=c('<=','<=','<=')
rhs=c(42,30,48)
lp('max',z,x,dir,rhs,compute.sens=TRUE)
## Success: the objective function is 168
lp('max',z,x,dir,rhs,compute.sens=TRUE)$solution
## [1] 14  0
lp('max',z,x,dir,rhs,compute.sens=TRUE)$sens.coef.from
## [1]  3.428571e+00 -1.000000e+30
lp('max',z,x,dir,rhs,compute.sens=TRUE)$sens.coef.to
## [1] 1.0e+30 2.8e+01
print("berikut hasil dual :")
## [1] "berikut hasil dual :"
lp('max',z,x,dir,rhs,compute.sens=TRUE)$duals
## [1]   4   0   0   0 -20