function sexagon(a0,b0,radius)
%
%
% displaying a circle.
%
x_first=a0+radius;
y_first=b0;
n_sides=6;
angle_read=pi/180*(360/n_sides);
for j=1:n_sides
  x_next=a0+radius*cos(j*angle_read);
  y_next=b0+radius*sin(j*angle_read);
  line([x_first x_next], [y_first y_next]);
  x_first=x_next;
  y_first=y_next;
end % j