theta = [0:0.08:2.5 * pi];

theta2 = [0:0.08:2.5 * pi];

%rand( 70, 1) * 2 * pi;

S2 = 0.25 * (rand( 30, 2 ) -0.5 );


S = 1 * [ cos(theta'), sin(theta') ];
S2 =  [ cos(theta2'), sin(theta2') ];

S2 = S2 .* ( 1 - 1.1 .* (1/(6*pi)) .* ( [theta2', theta2'] ) );
S = 0.8 .* S .* ( 1 - 1.2 .* (1/(5*pi)) .* ( [theta', theta'] ) );


S3 = [ S; S2 ];


[clusterLabels] = ngcluster( S3, 2, 0.044 );
