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