i,201); plot3(2*sin(aa),2*cos(aa),0*aa,'k:','LineWidth',2);
  hold on
  
  r = A*[0; 0; rw] + [x; y; 0];
  P = [r R{1}(:,1:M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),'r','EdgeColor','none');
  P = [r R{1}(:,M/2+1:3*M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),'r','EdgeColor','none');
  
  r = A*[rr*sin(psif); 0; rw+rr*cos(psif) ] + [x; y; 0];
  P = [r R{5}(:,1:M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),'b','EdgeColor','none');
  P = [r R{5}(:,M/2+1:3*M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),'b','EdgeColor','none');
  
  for j 