clc
clear all
close all
pd = makedist('Normal', 'mu' , 0.5 , 'sigma' , 0.3 );
A = zeros(1000,1);
for i=1:1000
    random1= random(pd);
    random2= rand;
    random3= rand;
if random1 < 0.25
    if random2 < 0.25
        if random3 < 0.25
            A(i)=1;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=2;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=3;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=4;
        end
    elseif (0.25<=random2) && (0.5>random2)
        if random3 < 0.25
            A(i)=5;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=6;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=7;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=8;
        end
    elseif (0.5<=random2) && (0.75>random2)
        if random3 < 0.25 
            A(i)=9;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=10;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=11;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=12;
        end
    elseif (0.75<=random2) && (1>random2)
            if random3 < 0.25 
            A(i)=13;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=14;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=15;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=16;
        end
    end
    
elseif (0.25<=random1) && (0.5>random1)
    if random2 < 0.25
        if random3 < 0.25
            A(i)=17;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=18;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=19;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=20;
        end
    elseif (0.25<=random2) && (0.5>random2)
        if random3 < 0.25
            A(i)=21;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=22;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=23;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=24;
        end
    elseif (0.5<=random2) && (0.75>random2)
        if random3 < 0.25 
            A(i)=25;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=26;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=27;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=28;
        end
    elseif (0.75<=random2) && (1>random2)
            if random3 < 0.25 
            A(i)=29;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=30;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=31;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=32;
        end
    end
elseif (0.5<=random1) && (0.75>random1)
    if random2 < 0.25
        if random3 < 0.25
            A(i)=33;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=34;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=35;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=36;
        end
    elseif (0.25<=random2) && (0.5>random2)
        if random3 < 0.25
            A(i)=37;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=38;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=39;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=40;
        end
    elseif (0.5<=random2) && (0.75>random2)
        if random3 < 0.25 
            A(i)=41;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=42;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=43;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=44;
        end
    elseif (0.75<=random2) && (1>random2)
            if random3 < 0.25 
            A(i)=45;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=46;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=47;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=48;
        end
    end
elseif (0.75<=random1) && (1>random1)
    if random2 < 0.25
        if random3 < 0.25
            A(i)=49;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=50;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=51;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=52;
        end
    elseif (0.25<=random2) && (0.5>random2)
        if random3 < 0.25
            A(i)=53;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=54;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=55;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=56;
        end
    elseif (0.5<=random2) && (0.75>random2)
        if random3 < 0.25 
            A(i)=57;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=58;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=59;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=60;
        end
    elseif (0.75<=random2) && (1>random2)
            if random3 < 0.25 
            A(i)=61;    
        elseif (0.25<=random3) && (0.5>random3)
            A(i)=62;  
        elseif (0.5<=random3) && (0.75>random3)
            A(i)=63;  
        elseif (0.75<=random3) && (1>random3)
            A(i)=64;
        end
 
    end
end
end 
figure;
cdfplot(A) 
% figure;
% plot(A)