function part2sim3bsigma(L,sigma)
  truecapacitiessigma = [];
  approxcapacitiessigma1 = [];
  approxcapacitiessigma2 = [];
  
  sigmavector=sigma:0.02:(sigma+2);
  for (sigmatry=sigmavector) 
    sigmatry
    N = L;
    [truecapacity,approx1,approx2] = part2trueapproxcapacity1(L,N,sigmatry);
    truecapacitiessigma = [ truecapacitiessigma truecapacity ];
    approxcapacitiessigma1 = [ approxcapacitiessigma1 approx1 ];
    approxcapacitiessigma2 = [ approxcapacitiessigma2 approx2 ];
  end
  
  plot(1./(sigmavector.^2),truecapacitiessigma,'k-',1./(sigmavector.^2),approxcapacitiessigma1,'k+',1./(sigmavector.^2),approxcapacitiessigma2,'ko');
  xlabel('\rho')
  legend('True capacitiy','Estimated capacity, 3 terms','Estimated capacity, 7 terms');