function [xx,yy]=fbez(x,y);
% Styckvisa kubiska bezierkurvor, riktningar med centraldiffkvot 
  t=(0.05:0.05:1)'; 
  F=[(1-t).^3   3*t.*(1-t).^2   3*t.^2.*(1-t)  t.^3];          
  n=length(x);  p=[x y];  dp=diff(p);  
  h=sqrt(dp(:,1).^2+dp(:,2).^2);  
  k=[];
  for i=2:n-1, ki=p(i+1,:)-p(i-1,:); k=[k; ki/norm(ki)]; end
  k1=2*dp(1,:)/h(1)-k(1,:);       k1=k1/norm(k1);  k=[k1; k];  
  kn=2*dp(n-1,:)/h(n-1)-k(n-1,:); kn=kn/norm(kn);  k=[k; kn];
  alfa=0.4*[h h];   gamma=alfa;               % 0.4 är lagom buktfaktor  
  b=p(1:n-1,:)+alfa.*k(1:n-1,:);  
  c=p(2:n,:)-gamma.*k(2:n,:);
  r=p(1,:);
  for i=1:n-1, xy=F*[p(i,:); b(i,:); c(i,:); p(i+1,:)];  r=[r; xy]; end
  xx=r(:,1);  yy=r(:,2);
 

  
