Msecante

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

function y=f(x)
 // y = tan(teta)*x - (g / 2*(vo^2)*(cos(teta)^2))*x^2 + (y - y0)
 y= tan(0.1862)*(347) - (2982.24/(2*(x^2)*cos(0.1862)^2))*347^2 + 55;
endfunction
function [raiz,iter]=secante(x0,x1,p,N)
 ra1 = x0; //raiz antiga 1
 rn = x1; //raiz nova
 cont = 0;
 while(%T)
 ra2 = ra1;
 ra1 = rn;
 rn = ra1 - f(ra1)*(ra2 - ra1)/(f(ra2)-f(ra1));
 cont = cont + 1;
 if (abs(rn - ra1)/abs(rn)<p) then
 break
 else if (cont>N) then
 break
 end
 end
 end
 raiz = rn;
 iter = cont;
endfunction
[raiz,iter] = secante(1100,1600,10^(-4),100)
disp(raiz,iter)
//iter = 5.
//v0 = 1242.801.