Formulas Electricidad

Hola a todos. Este es uno de mis primeros tiempos en realidad la programación en MATLAB, y mi maestro no ofrece mucha ayuda con la parte real de la programación. Tengo un algoritmo Bisecar que he creado, pero sólo se ejecuta el bucle una vez, y nunca llega a la respuesta final. También tengo ni idea de cómo representar los resultados con un barras de error. ¿Puede alguien ayudarme?

============================================= =========

función [rm, eqrm, err, eqrlow, eqrup] = Entradas % resistencia (func) nombre de la función % rlow dejó de punto final % rup derecho variable % delta tolerancia convergencia % Retorno solución% rm: la raíz solución eqrm%: el valor de la función estimación de error% err en la raíz % ------ -------------------------------------------------- ------------------ eqrlow = feval (func, rlow); eqrup = feval (func, RUP); rm = (rlow + RUP) / 2; eqrm = feval ( func, rm); errar = abs (RUP-rlow) / 2; si eqrlow * eqrup> 0, vacaciones, fin max1 = 1 + redonda ((log (RUP-rlow) -log (delta)) / log (2) ); para k = 1: max1, rm = (rlow + RUP) / 2; eqrm = feval (func, rm); si eqrm == 0, rlow = rm; RUP = rm; elseif eqrup * eqrm> 0, rup = rm; eqrup = eqrm; otra cosa rlow = rm; eqrlow = eqrm; final si rup-rlow <delta, descanso, extremo extremo rm = (rlow + RUP) / 2; eqrm = feval (func, rm); errar = abs (RUP-rlow) / 2; ========================================== ========== función eq = resistencia (r) r = 0; % Teniendo en cuenta: c = 0, L = 5, y (q / q0) = 0,01 eq = [exp (-0,005 * r) ] * [cos (sqrt (2000- (0,01 * (r ^ 2))) * 0,05)] - 0,01;

...