用meshgridr = -1:0.01:0.15;k = 5:0.1:100[R,K] = meshgrid(r,k);m = ((108-8*R.^3+12*(81-12*R.^3.*K).^(1/2)).*K^2).^(1/3)max(m)或者用智能寻优算法
