function x=gen(m,mu,sigma)
x=zeros(m,10);
for i=1:m
    e = 0.000001;
    cond=1;
    while cond
        % Original domain constraints:L2<=[x2,x3,x5,x6,x7,x9,x10]<=U2
        x(i,2)=normt1(mu(1,2),sigma(1,2),e,0.5);
        x(i,3)=normt1(mu(1,3),sigma(1,3),0.5,0.9);
        x(i,5)=normt1(mu(1,5),sigma(1,5),e,0.5);
        x(i,6)=normt1(mu(1,6),sigma(1,6),e,0.001);
        x(i,7)=normt1(mu(1,7),sigma(1,7),e,0.05);
        x(i,9)=normt1(mu(1,9),sigma(1,9),e,0.05);
        x(i,10)=normt1(mu(1,10),sigma(1,10),e,0.5);
        % New inequalities: L1<=[x1,x4,x8]<=U1
        x(i,1)=normt1(mu(1,1),sigma(1,1),e,0.05);
        x(i,4)=normt1(mu(1,4),sigma(1,4),e,0.005);
        x(i,8)=normt1(mu(1,8),sigma(1,8),e,0.1);
        % Original equalities:[x1,x4,x8]=b
        x(i,1) = 2 - (2*x(i,2) + 2*x(i,3) + x(i,6) +x(i,10));
        x(i,4) = 1 - (2*x(i,5) + x(i,6) + x(i,7));
        x(i,8) = 1 - (x(i,3) + x(i,7) + 2*x(i,9) + x(i,10));
        cond = any(x(i,:) < 0.000001);
    end   
    %x;
    %b1=x(i,1) + 2*x(i,2) + 2*x(i,3) + x(i,6) +x(i,10);
    %b2=x(i,4) + 2*x(i,5) + x(i,6) + x(i,7);
    %b3=x(i,8) + x(i,3) + x(i,7) + 2*x(i,9) + x(i,10);
end