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

  <span class="keyword">for</span> j = [1 4 5];
    plot3(R{j}(1,:),R{j}(2,:),R{j}(3,:),RR(j,:),<span class="string">'LineWidth'</span>,2)
  <span class="keyword">end</span>
  axis <span class="string">equal</span>; axis([-2 2 -2 2 0 1.5]);
  xlabel <span class="string">'x [m]'</span>;
  ylabel <span clas