CoCalc -- Calculations.ipynb

Views: ¹¹⁴

Kernel: SageMath 7.6

In [1]:

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

In [5]:

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

Out[5]:

In [15]:

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

In [23]:

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

In [24]:

sol

Out[24]:

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

In [25]:

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

Out[25]:

0.931677018633540

In [26]:

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

Out[26]:

3.10970498364025

In [0]: