% bez3h,  kubisk bezierkurva mellan p1 och p2   (sju varianter)
% En styrpunkt fix, den andra flyttas men riktningen är fix
  clear, clf
  t=(0:0.05:1)';  
  F=[(1-t).^3   3*t.*(1-t).^2  3*t.^2.*(1-t)  t.^3]; 
  p1=[0 5];  p2=[18 0];  p=[p1; p2]
  c=[20 10];                          % Styrpunkten c fix
  plot(p(:,1),p(:,2),':', p(:,1),p(:,2),'o')  
  axis([0 40 0 35]), axis('equal'), hold on
  title('Kubiska bezierkurvor')
  k1=[1 1];  k1=k1/norm(k1);          % Startriktning k1 
  alfa=0;     
  for kurvnr=1:7        
    alfa=alfa+10;         % Variera styrpunktsavståndet alfa
    b=p1+alfa*k1;  
    pbcp=[p1; b; c; p2];  r=F*pbcp;  
    plot(pbcp(:,1),pbcp(:,2),':', r(:,1),r(:,2))
    drawnow
  end

  
