% Métods de Aproximação em Engenharia - Mst. Eng. Construção
% Carlos Balsa - 2009/2010
% ficha nº 6 ex. 1

theta=5/18*pi;
g=9.81;
y0=1;
v0=25;

a=tan(theta);
b=g/(2*v0^2*cos(theta)^2);

x=0:1:60;
y=a.*x-b.*x.^2+y0;
dy=a-2.*b.*x;

plot(x,y,'b'); hold on;
plot(x,dy,'r');
legend('y','dy');
grid on;


g=@(x) a-2*b*x;
dg=@(x) -2*b;
x0=20;
maxit=10;
tol=0.01;
[x,res,k]=nle_newtraph(g,dg,x0,maxit,tol)
y=a*x-b*x^2+y0
