4x+y<= 9, x+2y >= 4, 2x-3y >= -6の時、2x+y、x^2+y^2の最大値と最小値「2010 京都大学 文系【２】」(その1) RとWolframAlphaとCVXPY(非線型計画法)でやってみた。

旧タイトル 2021年03月30日
RとWolframAlphaとCVXPYで(非線型計画法)「2010年度 京都大・文系数学第2問」 をやってみた。

オリジナル

・大学入試数学問題集成 様＞ 文系【２】 テキスト

・(参考) リンク切れ???
2010年度 京都大・文系数学第2問 問題と解答

https://www.densu.jp/kyoto/10kyotolpass.pdf#page=2

https://www.densu.jp/kyoto/10kyotolpass.pdf#page=7

vbaで

スタック・オーバーフローから削除されました。
https://ja.stackoverflow.com/questions/75328/excelのソルバーをvbaで自動化する方法-最小値をやってみた

Rで

(参考)

library(lpSolve)
f.obj <- c(2, 1)
f.con <- matrix (c(4, 1, 1, 2, 2, -3), ncol=2, byrow=TRUE) 
f.dir <- c("<=", ">=", ">=")
f.rhs <- c(9, 4, -6)
lp ("max", f.obj, f.con, f.dir, f.rhs)
lp ("max", f.obj, f.con, f.dir, f.rhs)$solution 
lp ("min", f.obj, f.con, f.dir, f.rhs)
lp ("min", f.obj, f.con, f.dir, f.rhs)$solution
# > lp ("max", f.obj, f.con, f.dir, f.rhs)
# Success: the objective function is 6 
# > lp ("max", f.obj, f.con, f.dir, f.rhs)$solution 
# [1] 1.5 3.0
# > lp ("min", f.obj, f.con, f.dir, f.rhs)
# Success: the objective function is 2 
#> lp ("min", f.obj, f.con, f.dir, f.rhs)$solution 
# [1] 0 2

WolframAlphaで

最大値

最小値

CVXPYで

最大値
最小値

いつもの? sympyの実行環境と 参考のおすすめです。

(テンプレート)

いつもと違うおすすめです。