% 未定乗数を含む式を最小化
% 未定乗数を含む式を最小化することで体積が一定で表面積が最小となる直方体を求める
qinit = [ 0; 0; 0; 0 ];  % initial value
options = optimset('Display','iter');
[xmin, fmin] = fminsearch(@area_multiplier_volume, qinit, options);
xmin
fmin
