function l=icesn(x) % 保存为 icesn.m 文件
d=x(1);h=x(2);
p=d/(tan((pi/4)+2*atan(d/h)-(pi/2)));
l=(p/2)*(d*sqrt(p^2+d^2)+p^2*log(sqrt(p^2+d^2)+d));
-------------------------------------------------------------------------------------------
运行以下:
[x fval exitflag]=fminsearch('icesn',[2 8])
% x 为最小值点
%fval为最小值
% exitflag=1表示函数收敛于解;反之,=0