% 体積が一定で表面積が最小となる直方体を求める

a = 5;

cond = @(q) cond_param(q,a);
qinit = [1;1;1];
[qopt,fopt] = multiplier_method(@area, cond, @ineq, qinit);

qopt, fopt
