Método da Secante em Scilab

-
Versão do Scilab: 5.5.2.

 clear 
clc
//METODO DA SECANTE
function [resultado] = equacao(x)
    resultado = %e^(x/2)-x^2             //aqui vai a equacao original
endfunction

 //INTERVALO CONTENDO AS S RAIZES
x0 = 1; //intervalo inicial
x1 = 2 ; //intervlo final

 //Erro Tolerável
erro=5*10^(-5)                 //erro
diferenca=2*erro;           //Condição de parada do 
            
contIteracao = 0;           //Contador de Iterações

 while diferenca > erro
    fx0 = equacao(x0) 
    fx1 = equacao(x1)  
    x = (x0*fx1 - x1*fx0)/(fx1 - fx0);
    //printf("%g",x)
    //printf("\n")
    fx = equacao(x)
    diferenca = abs(x-x1);
    x0 = x1
    x1 = x
                
    contIteracao = contIteracao +1;
end

 disp("Raiz aproximada da funcao 1: ")
disp(x0)
disp("Numero de iterações: ")
disp(contIteracao)





Comentários

Postar um comentário

Obrigado pelo contato, retornarei o mais breve possível!

Postagens mais visitadas deste blog

Método da Interpolação por Eliminação Gaussiana em Scilab

-
Versão do Scilab: 5.5.2 Observação: Só calcula equação linear e quadrática; clear clc //Eliminação Gaussiana //Variavel matriz é a matriz incompleta //Variavel matrizResposta é a matriz que contem as respostas dos sistemas lineares function  matrizCoeficientes= eliminacaoGauss(matriz, matrizResposta)     matrizCompleta=[matriz matrizResposta];     n=size(matriz,1);     disp("Resolução do Sistema")     disp("Matriz Ampliada = ")     disp(matrizCompleta)     for j=1:(n-1)         for i=(j+1):n            mij=matrizCompleta(i,j)/matrizCompleta(j,j)            matrizCompleta(i,:)=matrizCompleta(i,:)-mij*matrizCompleta(j,:)         end     end     disp("Matriz Escalonada")     disp(matrizCompleta)     //Fazer Retrosubstituição     matrizCoeficientes(n)=matr...
Leia mais

Cálculo do Erro para Interpolação em Scilab

-
Versão do Scilab: 5.5.2 clear clc function [orderValues] = getOrder(lastOrder, vectorX, orderNumber)     for i = 1 : size(lastOrder, 2) - 1 // size 1 = quantas linhas, size 2 = quantas colunas        orderValues(1,i) = (lastOrder(1,i+1) - lastOrder(1,i)) / (vectorX(1,orderNumber+i-1) - vectorX(1,i))     end endfunction function [f, polyF] = newtonInterpol(vectorY, vectorX)     lastOrder = vectorY     polyF = poly([vectorY(1,1)],"x", "coeff")     for i = 2 : size(vectorX, 2) // size 1 = quantas linhas, size 2 = quantas colunas         lastOrder = getOrder(lastOrder, vectorX, i)         polyA = poly([-vectorX(1,1) 1],"x","coeff")         for j = 2 : i-1             polyA = polyA * poly([-vectorX(1,j) 1],"x","coeff")         end         polyA = polyA * lastOrder(1,1) ...
Leia mais

Método de Interpolação de Lagrange em Scilab

-
Versão do Scilab : 5.5.2 clc  clear //funcao que retorna o polinomio de lagrange function y=polinomioLagrange(tab); // calcula o polinomio interpolador pelo metodo de lagrange;     [l,n]=size(tab);     x=tab(1,:);     for i=1:n           if i>size(x) then            error(" i incompativel com dimensao da tabela");         end;         [l,n]=size(x);         for k=1:i-1             r(k)=x(k)         end;         for k=i:n-1              r(k)=x(k+1)         end;         P=1;         for k=1:(n-1)              P=P*(x(i)-r(k))         end;         if P==0 then    ...
Leia mais