function res = h3d_test03(axis)

angular_cell_size_deg = 8; %7;
rho_cell_size = 0.3;
rho_min = -1.5;
rho_max = 1.5;

base = h3d_create_struct(angular_cell_size_deg, rho_cell_size, rho_min, rho_max);

points = h3d_test_shape('pyramid');

hts = {};

delta = 9;
angles = 0:delta:360-delta

for i=1:numel(angles)
	angle = angles(i);
	fprintf('angle: %d\n', angle);

	R1 = rot_axis_angle(axis, deg2rad(angle) );
	
	s1 = h3d_compute_cube_ht(base, R1*points);

	hts{i} = s1;
end

res.hts = hts;
res.axis = axis;