fid = fopen('test_001.img','w');
r=random('Normal',100,0,64,64);
z=random('Uniform',0,5,64,64);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a=14;
b=16;
for i=15:18
    for j=a:b
        r(i,j)=random('Normal',102,0,1,1);
    end
    a = a-1;
    b = b+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=15:18
    for j=31:34
        r(i,j)=random('Normal',102,0,1,1);
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=14:19
    for j=47:48
         r(i,j)=random('Normal',102,0,1,1);
    end
end

for i=16:17
    for j=45:50
        r(i,j)=random('Normal',102,0,1,1);
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

a=13;
b=17;
for i=30:35
    for j=a:b
       r(i,j)=random('Normal',102,0,1,1); 
    end
    
    a=a-1;
    b=b+1;
    
end 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=30:35
    for j=30:35
        r(i,j)=random('Normal',102,0,1,1); 
    end    
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i =29:35
    for j=47:49
         r(i,j)=random('Normal',102,0,1,1);
    end 
end 

for i=31:33
    for j=45:51
        r(i,j)=random('Normal',102,0,1,1);
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

a=13;
b=18;
for i=46:53
    for j=a:b
        r(i,j)=random('Normal',102,0,1,1);
    end
    a=a-1;
    b=b+1;
end 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=45:53
    for j=29:37
        r(i,j)=random('Normal',102,0,1,1);
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=44:53
    for j=46:49
        r(i,j)=random('Normal',102,0,1,1);
    end
end 

for i=47:50
   for j=43:52
        r(i,j)=random('Normal',102,0,1,1);
   end
end 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=r+z;count =fwrite(fid,r);
fclose('all');

fid = fopen('test_002.img','w');
s=random('Normal',100,0,64,64);z=random('Uniform',0,5,64,64);
for i=8:12
for j=8:12
s(i,j)=random('Normal',102,0,1,1);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
s=s+z;count=fwrite(fid,s);
fclose('all');