% cosbez,  cos(x) approx. av kvadratisk och kubisk bezierkurva
  clear, clf
  t=(0:0.05:1)';  
  F2=[(1-t).^2  2*t.*(1-t)  t.^2];                   % kvadr. bezier
  F3=[(1-t).^3  3*t.*(1-t).^2  3*(1-t).*t.^2  t.^3]; % kubisk bezier
  v=-pi/2:pi/100:pi/2;  
  plot([-pi/2 pi/2],[0 0],'o', v,cos(v),'y:'),
  text(-1.5,1,'Cosinuskurvan prickad') 
  axis([-1.6 1.6 0 1.1]), hold on, drawnow
  p1=[-pi/2 0];  p2=[pi/2 0]; 
  b=[0 2];  pbp=[p1; b; p2];                % styrpunkt kvadr.bez.
  r=F2*pbp;  plot(r(:,1),r(:,2),'-.')
  b=[-pi/2+4/3  4/3];  c=[pi/2-4/3  4/3];   % styrpunkter kub.bez.
  pbcp=[p1; b; c; p2];
  r=F3*pbcp;  plot(r(:,1),r(:,2))
  title('Kvadratisk bezierkurva streckprickad')
  xlabel('Kubisk bezierkurva heldragen')


plot(-pi/2*(1-2*t),4*t.*(1-t),'g')
