b1I1 / b2I2

 

[Previously saved workspace restored]

> b1 = 1

> b1 = 2

> b2= 1.8

> I1 = c(0.005, 0.5, 0.9)

> I2 = c(0.8, 0.5, 0.1)

> b1*I1 / b2 * I2

[1] 0.004444444 0.277777778 0.100000000

> ratio = b1*I1 / b2 * I2

> t = 1:3

> plot( ratio ~ t)

> t = 1:5

> t

[1] 1 2 3 4 5

> I1 = 0.005*2^5

> I1

[1] 0.16

> I1 = 0.005*2^t

> I1

[1] 0.01 0.02 0.04 0.08 0.16

> I2 = 1 - I1

> ratio = b1*I1 / b2 * I2

> plot( ratio ~ t)

> plot( log(ratio) ~ t)