% bez3fler, program till figur 4 i bezierhäftet
% Styckvis kubiska bezierkurvor mellan fem punkter 
  clear, clf
  dt=0.05;  t=(dt:dt:1)';                 
  F= [(1-t).^3  3*t.*(1-t).^2  3*t.^2.*(1-t)  t.^3];
  n=5; 
  p=[0 5; 8 10; 15 20; 24 10; 16 4];      % givna punkter
  plot(p(:,1),p(:,2),'o', p(:,1),p(:,2),':')
  axis equal, hold on 
  k=[   1         0           % givna normerade riktningar 
    cos(pi/3)   sin(pi/3) 
        1         0 
        0        -1 
    cos(3*pi/4) sin(3*pi/4)];         
  faktor=0.4;                             % lagom buktning
  rr=p(1,:);                              % första punkten
  for i=1:n-1 
    alfa=faktor*norm(p(i+1,:)-p(i,:)); 
    gamma=alfa;   
    b=p(i,:)+alfa.*k(i,:);        
    c=p(i+1,:)-gamma.*k(i+1,:);    
    r=F*[p(i,:); b; c; p(i+1,:)];         % i:te kurvan
    rr=[rr; r];                           % bygg på totalkurvan
  end
  plot(rr(:,1),rr(:,2))

  

