Programação

vanirialima
(usa Outra)

Enviado em 04/10/2017 - 17:11h

Boa tarde pessoal estou precisando de ajuda urgente para resolver este algoritmo,ficaria muito agradecida se alguém poder me ajudar,linguagem gnu/octave.
Aqui vai o link da questão: https://pt-static.z-dn.net/files/d99/9b89abdab3ad39f5acdcd208f543ed60.jpg
a maneira como fiz o programa no octave foi assim:
obs: metodo de newton
erro:0,0001 e x1:0.8

f(x)
function fv = f (x)
mf=70;r0=1030;r=0.45;
fv=mf-r0*pi*(1/3)*(2*r-h)^2*(r+h);
endfunction

f'(x)
function dfv = df (x)
mf=70;r0=1030;r=0.45;
dfv= r0*pi*(2*r-h)*h;
endfunction

e o programa para rodar:

function newton(xo,Niter,tol)
disp("")
disp ("Output for the Newton method")
disp("")
disp (" n x err f(x)")
for i=1:Niter
x=xo-f(xo)/df(xo);
if (f(xo)==0 |abs(x-xo)<tol)
disp("")
disp ("The method completed successfully!")
disp("")
return;
else
epsilon=abs(x-xo);
disp ([i-1, xo, epsilon, f(xo)]);
xo=x;
endif
endfor
disp("")
disp ("The method failed after (Niter)")
disp (Niter)
disp ("iterations")
disp("")
endfunction