D1 = (min1:0.001:max1)/2;D2 = (min2:0.001:max2)/2;V = [];for i = 1:size(D1,2) for j = 1:size(D2,2) V = [Vm πL(D1^2＋D1*D2＋D2^2)/3]; endendVmin = min(min(V))
不是有函数嘛
