ALGORITMO PARA SECANTE
r=input('Digite el radio de la bala:');
x0=input('Digite el extremo izquierdo del intervalo:');
x1=input('Digite el extremo derecho del intervalo:');
    tol=0.00001;
    j=2;
    i=1;
    x(1)=x0;
    x(2)=x1;
    ea(i)=100;
    while abs(ea(i))>=tol,
       x(j+1)=x(j)-f(x(j),r)*(x(j-1)-x(j))/(f(x(j-1),r)-f(x(j),r));
       ea(i+1)=(x(j+1)-x(j))/x(j+1)*100;
       j=j+1;
       i=i+1;
    end
    fprintf(' i     x(i)         Error aprox (i) \n');
    fprintf('%2d \t %11.7f \t \n',0,x(1));
    for k=2:j;
        fprintf('%2d \t %11.7f \t %7.3f \n',k-1,x(k),ea(k-1));
    end
disp(['Velocidad terminal=',num2str(x(j))])

Función Usada Dentro De Secante
function y=f(a,r)
    y=(600.1*(r^2)*nthroot(a^5,2))+(0.942*(a^2)*(r^2))-0.0196;
end