clear all
  close all
 
%
% Input number of random points
%

  N = input('N = ')

%
%  Compute uniformly distributed points and Sobol points.
%

  U = rand(N,2);
  p = sobolset(2,'Skip',1e3,'Leap',1e2);
  S = net(p,N);

%
% Plot the uniformly distributed and Sobol points.
%

  figure(1)
  box on
  axis([0 1 0 1])
  hold on
  scatter(U(:,1),U(:,2),'linewidth',2);
  title('Uniformly Distributed Points')
  set(gca,'Fontsize',[22]);
  hold off

  figure(2)
  box on
  axis([0 1 0 1])
  hold on
  scatter(S(:,1),S(:,2),'linewidth',2);
  title('Sobol Points')
  set(gca,'Fontsize',[22]);
  hold off