function out=var_data(a,N,x0)
%x0 ... Anfangsbedingung in Spaltenvektor
%erzeugt Daten durch VAR-Model 
%Länge des Datensatzes n
[n1,n2]=size(a);
noise=randn(n2,N);
data=zeros(n2,N);
x=x0;
for j=1:N
    x=a*x+noise(:,j);
    data(:,j)=x;
end
out=data';
