gistfile1.py

#!/usr/bin/env python

import matplotlib.pyplot as plt
from numpy import vectorize, sqrt, linspace

R = 1
C = 0.2
w = 1


@vectorize
def impedance(l):
    return sqrt(R**2 + (w*l - 1.0/(w*C))**2)



L = linspace(0.0001, 40.0, 200)

I = impedance(L)

plt.plot(L, I)
plt.xlabel('$L$')
plt.ylabel('$Z_\mathrm{eff}$')
plt.show()