0
|
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 |
|