function rlc_f_bisection C1 = 2e-6; C2 = 7e-6; a = C1; b = C2; for i=1:100 c =(a+b)/2; if f(a) * f(c) < 0 b = c; else a = c; end if((abs(f(a)) < 1e-7) || (abs(f(b)) < 1e-7)) break; end a b i end end function y = f(x) f0 = 500; L = 20e-3; y = f0 - 1./(2*pi.*sqrt(L*x)); end