function [  ] = zeichnegeschlossenekurve( x,y )
% Zeichnung einer geschlossenen Kurve durch die Datenpunkte (x(i),y(i))
% Erster Punkt muss mit letztem übereinstimmen
%
n=length(x);
%
par=3
switch par
    case 1
        % Parametrisierung 1
        for i = 1:n
            tau(i)=i;
        end;
%
    case 2
        % Parametrisierung 2
        tau(1)=0;
        for i = 2:n
            tau(i)=tau(i-1)+sqrt((x(i)-x(i-1))*(x(i)-x(i-1))+(y(i)-y(i-1))*(y(i)-y(i-1)));
        end;
    case 3
        % Parametrisierung 3
        tau(1)=0;
        for i = 2:n
            tau(i)=tau(i-1)+((x(i)-x(i-1))*(x(i)-x(i-1))+(y(i)-y(i-1))*(y(i)-y(i-1)));
        end;
    otherwise
        fprintf('Parametrisierung %4i nichtdefiniert in <zeichnegeschliossenekurve>\n',par)
end;
%
cx=periodischerspline(tau,x);
cy=periodischerspline(tau,y);
%
t=linspace(0,1,10);
psi1=1-t;
psi2=t;
psi3=-1/6*t.*t.*t+1/2*t.*t-1/3*t;
psi4=1/6*t.*t.*t-1/6.*t;
%
clf;
figure(1);
axis equal;
hold on;
plot(x,y,'o');
%
for i = 1:n-1
    h=tau(i+1)-tau(i);
    p=x(i)*psi1+x(i+1)*psi2+cx(i)*h*h*psi3+cx(i+1)*h*h*psi4;
    q=y(i)*psi1+y(i+1)*psi2+cy(i)*h*h*psi3+cy(i+1)*h*h*psi4;
    plot(p,q);
end;
end

