Optimasi

Nama             : Sherly Taurin Siridion
NIM               : 20194920011
Email             : sherlytaurinsiri@gmail.com
RPubs            : https://rpubs.com/sherlytaurin/
Github           : https://github.com/sherlytaurin/
Jurusan          : Statistika
Address         : ARA Center, Matana University Tower
                         Jl. CBD Barat Kav, RT.1, Curug Sangereng, Kelapa Dua, Tangerang, Banten 15810.

Sebuah perusahaan memproduksi dua model kursi: A dan B. Model A membutuhkan 4 balok, 1 dudukan kursi, dan 1 sandaran kursi. Sedangkan model B membutuhkan 3 balok dan 1 dudukan kursi. Perusahaan memiliki stok awal 200 balok, 500 dudukan kursi, dan 100 sandaran kursi. Jika perusahaan membutuhkan lebih banyak balok, dudukan kursi, dan sandaran kursi, mereka dapat membeli kayu balok standar, yang harganya 80 euro per balok. Perusahaan dapat memproduksi 20 balok, 10 dudukan kursi, dan 2 sandaran kursi dari balok kayu standar. Biaya produksi model A adalah 30 euro/kursi, sedangkan biaya model B adalah 40 euro/kursi. Jika perusahaan menginformasikan bahwa jumlah minimum kursi yang harus diproduksi adalah 1000 unit per bulan. Tentukan model program linier, yang dapat meminimalkan total biaya.

Selesaikan permasalahan optimasi diatas dengan cara:

1. Grafik
2. Manual
3. Menggunakan R
4. Menggunakan Python

library(lpSolve)

# menentukan fungsi tujuan dan kendala
ftujuan <- c(30,40,80)

kendala_kiri <- matrix(c(4,3,-20,
                         1,1,-10,
                         1,0,-2,
                         1,1,0),nrow = 4, byrow = TRUE)

kendala_kanan <- c(200,500,100,1000)
arah_kendala <- c("<=","<=","<=",">=")

# mencari solusinya dengan lpsolve

solusi_r <- lp( direction = "min",
               objective.in = ftujuan,
               const.mat = kendala_kiri,
               const.rhs = kendala_kanan,
               const.dir = arah_kendala,
               all.int = T,
               compute.sens = TRUE)

solusi_terbaik <- solusi_r$solution

#Menampilkan nilai A,B,c
names(solusi_terbaik) <- c('A' , 'B', 'c')

print(solusi_terbaik)

##   A   B   c 
## 420 580 161

# menampilkan total cost
print(paste("Total cost:", solusi_r$objval, sep=""))

## [1] "Total cost:48680"