Resolve o problema de curva de perseguição

// Resolve o problema de curva de perseguição

// A posição do alvo é (x_a,y_a) e do perseguidor (x_p,y_p)

// Considera que o alvo encontra-se em movimento circular uniforme com

// velocidade angular omega=10 e raio=10

// O perseguidor está a uma velocidade V=50

// As posições iniciais são: x_a=10, y_a=0, x_p=100, y_p=0

// Associa t, x_a, y_a, x_p, y_p aos elementos 1, 2, 3, 4 e 5 do vetor x

X0=[0 ; 10 ; 0 ; 100 ; 0]

Y=rk4(curva_perseguicao,X0,0.001,100)

clf

fig0=figure(0)

clf

fig0.background=8

plot(Y(:,2),Y(:,3),'r-',Y(:,4),Y(:,5),'b-')

mtlb_axis('equal')

xlabel('x')

ylabel('y')

set(gca(),"grid",[1 1])

fig1=figure(1)

clf

fig1.background=8

plot(Y(:,1),Y(:,2)-Y(:,4),'r-',Y(:,1),Y(:,3)-Y(:,5),'g-')

legend('x_a-x_p','y_a-y_p')

xlabel('t')

ylabel('diferença de coordenadas')

set(gca(),"grid",[1 1])