~dinko-metalac/calculus-app2/trunk

        r2 = jtheta(1, z, q)

    assert r1.ae(r2)

    mp.dps = 15

    assert jtheta(3, 4.5, 0.25, 9).ae(1359.04892680683)

    assert jtheta(3, 4.5, 0.25, 50).ae(-6.14832772630905e+33)

    mp.dps = 50

    res = mpf('0.9762691700944007312693721148331') - \

          mpf('0.0307203994181623243583169154824')*j

    r = jcn(u, m)

    assert(mpc_ae(r, res))

 

    #N[JacobiDN[1/4 + I/8, 1/3 + I/7], 35]

        assert elliprg(0, mpf('0.0796'), 4).ae('1.0284758090288040009838871385180217366569777284430590125081211090574701293154645750017813190805144572673802094')

    mp.dps = 15

 

    assert isnan(qfrom(m=nan))