Vx = j(nrow(jari),1); /*matrix initialization*/
Vy = j(nrow(jari),1); /*matrix initialization*/
do i = 1 to nrow(jari);
    /* velocity vector (vx,vy)*/
  Vx[i]= velo[i] * cosm[i];
    Vy[i]= velo[i] * sinm[i];
end;
print Vx, Vy;

/* intersection point*/
Xpcij = j(nrow(m), nrow(m));
Ypcij = j(nrow(m), nrow(m));
do i = 1 to nrow(m);
  do j=1  to nrow(m);
      Xpcij[i,j] = (koordinatY[j,]-((Vy[j,]/Vx[j,])*koordinatX[j,])-koordinatY[i,]-((Vy[i,]/Vx[i,])*koordinatX[i,]))/((Vy[i,]/Vx[i,])-(Vy[j,]/Vx[j,]));
      Ypcij[i,j] = ((Vy[i,]/Vx[i,])*(koordinatY[j,]-((Vy[j,]/Vx[j,])*koordinatX[j,]))-(Vy[j,]/Vx[j,])*(koordinatY[i,]-((Vy[i,]/Vx[i,])*koordinatX[i,])))/((Vy[i,]/Vx[i,])-(Vy[j,]/Vx[j,]));
  end;
end;
print Xpcij;
print Ypcij;

/* false conflict criteria*/
d1ij = j(nrow(m), nrow(m));
d1ji = j(nrow(m), nrow(m));
d2ij = j(nrow(m), nrow(m));
d2ji = j(nrow(m), nrow(m));
do i = 1 to nrow(m);
  do j=1  to nrow(m);
      d1ij[i,j] = sqrt(((koordinatX[i,]-Xpcij[i,j])**2)+((koordinatY[i,]-Ypcij[i,j])**2));
      d1ji[i,j] = sqrt(((koordinatX[j,]-Xpcij[i,j])**2)+((koordinatY[j,]-Ypcij[i,j])**2));
      d2ij[i,j] = sqrt((((koordinatX[i,]+cosm[i,])-Xpcij[i,j])**2)+(((koordinatY[i,]+sinm[i,])-Ypcij[i,j])**2));
      d2ji[i,j] = sqrt((((koordinatX[j,]+cosm[j,])-Xpcij[i,j])**2)+(((koordinatY[j,]+sinm[j,])-Ypcij[i,j])**2));
  end;
end;
print d1ij;
print d1ji;
print d2ij;
print d2ji;

di= d2ij-d1ij;
dj= d2ji-d1ji;
print di,dj;

/*False Conflict*/
FC = j(nrow(m), nrow(m));
do i = 1 to nrow(m);
  do j=1  to nrow(m);
        if di[i,j] ^=missing(di[i,j]) | dj[i,j] ^= missing(dj[i,j]) then
        do;
           if di[i,j]>0 & dj[i,j]>0 then FC[i,j]= 1;
             else FC[i,j] =0;
          end;
          else
        do;
             FC[i,j] =0;
        end;
  end;
end;
print FC;