def METHOD_NAME():
    g = h.Graph()
    g.size(0, 10, 0, 1)
    g.color(1)
    g.label(0.2, 0.95, "ik vs t (V = 20 mV)")
    g.label(0.8, 1, "")
    s = [h.Section(name="soma_" + str(i)) for i in range(2)]
    s[0].insert("cagk")
    s[1].insert("cagkftab")
    for sec in s:
        sec.diam = 10
        sec.L = 10
    seclamps = [h.SEClamp(i(0.5)) for i in s]
    for se in seclamps:
        se.dur1 = 10
        se.amp1 = 20
        se.dur2 = 1e9
        se.amp2 = -65
    tvec = h.Vector().record(h._ref_t, sec=s[0])
    ikvecs = [h.Vector().record(sec(0.5)._ref_ik, sec=sec) for sec in s]
    def plt(i):
        for j, line in enumerate(ikvecs):
            line.line(g, tvec, [8, i][j], [4, 1][j])
        assert ikvecs[0].c().sub(ikvecs[1]).abs().max() < 1e-4
    def runvc(i, cai):
        for sec in s:
            sec.cai = cai
        h.finitialize(-65)
        while h.t < 15.0:
            h.fadvance()
        plt(i)
    for i, cai in enumerate([1e-3, 1e-2, 1e-1]):
        g.color(i + 1)
        g.label("cai=%g" % cai)
        runvc(i + 1, cai)
    return g, tvec, ikvecs