Bissecção

Conteúdo do arquivo

%             -------------                             ------------
%             function file                             bisseccao.m
%             -------------                             ------------

        function  [x] = bisseccao(a,b,epsilon)
        %=================================================
        format long
        max=50;
        fa = (a)^5 + (-10/9)*a^3 +(5/21)*a;

        for k=1:max
            
            x(k)=(b+a)*0.5;
            
            fxk = (x(k))^5 + (-10/9)*x(k)^3 +(5/21)*x(k);
            
            errox=b-a;
         
            errof=abs(fxk);
            
            if (errox < epsilon) & (errof < epsilon)
               break
            end    
            
            if ( (fxk*fa) < 0.0 )
               b = x(k);
            else 
               a = x(k);
               fa = fxk;
            end 
                       
        end
        %=================================================