casteljau2.m

function casteljau2(t,b,a)
y=0:0.01:1;
B0=(1-y).^2;
B1=2*(1-y).*y;
B2=y.^2;
B=[B0;B1;B2];
%pol bernstein in camp de vect
f=b*B; %curba bezier
%repr

b1=zeros(2,2); b2=zeros(2,1);
for i=1:2
b1(:,i)=b(:,i)*(1-t)+b(:,i+1)*t;
end
b1;
b2(:,1)=b1(:,1)*(1-t)+b1(:,2)*t;
if a==1
hold on


plot(f(1,:),f(2,:),'-k','LineWidth',2)
plot(b(1,:),b(2,:),'-ko','MarkerSize',5)
plot(b(1,1),b(2,1),'-ko','MarkerSize',5,'MarkerFaceColor','k')
plot(b(1,end),b(2,end),'-ko','MarkerSize',5,'MarkerFaceColor','k')
% plot(b1(1,:),b1(2,:),'g*')
plot(b1(1,:),b1(2,:),'k-')
% plot(b2(1,:),b2(2,:),'k*')
plot(b2(1,:),b2(2,:),'k-')
% plot(b3(1,:),b3(2,:),'o','MarkerSize',2,'MarkerEdgeColor','k','MarkerFaceColor','r')
hold off 
elseif a==0
    plot(f(1,:),f(2,:),'-k','LineWidth',2)
else
eror('parametrul trebuie sa fie 1 (cu poligon de ontrol) sau 0 (fara poligon de control)')
end