equal
deleted
inserted
replaced
|
1 #!/usr/bin/env python |
|
2 # 11.5 |
|
3 import os, sys |
|
4 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"] |
|
5 |
|
6 import pylab as pl |
|
7 from xdync import xdync |
|
8 from polyfuncs import polmul |
|
9 |
|
10 def mv(A, dA, B, dB, C, dC, k, int1): |
|
11 zk = pl.zeros(k+1) |
|
12 zk[-1] = 1 |
|
13 if int1 >= 1: |
|
14 A, dA = polmul([1, -1], 1, A, dA) |
|
15 Fj, dFj, Ej, dEj = xdync(zk, k, A, dA, C, dC)[:4] |
|
16 Gk, dGk = polmul(Ek, dEk, B, dB) |
|
17 S, dS = Fk, dFk |
|
18 R, dR = Gk, dGk |
|
19 return S, dS, R, dR |
|
20 |