x=[     0.38   1.45
       -1.08   2.81
        1.19   0.79
        0.05   2.25
        1.48   2.24
       -0.21   2.95
        1.53   0.43
       -1.51   1.09
        0.43   2.20
       -0.32   2.38
        1.09   2.39
       -1.03   1.40];


y=zeros(200,5);
k=1;
p=pi^(0.5);

for i=1:12,
for j=i+1:12,
     m1=x(i,1);  v1=abs(x(i,2));
     m2=x(j,1);  v2=abs( x(j,2));
     y(k,1)=x(i,1);
     y(k,2)=x(j,1);
     y(k,3)=x(i,2);
     y(k,4)=x(j,2);
if(m1<m2)
    t=m1;
    m1=m2;
    m2=t;
    t=v1;
    v1=v2;
    v2=t;
end

if((m1-m2)>(v1+v2)*p)
E4=0;
y(k,5)=0;  z(k)=4;
k=k+1;
end

 if((m1-m2)<abs((v2-v1))*p)
      if(v2>v1)
          c1=(v1*(m2-m1)+v1*(v1+v2)*p)/(v1+v2);
          c2=(v1*(m2-m1)+v1*(v2-v1)*p)/(v2-v1);
          c3=2*v1*p-(c1+c2);
          h1=((m2-m1)+(v1+v2)*p)/((v1+v2)*p);
          h2=((m2-m1)+(v2-v1)*p)/((v2-v1)*p);
          h3=h1+h2;
          E31=(c1*h1+c2*h2+c3*h3)/(2*(v1+v2)*p-(c1*h1+c2*h2+c3*h3));

            y(k,5)=E31; z(k)=31;
            k=k+1;
          end
      if(v1>v2)
          p=pi^(0.5);
          c2=(v2*(m2-m1)+v2*(v1+v2)*p)/(v1+v2);
          c1=(v2*(m2-m1)+v2*(v1-v2)*p)/(v1-v2);
          c3=2*v2*p-(c1+c2);
          h2=((m2-m1)+(v1+v2)*p)/((v1+v2)*p);
          h1=((m2-m1)+(v1-v2)*p)/((v1-v2)*p);
          h3=h1+h2;
          E32=(c1*h1+c2*h2+c3*h3)/(2*(v1+v2)*p-(c1*h1+c2*h2+c3*h3)) ;

            y(k,5)=E32; z(k)=32;
            k=k+1;
      end
end

if( (m1-m2)>p*abs(v1-v2) )
if( (m1-m2)<(v1+v2)*p )
     c1=(v1*(m2-m1)+v1*(v1+v2)*p)/(v1+v2);
     c2=(v2*(m2-m1)+v2*(v2+v1)*p)/(v2+v1);
     h1=((m2-m1)+(v1+v2)*p)/((v1+v2)*p);
     E2=(c1*h1+c2*h1)/(2*(v1+v2)*p-(c1*h1+c2*h1));

       y(k,5)=E2; z(k)=2;
       k=k+1;
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end
end

h=1;
for k=1:length(y),
   if(y(k,5)>0.8),
     r4(h,1)=y(k,1);
     r4(h,2)=y(k,2);
     r4(h,3)=y(k,3);
     r4(h,4)=y(k,4);
     r4(h,5)=y(k,5);
     r4(h,8)=z(k);
     h=h+1;
   end
end

for i=1:length(r4)
   if(r4(i,8)==31)
   r4(i,6)=(r4(i,1)+r4(i,2))/2;
   r4(i,7)=(r4(i,3)+r4(i,4))/2;
   end
   if(r4(i,8)==32)
   r4(i,6)=(r4(i,1)+r4(i,2))/2;
   r4(i,7)=(r4(i,3)+r4(i,4))/2;
   end
   if(r4(i,8)==2)
   r4(i,6)=(r4(i,1)+r4(i,2)+ p*(r4(i,4)-r4(i,3)))/2;
   r4(i,7)=(-r4(i,1)+r4(i,2)+ p*(r4(i,4)+r4(i,3)))/(2*p);
   end;
end