x=[
   1.0839     3.5188
  -1.0303     7.4276
  -1.0917     8.4593
  -5.6955     3.3767
  -5.5765    -8.3346
  -2.0102     5.2161
  -2.3441     7.0289
  -2.7482     6.3115
   1.9131   -14.1712
   4.1264     2.9529
   0.8063   109.0353
   1.5041     0.1678
   2.5590     2.6079
   0.1645     6.8001
   1.4596     3.6431
   1.3348     0.3024
   1.7751     2.2122
   0.4185     0.6944
  -0.2041     0.9345
   4.1705  2485.6069
  -0.0245     8.7036
  -5.2400     3.9237
   0.3409     2.4973
  -2.2344     0.8467
   0.7570     1.0680 ];

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

for i=1:25,
for j=i+1:25,
     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=v1;
end

if((m1-m2)>(v1+v2)*p)
E4=0;
y(k,5)=0;
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;
            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;
            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;
       k=k+1;
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end
end

r=zeros(50,3);
h=1;
for k=1:length(y),
   if(y(k,5)>0.8),
     r(h,1)=y(k,1);
     r(h,2)=y(k,2);
     r(h,3)=y(k,3);
     r(h,4)=y(k,4);
     r(h,5)=y(k,5);
     h=h+1;
   end
end