método da secante

Para responder essa pergunta devemos colocar em prática nosso conhecimento sobre Cálculo Numérico.

% Dado uma funçao f(x)=f e o C.P.=<epsolon e gama calcule(aproxime o zero(raiz))de f(x):

f=input ('Digite a expressao f(x):' );

x0=input ('Digite o intervalo inicial: ');

x1=input ('Digite o intervalo final: ');

erro=input('Digite o valor do Erro: ');

f1=inline(f);

x2=x1-(f1(x1)*(x1-x0)/(f1(x1)-f1(x0)));

m=abs(x1-x0);

while (m>=erro)

x0=x1;

x1=x2;

x2=x1-(f1(x1)*(x1-x0)/(f1(x1)-f1(x0)));

m=abs(x1-x0);

end

fprintf('%d',x2);

%plotagem do grafico

vi=input('Digite intervalo inicial para o grafico: ')

vf=input('Digite o intervalo final para o grafico: ')

y=[vi: 0.1 :vf];

y1=inline(f);

fplot(f,[vi vf]);

axis([vi vf vi vf]);

grid('on');

zoom(1);

legend (f)

xlabel ('Eixo X')

ylabel ('Eixo Y')

Portanto, o algoritmo apresentado é o código para resolução do método da secante para matlab, assim como solicitado.

interessante

Para responder essa pergunta devemos colocar em prática nosso conhecimento sobre Cálculo Numérico.

% Dado uma funçao f(x)=f e o C.P.=<epsolon e gama calcule(aproxime o zero(raiz))de f(x):

f=input ('Digite a expressao f(x):' );

x0=input ('Digite o intervalo inicial: ');

x1=input ('Digite o intervalo final: ');

erro=input('Digite o valor do Erro: ');

f1=inline(f);

x2=x1-(f1(x1)*(x1-x0)/(f1(x1)-f1(x0)));

m=abs(x1-x0);

while (m>=erro)

x0=x1;

x1=x2;

x2=x1-(f1(x1)*(x1-x0)/(f1(x1)-f1(x0)));

m=abs(x1-x0);

end

fprintf('%d',x2);

%plotagem do grafico

vi=input('Digite intervalo inicial para o grafico: ')

vf=input('Digite o intervalo final para o grafico: ')

y=[vi: 0.1 :vf];

y1=inline(f);

fplot(f,[vi vf]);

axis([vi vf vi vf]);

grid('on');

zoom(1);

legend (f)

xlabel ('Eixo X')

ylabel ('Eixo Y')

Portanto, o algoritmo apresentado é o código para resolução do método da secante para matlab, assim como solicitado.