function risi_bez(b,u,nt,ink);
  k=length(u)-1;
  hold on;
  format=ink;
  plot(b(1,:),b(2,:),'k');
  plot(b(1,1:3:3*k+1,:),b(2,1:3:3*k+1),'o');
  for i=1:k;
    ui=linspace(u(i),u(i+1),nt);
    ti=(ui-u(i))/(u(i+1)-u(i));
    %Bi=b(:,3*i-2:3*i+1)*[nchoosek(3,0)*(1-ti).^3;
    %	nchoosek(3,1)*ti.*(1-ti).^2;
    %	nchoosek(3,2)*ti.^2.*(1-ti);
    %	nchoosek(3,0)*ti.^3];
    Bi=bezier(b(:,3*i-2:3*i+1),ti);
    plot(Bi(1,:),Bi(2,:),format)
  end;
  axis equal;
  hold off
