function [out]=soft_decision(data1,data2);
d(1)=(abs(data1-1/sqrt(10)*(-3-j*3)))^2+(abs(data2-1/sqrt(10)*(1+j)))^2;
d(2)=(abs(data1-1/sqrt(10)*(-3-j)))^2+(abs(data2-1/sqrt(10)*(1-j*3)))^2;
d(3)=(abs(data1-1/sqrt(10)*(-3+j)))^2+(abs(data2-1/sqrt(10)*(1+j*3)))^2;
d(4)=(abs(data1-1/sqrt(10)*(-3+j*3)))^2+(abs(data2-1/sqrt(10)*(1-j)))^2;
d(5)=(abs(data1-1/sqrt(10)*(-1-j*3)))^2+(abs(data2-1/sqrt(10)*(-3+j)))^2;
d(6)=(abs(data1-1/sqrt(10)*(-1-j)))^2+(abs(data2-1/sqrt(10)*(-3-j*3)))^2;
d(7)=(abs(data1-1/sqrt(10)*(-1+j)))^2+(abs(data2-1/sqrt(10)*(-3+j*3)))^2;
d(8)=(abs(data1-1/sqrt(10)*(-1+j*3)))^2+(abs(data2-1/sqrt(10)*(-3-j)))^2;
d(9)=(abs(data1-1/sqrt(10)*(1-j*3)))^2+(abs(data2-1/sqrt(10)*(3+j)))^2;
d(10)=(abs(data1-1/sqrt(10)*(1-j)))^2+(abs(data2-1/sqrt(10)*(3-j*3)))^2;
d(11)=(abs(data1-1/sqrt(10)*(1+j)))^2+(abs(data2-1/sqrt(10)*(3+j*3)))^2;
d(12)=(abs(data1-1/sqrt(10)*(1+j*3)))^2+(abs(data2-1/sqrt(10)*(3-j)))^2;
d(13)=(abs(data1-1/sqrt(10)*(3-j*3)))^2+(abs(data2-1/sqrt(10)*(-1+j)))^2;
d(14)=(abs(data1-1/sqrt(10)*(3-j)))^2+(abs(data2-1/sqrt(10)*(-1-j*3)))^2;
d(15)=(abs(data1-1/sqrt(10)*(3+j)))^2+(abs(data2-1/sqrt(10)*(-1+j*3)))^2;
d(16)=(abs(data1-1/sqrt(10)*(3+j*3)))^2+(abs(data2-1/sqrt(10)*(-1-j)))^2;

[distance,position]=min(d);

switch position
    case 1
        temp=[0 0 0 0];
    case 2
        temp=[0 0 0 1];
    case 3
        temp=[0 0 1 0];
    case 4
        temp=[0 0 1 1];
    case 5
        temp=[0 1 0 0];
    case 6
        temp=[0 1 0 1];
    case 7
        temp=[0 1 1 0];
    case 8
        temp=[0 1 1 1];
    case 9
        temp=[1 0 0 0];
    case 10
        temp=[1 0 0 1];
    case 11
        temp=[1 0 1 0];
    case 12
        temp=[1 0 1 1];
    case 13
        temp=[1 1 0 0];
    case 14
        temp=[1 1 0 1];
    case 15
        temp=[1 1 1 0];
    case 16
        temp=[1 1 1 1];
end
out=temp;