import math

RTT1 = 0.01
RTT2 =0.0442857142857
RTT3 =0.0785714285714
RTT4 =0.112857142857
RTT5 =0.147142857143
RTT6 =0.181428571429
RTT7 =0.215714285714
RTT8 =0.25

alpha1 = 1
alpha2 = 1
alpha3 = 1
alpha4 = 1
alpha5 = 1
alpha6 = 1
alpha7 = 1
alpha8 = 1

beta1 = 0.50/RTT1*RTT1
beta2 = 0.50/RTT2*RTT2
beta3 = 0.50/RTT3*RTT3
beta4 = 0.50/RTT4*RTT4
beta5 = 0.50/RTT5*RTT5
beta6 = 0.50/RTT6*RTT6
beta7 = 0.50/RTT7*RTT7
beta8 = 0.50/RTT8*RTT8

drop_prob = 0.589985218787

W1=W2=W3=W4=W5=W6=W7=W8=0
simtime =0
delta = 0.01
drift1 = drift2 = drift3 = drift4 = drift5 = drift6 = drift7= drift8= 1
totWinSize = 0


while (1):

  simtime+=delta

  drift1 = alpha1/RTT1 - ( (beta1* drop_prob *  W1*W1) / (RTT1)  )
  W1 = W1 + delta * drift1

  drift2 = alpha2/RTT2 - ( (beta2* drop_prob *  W2*W2) / (RTT2)  )
  W2 = W2 + delta * drift2

  drift3 = alpha3/RTT3 - ( (beta3* drop_prob *  W3*W3) / (RTT3)  )
  W3 = W3 + delta * drift3

  drift4 = alpha4/RTT4 - ( (beta4* drop_prob *  W4*W4) / (RTT4)  )
  W4 = W4 + delta * drift4

  drift5 = alpha5/RTT5 - ( (beta5* drop_prob *  W5*W5) / (RTT5)  )
  W5 = W5 + delta * drift5

  drift6 = alpha6/RTT6 - ( (beta6* drop_prob *  W6*W6) / (RTT6)  )
  W6 = W6 + delta * drift6

  drift7 = alpha7/RTT7 - ( (beta7* drop_prob *  W7*W7) / (RTT7)  )
  W7 = W7 + delta * drift7

  drift8 = alpha8/RTT8 - ( (beta8* drop_prob *  W8*W8) / (RTT8)  )
  W8 = W8 + delta * drift8

  totWinSize = (4*W1) + (4*W2) + (4*W3) + (4*W4) + (4*W5) +(4*W6) + (4*W7) +(4*W8)
  print totWinSize ,simtime