close all;
clear all;



% Definition der Konstanten
dim    = 2;       % Dimension der Datenpunkte
m      = 5000;   % Anzahl der Datenpunkte
phi    = 135;     % Drehwinkel in Grad
length = [1 100];  % Anzahl der Komponenten = dim !
mean   = [0  0];  % Anzahl der Komponenten = dim ! 


% Erzeugung der normalverteilten Datenpunkte
samples = randn(m, dim);


% Drehung und Streckung der Daten 
trans   =  diag(sqrt(length));

if dim == 2
    alpha   = phi/180 * pi;
    r       = [ cos(alpha)  sin(alpha) ; 
               -sin(alpha)  cos(alpha) ];
    trans   = trans * r;
end;

samples =  samples * trans;


samples =  samples +  ones(m,1) * mean;


% Ausgabe
if dim == 2
  plot(samples(:,1),samples(:,2) ,'.');
  title('data');
  print samples.ps -dps;
end;

save samples.dat samples -ascii;
print wolke.eps -deps


