%?the?script?to?compute?the?system?throughput?following?Bianchi‘s?analysis
%?G.?Bianchi?etc.?“Remarks?on?IEEE?802.11?DCF?Performance?Analysis“?IEEE
%?Comm.?Letters?Aug.?2005.

for?k?=?1:100
N?=?k;??%?number?of?nodes

payload?=?1500;?%?in?unit?of?bytes?assuming?fixed-size?payload?here

R?=?7;??%?retransmission?limit

%?802.11b
DIFS?=?50;??%?us
SIFS?=?10;??%?us
slot?=?20;??%?us
dataRate?=?11;?%?Mbps
CWmin?=?31;?%?we?assume?the?range?is?[0?CWmin]?so?it?is?31?rather?than?32?
ACK?=?248;?%?us
MAC?=?28;?%?bytes
PHY?=?192;?%?us

precision?=?10000;
%?compute?p?and?tau
for?i?=?1:precision+1
????p1?=?（i-1）*1.0/precision;
????
????%?tau?in?terms?of?p
????tau_down?=?0;
????for?j?=?1:R+1
????????tau_down?=?tau_down?+?p1^（j-1）*（（CWmin+1）*2^（j-1）?-?1）?/?2;
????end
????tau_down?=?1?+?tau_down*（1-p1）/（1-p1^（R+1））;
????tau