function part2sim3bL(L,sigma)
  truecapacitiesL = [];
  approxcapacitiesL1 = [];
  approxcapacitiesL2 = [];
  
  Lvector = 5:2:(2*L);
  for (Ltry=Lvector)
    Ntry = Ltry
    [truecapacity,approx1,approx2] = part2trueapproxcapacity1(Ltry,Ntry,sigma);
    truecapacitiesL = [ truecapacitiesL truecapacity ];
    approxcapacitiesL1 = [ approxcapacitiesL1 approx1 ];
    approxcapacitiesL2 = [ approxcapacitiesL2 approx2 ];
  end
  
  approxcapacitiesL1 = abs(approxcapacitiesL1-truecapacitiesL);
  approxcapacitiesL2 = abs(approxcapacitiesL2-truecapacitiesL);
  plot(Lvector,approxcapacitiesL1,'k+',Lvector,approxcapacitiesL2,'ko');
  xlabel('L');
  legend('Deviation from true capacity for 3 terms estimated capacity','Deviation from true capacity for 7 terms estimated capacity');