metodo de newton

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clc
N=10
x0=3
erro = 10^-3
function y=f(x) 
 // y=log(x-1)+cos(x-1)
 
 // y=exp(x)+2^-x+2,cos(x)-6
 
 // log(x-1)+cos(x-1)
 //y= 5*x^2+x-5
 y= x^3-x-4
endfunction
function y= fd(x)
 
 // y= 1/log(10)(x-1)
//y= exp(x)-log(2),2^x-2,sin(x)
 
 //y= 1/log(10),(x-1)-sin(x-1)
 // y= 10*x-1
 y=(3*x^2)-1
endfunction
xn = x0
for n= 1:N
 xn1 = xn - (f(xn)/ fd(xn))
 if abs((xn1-xn)/xn1) < erro then
 disp("Raiz é",xn1)
 return 
 end
 xn = xn1
 
end
printf("nao converge em n=%f interaçoes",N)