function gps_order_growth_plot ( )

  order = zeros(33,1);

  for level = 0 : 32
    if ( level == 0 )
      order(level+1) = 1;
      prec(level+1)=1;
    else
      p = 5;
      o = 3;
      while ( p < 2 * level + 1 )
        p = 2 * p + 1;
        o = 2 * o + 1 ;        
      end
      order(level+1) = o;
      prec(level+1) = p;
    end
  end

  bar ( 0:32, prec, 'r' )
  hold on
  bar ( 0:32, order, 'b' )
  plot ( [0,32], [0,65], 'w' )
  xlabel ( 'Level' )
  ylabel ( 'GPS Order and Accuracy' )
  title ( 'GPS Order (blue) and Accuracy (red) as Function of Level' )
  hold off

  return
end

