function C = DeCasteljau(P,T);

[n,d] = size(P);
m = length(T);

for k=1:m
  p = P;
  for i=2:n
    p = (1-T(k)).*p(1:end-1,:) + T(k).*p(2:end,:);
  end
  C(k,:) = p;
end