ejercicio1

Contents

apartado c

% He reescrito las funciones factr y combina en este fichero para evitar
% problemas con el directorio. Pero sé que podría usar los ficheros que las
% contienen de los apartados a y b.

t=linspace(0,1);
n=3;
for i=0:n
    b=bernstein(n,i,t);
    plot(t,b);
    hold on;
end
xlabel('t')
ylabel('Polinomio de Bernstein');
title('Polinomios de Bernstein de grado 3');
legend('B_3_,_0','B_3_,_1','B_3_,_2','B_3_,_3');

function b = bernstein(n,i,t)
b=combina(n,i).*t.^i.*(1-t).^(n-i);
end

function f = combina(n,i)
f = factr(n) / (factr(i) * factr(n-i));
end

function f =  factr (n)
if n == 0;
   f = 1;
else;
   f = n *  factr (n-1);
end
end

Published with MATLAB® R2017a