minv/mv.py
changeset 0 0efde00f9229
equal deleted inserted replaced
-1:000000000000 0:0efde00f9229
       
     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