Calculations.ipynb

Kernel: SageMath 7.6

a,b,i0,v0,C = var('a','b','i0','v0','C')

show(solve([i0 == (1-1/(a*b))/((1/b)+C), v0 == a*i0/(1+a*i0)],a,b))

i0 = 0.46
v0 = 0.3
C = 1+0.1*(1+0.02)

sol = solve([i0 == (1-1/(a*b))/((1/b)+C), v0 == a*i0/(1+a*i0)],a,b)

sol

[[a == (150/161), b == (115000/36981)]]

n(sol[0][0].rhs())

0.931677018633540

n(sol[0][1].rhs())

3.10970498364025