新建函数文件ff.m：
function y = ff(x) f = [.0187 .0474 inf]; A = [1.05e-5 3.56e-3 1.98e-2]; K = [3.4 1.93 1.36]; t = arrayfun(@(x)find(x y = A(t).*x.^-K(t);end
主程序：
t = linspace(.01, .153, 100);plot(t, ff(t))
得出结果：
你所画的图有错，f^-k是凹函数。
