1 #!/usr/bin/python

     2 import os, sys

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

     4 

     5 

     6 import pylab as pl

     7 from xdync import xdync

     8 

     9 N = pl.array([[0,  4, 0, 1],

    10               [-1, 8, 0, 3]])

    11 dN = 1

    12 

    13 D = pl.array([[0, 0,  1, 4, 0, 1],

    14               [0, 0, -1, 0, 0, 0]])

    15 dD = 2

    16 C = pl.array([[1, 0, 1, 1],

    17               [0, 2, 0, 1]])

    18 dC = 1

    19 

    20 Y, dY, X, dX, B, dB, A, dA = xdync(N, dN, D, dD, C, dC)

    21 

    22 print "Y =", Y

    23 print "dY =", dY

    24 print "X =", X

    25 print "dX =", dX

    26 print "B =", B

    27 print "dB =", dB

    28 print "A =", A

    29 print "dA =", dA

    30 

    31 

    32