function y = Bern (i,n,t)
%
% Bern(i,n,t) = B_{i}^n(t)= nchoosek(n,i-1)*t.^(i-1).*(1-t).^(n-i+1);
%
% Polinomio di Bernstein di grado n
%
if ((i == 0) || (i > n + 1))
  y = zeros (size (t));
elseif ((n == 0) && (i == 1))
  y = ones (size (t));
else
  y = (1 - t) .* Bern (i, n - 1, t) + t .* Bern (i - 1, n - 1, t);
end
%!demo
%! t = linspace(0,1,101);
%! i = 2;
%! n = 2;
%! plot(t,Bern(i,n,t),'*',t,nchoosek(n,i-1)*t.^(i-1).*(1-t).^(n-i+1))
%! legend('Bern(2,2,t) per ricorrenza','Bern(2,2,t) da definizione')
