MT.py 文件源码

def run(n, plotIt=True):
    # something to make a plot

    F = frange(-5., 5., 20)
    H = thick(50., 100., n)
    sign = sig(-5., 0., n)
    mun = mu(1., 2., n)
    epsn = eps(1., 9., n)
    chg = np.zeros_like(sign)
    taux = np.zeros_like(sign)
    c = np.zeros_like(sign)

    Res, Phase = appres(F, H, sign, chg, taux, c, mun, epsn, n)

    if plotIt:

        PlotAppRes(F, H, sign, chg, taux, c, mun, epsn, n, fenvelope=1000., PlotEnvelope=True)

    return Res, Phase