function part2sim3cL(L,sigma)
  truecapacitiesL = [];
  approxcapacitiesL = [];
  
  Lvector = 10:10:(2*L);
  for (Ltry=Lvector)
    Ntry = Ltry
    [truecapacity,approx] = part2trueapproxcapacity2(Ltry,Ntry,sigma,200);
    truecapacitiesL = [ truecapacitiesL truecapacity ];
    approxcapacitiesL = [ approxcapacitiesL approx ];
  end
  
  plot(Lvector,truecapacitiesL,'k+',Lvector,approxcapacitiesL,'ko');
  xlabel('L');
  legend('True capacitiy','Estimated capacity');