MATLAB 求取離散點的曲率最大值

x0 = linspace(0.1,2,100);%x0,y0驗證函數離散點，可以非等間隔
y0 = 1./x0;
h1 = abs(diff(x0)) ;
h = [h1 h1(end)];
ht = h;
yapp1 = gradient(y0)./ht; %matlab數值近似
yapp2 = del2(y0)./ht; %matlab數值近似
k2 = abs(yapp2)./(1+yapp1.^2).^(3/2);
figure
plot(k2)
title('曲率曲線')
[~,maxFlag] = max(k2);%曲率最大位置
x_max = x0(maxFlag);
y_max = y0(maxFlag);
%畫出圖像 標註曲率最大點
figure
plot(x0,y0,'.-');
hold on;
plot(x_max,y_max,'rp')
title('標註最大麴率點')
xlabel('log10((norm(B*Xk-L)))')
ylabel('log10((norm(Xk)))')