

  clear all
%close all

  n = 101;
  p = 1e+5;
  ind = 1e+3;
  ell = 75;
  h = 1/(n-1);
  t = 0:h:1;


  [I,T] = meshgrid(1:p,t);
  A = sin(2*pi*I.*T);


% RandSVD on A (using an oversampling of 10)

  k = 65; 
  tic
  [U1,S1,V1] = randsvd_case_iii(A,k,10,ind);
  toc

  figure(1)
  semilogy(1:k,diag(S1), 'ob', 'linewidth',3)
  axis([0 k 1e-7 1e+5])
  set(gca,'Fontsize',[22]);
  xlabel('Index')
  ylabel('Singular Values')

