function tmax = show_raptorchase(theta)

[tmax,sol] = raptorchase(theta);

positions = deval(sol,sol.x);

nraptors = size(positions,1)/2 - 1;
assert(ceil(nraptors) == nraptors) % integer operations are exact

cla;
hold on;
axis off;
plot(positions(1,:),positions(2,:),'k-', 'LineWidth', 3);
for i=1:nraptors
    plot(positions(2*i+1,:),positions(2*i+2,:),'r-', 'LineWidth', 3);
end
axis square;
hold off;

title(sprintf('Surival time = %.3f sec', tmax),'FontSize',20);
