
Vt = zeros(500,500);

%Vt(1:50,1:50) = ones(50,50);

rMax = 250;
origin = rMax;

theta2 = 3*pi/4;
theta1 = 3*pi/4 + pi/12;

theta2 = pi/4;
theta1 = pi/4 + pi/12;

theta2 = pi/4 - 2*pi/12;
theta1 = pi/4 - pi/12;


if( 0 )
  for( y = 1:rMax )
    for( x = round(y/tan(theta1)):round(y/tan(theta2)) )
      Vt(y+origin+1,x+origin+1) = 1;
    end
  end
end

if(0)
  theta = 0
  for( r = 1:rMax )
    Vt(origin + round(r*cos(theta)),  origin + round(r*sin(theta))) = 1;
  end
end


theta2 = pi/4 - pi/12;
theta1 = theta2+pi/6;



for( az = theta2:2*pi/2048:theta1 )
  for( r = 1:rMax )
    x = round( r*cos(az) );
    y = round( r*sin(az) );
    Vt(y+origin+1,x+origin+1)=1;
  end
end


imshow(Vt)
set(gca,'ydir','normal','visible','on','fontsize',22)
