1 #!/usr/bin/env python

     2 # 11.3

     3 import os, sys

     4 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

     5 

     6 import pylab as pl

     7 from xdync import xdync

     8 

     9 C, dC = pl.array([1, 0.5]), 1

    10 A, dA = pl.array([1, -0.6, -0.16]), 2

    11 j = 2

    12 

    13 zj = pl.zeros(j+1)

    14 zj[-1] = 1

    15 Fj, dFj, Ej, dEj = xdync(zj, j, A, dA, C, dC)[:4]

    16 print Fj, dFj, Ej, dEj 

    17 print "Fj, dFj, Ej, dEj "

    18