slut=4*pi;antal=1000;steg=slut/antal;
t=0:steg:slut-steg;
A=[0 1;-1 0];
y(:,1)=[0,1]';
for k=1:(antal-1)
  y(:,k+1)=y(:,k)+steg*A*y(:,k);
end
