r to scale the axes</p><pre class="language-matlab">colormap <span class="string">gray</span>
D = M(round(1:info.dime.dim(2)/2),:,:,1); <span class="comment">% split along sagittal plane</span>
Ds = smooth3(D);
hiso = patch(isosurface(Ds,0),<span class="keyword">...</span>
  <span class="string">'FaceColor'</span>,[.5,.5,.8],<span class="keyword">...</span>
  <span class="string">'EdgeColor'</span>,<span class="string">'none'</span>);
isonormals(Ds,hiso)
hcap = patch(isocaps(D,0),<span class="keyword">...</