プログラムβ版

%初期値設定
n=3;
Px=[1 2 4 3];
Py=[1 3 3 1];

%制御点のプロット
plot(Px, Py);
hold on;

%ベジェ曲線の計算とプロット
for t=[0 0.15 0.35 0.5 0.65 0.85 1.0];
  rx=0;
  ry=0;
  for i=0:n
    rx=rx+factorial(n)/factorial(n-i)/factorial(i)*((1-t)^(n-i))*(t^(i))*Px(i+1);
    ry=ry+factorial(n)/factorial(n-i)/factorial(i)*((1-t)^(n-i))*(t^(i))*Py(i+1);
  end
  plot(rx, ry,'r');
end

時間無くて数値入力タイプのはうまくできなかったので、
とりあえずこれをＵＰ