function [] = WorstCase(T,n);

%-------------------------------------------
% Bestimmung der Greville Abszisse GA=Mue:
%-------------------------------------------
m=length(T)-n;
for j=1:m
    GA(j)=0;
    for k=1:n-1
        GA(j)=GA(j)+T(k+j);
    end
    GA(j)=GA(j)/(n-1);
end

%--------------------------------------------------------------
% Bestimmung des t Vektors, dieser wird zwischen Mue definiert:
%--------------------------------------------------------------
t=linspace(max(GA(1),T(n)),min(GA(m),T(m+1)),300);

%---------------------
%Bestimmung von sigma²
%---------------------
for j=1:m
    Sigq(j)=0;
    for k=1:(n-1)
        Sigq(j)=Sigq(j)+(T(k+j)-GA(j))^2;
    end
    Sigq(j)=Sigq(j)/(n-2);
end
%-------------------------------
% Bestimmung des P* -Vektors =P:
%-------------------------------
for j=1:m
    P(j)=(GA(j)^2-Sigq(j)/(n-1))/2;
end
%------------------------------------------------
% Neuer Knotenvektor über Mue definieren T'=Tau :
%------------------------------------------------
Tau=GA([1 1:end end]);
%------------------------------------------------
% DeBoor Algorithmus anwenden (B(t)P*)=f(t):
%------------------------------------------------
f=seval(P,T,t);	% = t.^2/2
h=seval(P,Tau,t);
d=f-h;
%-----------------------------------------------------------
% Ergebnis visualisieren:
% ----------------------------------------------------------
figure(1)
clf reset
plot(t,h,'r','LineWidth',2)

hold on
plot(t,f,'b','LineWidth',2)

% FehlerPlot
figure(2)
plot(t,d,'b','LineWidth',2);