(:,1:M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),<span class="string">'r'</span>,<span class="string">'EdgeColor'</span>,<span class="string">'none'</span>);
  P = [r R{1}(:,M/2+1:3*M/4+1) r]; fill3(P(1,:),P(2,:),P(3,:),<span class="string">'r'</span>,<span class="string">'EdgeColor'</span>,<span class="string">'none'</span>);

  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,:),<span class="string">'b'</span>,<span class="string">'EdgeColor'</span>,<sp