#!/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