matlab用二分法求E(i) = Q1/4piAr^2-Q2/4piA(1-r)^2零点的代码

很高兴能为您解答问题。下面是用二分法求解E(i)函数零点的MATLAB代码：

function [r] = find_root(Q1, Q2, A)
% Q1: Q1/4piA
% Q2: Q2/4piA
% A: Area of cylindrical electrode
low = 0;
high = 1;
while (high - low) > 1e-8
    mid = (low + high)/2;
    if (Q1/(4*pi*A*mid^2) - Q2/(4*pi*A*(1-mid)^2)) > 0
        high = mid;
    else
        low = mid;
    end
end
r = (low + high)/2;
end

希望上面的代码对您有所帮助。如果有其他问题，欢迎随时提问。