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