Falsa Posição

Conteúdo do arquivo

%             -------------                             ------------
%             function file                             falsaPosicao.m
%             -------------                             ------------

        function  [x] = falsaPosicao(a,b,epsilon)
        %====================================================
        format long
        max =  50;
        fa = (a)^5 + (-10/9)*a^3 +(5/21)*a;
        fb = (b)^5 + (-10/9)*b^3 +(5/21)*b;
                
        for k=1:max
                       
            x(k)=(a*fb - b*fa)/(fb-fa);
                        
            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);
               fb=fxk;
            else 
               a = x(k);
               fa = fxk;
            end 
                       
        end
        %====================================================