
Sub Pulsante1_Click() Ndays = 100 Ntrials = 100 Randomize Pi = 3.14 For T = 1 To Ndays For j = 1 To Ntrials u1 = Rnd u2 = Rnd C = Cells(10 + T, 1).Value N = 0.23 + 0.15 * (-Log(u1)) * Cos(2 * Pi * u2) k = 0.95 g = 123 kRT = Exp(323 * ((1 / 15) - (1 / 22))) kcl = 0.2 r0 = 123 ro = r0 * k * g * kcl * (365 / 28) ^ N m0 = 875 k = 2.63 Ccr = 0.04 P = 1 + k * (C - Ccr) V = m0 * P / ro Vd = V * 0.32 Next j Next T For T = 1 To Ndays Cells(T + 10, 2).Value = Vd Next T End Sub |
