% femuddsnurr
  clear, clf
  x=5+[0 5 0 3.5 3 0]; y=[0 2.8 3.3 0 5 0]; 
  plot(0,0,'o', x,y), axis equal, hold on, axis off
  n=10; v=2*pi/n; 
  S=[cos(v) -sin(v); sin(v) cos(v)];
  for k=1:n, P=S*[x;y]; x=P(1,:); y=P(2,:); fill(x,y,'y'), end
