% Testroutine zum Programm 
% BezRoot, Blatt 5

% Datenpunkte generieren
P = [1 -2.9 4 -4 2.9 -0.86];

% Aufruf des Programms
R1 = BezRoot(P,0.05);
R2 = BezRoot(P,0.01);

% Ergebnis plotten
t = linspace(0,1);
n = length(P);
figure(1)
clf reset
plot(t,DeCasteljau(P',t))
hold on
plot([0 1],[0,0],'r')
plot(R1(:),0*R1(:),'r+')

figure(2)
clf reset
plot(t,DeCasteljau(P',t))
hold on
plot([0 1],[0,0],'r')
plot(R2(:),0*R2(:),'r+')