M=[1-10i 1; 1 1-10i];


R = 1;
I = 1000;
%T(I/2)=1+1i;
%for j = 1:I/2;
%T(I/2+j) = MobiusOnPoint(M,T(I/2+j-1));
%end

%%%T(I)=1+1i;
for k = 1:2

if k > 1;
		M=inv(M);
end
 
T = 1+i;
for j = 2:I
T(j) = MobiusOnPoint((M),T(j-1));
end


D = [real(T); imag(T)];
[u,s,v] = svd(D,0);

D = D + repmat(u(:,1),1,size(D,2));

[p,N] = StereographicProj(D,0);

P{k} = p;

if k > 1; 
		P{k} = p(:,end:-1:1);
end

end

% Q is a set of points lying in a spiral

Q=[P{2} P{1}];
plot3(Q(1,:),Q(2,:),Q(3,:),'-o','LineWidth',1);

