clear all; close all; clc
%function m=Map;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  Start @ location (12,12). Apply controls to reach (0,0)  %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xw=-14; xe=14; ys=-14; yn=14; %map dimensions, west, east, south, north
mrx=0.05; mry=0.05; % map resolution {x,y}
[xm ym]=meshgrid(xw:mrx:xe,ys:mry:yn);

m=0*xm.*ym;  % initialize the map m

for i=1:size(xm,2)
    for j=1:size(ym,1)
        r=sqrt(xm(1,i)^2+ym(j,1)^2);
        ang=atan2(ym(j,1),xm(1,i))+pi;
        if 10<=r && r<=11 && (pi/2<=ang && ang<=2*pi) 
            m(j,i)=1;
        elseif  6<=r && r<=7 && ((0<=ang && ang<=7*pi/4) || (13*pi/7<=ang && ang<=2*pi))
            m(j,i)=1;
        elseif  3<=r && r<=4 && ~(pi/2<=ang && ang<=pi)
            m(j,i)=1;
        elseif (abs(xm(1,i)+ym(j,1)))<=0.5 && -4<ym(j,1)*xm(1,i) && 0<xm(1,i) && ym(j,1)<0
            m(j,i)=1;
        elseif xm(1,i)<0.98*xw || 0.98*xe<xm(1,i) || ym(j,1)<0.98*ys || 0.98*yn<ym(j,1)
            m(j,i)=1;
        end
    end
end

mesh(xm,ym,m);    axis image;  view(0, 90); xlabel('x'); ylabel('y');