bissecao

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[x,inter]=bissecao(a,b,p,N)
 erro=1;
 inter=0;
 x=a;
 while(erro>p & cont<N)
 xa=x;
 x=(a+b)/2;
 if f(x)*f(a)<0 
 b=x;
 else
 a=x
 end
 erro=abs((x-xa)/x);
 inter=inter+1;
 end
endfunction
[x,inter]=bissecao(1100,1600,10^(-4),100)
// inter = 12;
// raiz = 1242.7979