1 #!/usr/bin/python

     2 # 7.9 

     3 import os, sys

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

     5 

     6 import pylab as pl

     7 from left_prm import left_prm

     8 

     9 D = pl.array([[1, 0, 0, 0, 0, 0],

    10               [0, 1, 0, 1, 0, 0],

    11               [0, 0, 1, 1, 1, 0]])

    12 

    13 N = pl.array([[1, 0, 0],

    14               [0, 1, 0],

    15               [0, 0, 1]])

    16 

    17 dD, dN = 1, 0

    18 

    19 B, dB, A, dA = left_prm(N, dN, D, dD)

    20 

    21 print "A =", A

    22 print "dA =", dA

    23 

    24 print "B =", B

    25 print "dB =", dB

    26