Fuzing_Conditions_Dynamic_Case_I.m

function [detonate , missDistance ] = fuze( xx )
% returns detonate = 1 and miss distance -----> fuzing conditions have been met.
% returns detonate = 0 and miss distance = Nan ------> fuzing conditions haven't been met

R = sqrt( ( xx(:,7)-xx(:,3) ).^2 + ( xx(:,8)-xx(:,4) ).^2  );
temp1a = ( abs(R(2:end,:)) < 10); %Is R<10 in intersample?
temp1b = ( R(2:end,:).*R(1:end-1,:) < 0) ; % Is R crossing 0 in intersample?
temp1c = ( R(2:end,:)-R(1:end-1,:) > 0) ; % Is R increasing in intersample?
temp1bc = any( [temp1b temp1c] , 2 );
temp1 = all( [temp1a temp1bc] , 2 ); %Is R<10 AND R crossing 0 or Rdot>0 in intersample?
temp2 = xx(2:end,3) < 0; %Did pursuer hit the ground in intersample?
temp3 = xx(2:end,7) < 0; %Did evader hit the ground in intersample?  
temp4 = xx(2:end,1) < 0; %Did pursuer slow to 0 speed in intersample?
temp5 = xx(2:end,5) < 0; %Did evader slow to 0 speed in intersample?
temp = any( [temp1 temp2 temp3 temp4 temp5] , 2); %Is any fuzing condition satisfied in intersample?
if any(temp,1)
    fuze_index = 0;
    for kk = 1:length(temp)
        if temp(kk) == 1
            fuze_index = kk;
            break
        end
    end
    missDistance = R(fuze_index,:);
    detonate = 1;
else
    missDistance = NaN;
    detonate = 0;
end

end