// Mtodo de Newton-Raphson - NR
// Encontra uma raiz aproximada de uma determinada funo
deff('y=f(x)', 'y=x^3-5*x^2+x+3');
// funo que deseja encontrar a raiz
deff('y=f1(x)', 'y=3*x^2-10*x+1');
// derivada da funo
deff('y=fi(x)','y=x-(f(x)/f1(x))');
//funo de iterao
x0=4.2;
//valor inicial de x
epsilon=0.0000000001;
//tolerncia (preciso desejada)
Nmax=100;
//nmero mximo de iteraes
k=0;
//contador de iteraes
x=fi(x0);
while (abs(x-x0)>epsilon) & (k<Nmax)
x0=x;
k=k+1;
x=fi(x0);
disp(x0)
disp(x)
end;
disp(x) //Comando para exibir o valor de x
disp(k+1) //Comando para exibir o nmero de iteraes executadas
//Mtodo da Secante
//Encontra uma raiz aproximada de uma determinada funo
deff('y=f(x)', 'y=log(x)-sin(x)') //funo da qual se deseja encontrar a raiz
xa=1.5; //valor inicial de xa
xb=1.7; // valor inicial de xb
epsilon=0.0001; // tolerncia (preciso desejada)
Nmax=100; // Nmero mximo de iteraes
k=0;
while (abs(xb-xa/xb)>epsilon)&(k<Nmax)
k=k+1;
xm=(xa*f(xb)-xb*f(xa))/(f(xb)-f(xa));
xa=xb;
xb=xm;
disp(xm)
end
disp(k)