易拉罐模型的MATLAB优化

数学模型第五版第三章第九节的易拉罐模型求解
在满足b=0.103,kb=0.3,k1b=0.306,V2=364.8的情况下优化非线性规划,得到使得SV2取到最小值的情况下的r,r1,h,h1。

img

以下是我自己在MATLAB编的代码
r=optimvar('r','LowerBound',0);
h=optimvar('h','LowerBound',0);
r1=optimvar('r1','LowerBound',0);
h1=optimvar('h1','LowerBound',0);
b=0.103;
k=0.3/b;
k1=0.306/b;
prob=optimproblem;
prob.Objective=2pi.hb+pir1.^2kb+pi.r1.^2k1*b+pi.sqrt((r-r1).^2+h1.^2).(r+r1).*b;
prob.Constraints.cons1=pi.*r1.^2.*h+pi/3.h1.(r.^2+r1.^2+r.*r1)==364.8;
initialGuess.r=33.05;
initialGuess.r1=30.05;
initialGuess.h=110.5;
initialGuess.h1=5.05;
[sol,optval]=solve(prob,initialGuess);
sol.r;
sol.h;
sol.r1;
sol.h1;

但程序结果如下

img

求指点下,万分感谢