function n = hex_grid_angle_01_size ( center, angle, h )

%% HEX_GRID_ANGLE_01_SIZE counts the points in an angled hex grid in the unit box.
%
%  Modified:
%
%    24 October 2006
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, real CENTER(2), the center of the grid.
%    This point must be inside the unit square.
%
%    Input, real ANGLE, the angle, in degrees, of the grid.
%    Normally, 0 <= ANGLE <= 180, but any value is allowed.
%
%    Input, real H, the spacing between neighboring
%    points on a grid line.
%
%    Output, integer N, the number of points of the angled hex grid
%    that are within the unit square.
%
  dim_num = 2;
%
%  Ninny checks.
%
  if ( ~box_01_contains_point_2d ( center ) )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'HEX_GRID_ANGLE_01_SIZE - Fatal error!\n' );
    fprintf ( 1, '  The center point of the grid is not\n' );
    fprintf ( 1, '  inside the unit square.\n' );
    fprintf ( 1, '  CENTER = (%f, %f)\n', center(1:dim_num) );
    error ( 'HEX_GRID_ANGLE_01_SIZE - Fatal error!' );
  end

  if ( h == 0.0 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'HEX_GRID_ANGLE_01_SIZE - Fatal error!\n' );
    fprintf ( 1, '  The grid spacing must be nonzero.\n' );
    fprintf ( 1, '  H = %f\n', h );
    error ( 'HEX_GRID_ANGLE_01_SIZE - Fatal error!' );
  end

  n = 0;

  layer = 0;
  point(1:dim_num) = center(1:dim_num);

  n = 1;

  while ( 1 )

    layer = layer + 1;

    layer_size = 0;

    angle2 = angle;
%
%  Compute the first point on the new layer.
%
%  Note that some MATLAB ninnyosity prevents me from
%   point(1:dim_num) = point(1:dim_num) ...
%     + h * [ cos_deg ( angle2 ), sin_deg ( angle2 ) ];
%  
    point(1) = point(1) + h * cos_deg ( angle2 );
    point(2) = point(2) + h * sin_deg ( angle2 );

    angle2 = r8_modp ( angle2 + 60.0, 360.0 );

    for i = 1 : 6

      angle2 = r8_modp ( angle2 + 60.0, 360.0 );

      for j = 1 : layer

        point(1) = point(1) + h * cos_deg ( angle2 );
        point(2) = point(2) + h * sin_deg ( angle2 );

        if ( box_01_contains_point_2d ( point ) )
          layer_size = layer_size + 1;
          n = n + 1;
        end

      end

    end

    if ( layer_size == 0 )
      break
    end

  end