1from lcapy import Circuit
2cct = Circuit()
3
4cct.add('V 1 0 step 10')
5cct.add('L 1 2 1e-3')
6cct.add('C 2 3 1e-4')
7cct.add('R 3 0 1')
8
9import numpy as np
10t = np.linspace(0, 0.01, 1000)
11vr = cct.R.v.evaluate(t)
12
13from matplotlib.pyplot import figure, savefig
14fig = figure()
15ax = fig.add_subplot(111)
16ax.plot(t, vr, linewidth=2)
17ax.set_xlabel('Time (s)')
18ax.set_ylabel('Resistor voltage (V)')
19ax.grid(True)
20
21savefig('circuit-VRLC2-vr.png')
22