% bsplejip, B-splines med ginput
  clear, clf
  t=(0:0.1:1)';
  F=[(1-t).^3  3*t.^3-6*t.^2+4  -3*t.^3+3*t.^2+3*t+1  t.^3]/6;
  axis([0 100 0 80]), hold on
  disp(['Klicka punkter (bryt med musens mittknapp)'])
  [x,y]=ginput(1); plot(x,y,'o'), p1=[x y];  
  [x,y]=ginput(1); plotus(x,y,'o'), plotus([p1(1) x],[p1(2) y],'m:')
  p=[x y];  P=[p1; p1; p1; p];      % 2 kopior av första punkten behövs 
  r=F*P; plotus(r(:,1),r(:,2),'r-')
  for i=1:48                        % maxantal punkter är 50
    [x,y,button]=ginput(1);  
    if button==2, break
    else plotus(x,y,'o'), plotus([p(1) x],[p(2) y],'m:')
      p=[x y];  P=[P(2:4,:); p];    % utnyttja de fyra senaste punkterna
      r=F*P; plotus(r(:,1),r(:,2),'r-')
    end
  end
  for i=1:2                         % avsluta med 2 kopior av slutpunkten
    P=[P(2:4,:); p]; r=F*P; plotus(r(:,1),r(:,2),'r-')
  end






