Initial commit.
#!/usr/bin/python
import os, sys
sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
import pylab as pl
from xdync import xdync
N = pl.array([[0, 4, 0, 1],
[-1, 8, 0, 3]])
dN = 1
D = pl.array([[0, 0, 1, 4, 0, 1],
[0, 0, -1, 0, 0, 0]])
dD = 2
C = pl.array([[1, 0, 1, 1],
[0, 2, 0, 1]])
dC = 1
Y, dY, X, dX, B, dB, A, dA = xdync(N, dN, D, dD, C, dC)
print "Y =", Y
print "dY =", dY
print "X =", X
print "dX =", dX
print "B =", B
print "dB =", dB
print "A =", A
print "dA =", dA