minv/mv.py
changeset 0 0efde00f9229
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/minv/mv.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,20 @@
+#!/usr/bin/env python
+# 11.5
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from xdync import xdync
+from polyfuncs import polmul
+
+def mv(A, dA, B, dB, C, dC, k, int1):
+    zk = pl.zeros(k+1)
+    zk[-1] = 1
+    if int1 >= 1:
+        A, dA = polmul([1, -1], 1, A, dA)
+    Fj, dFj, Ej, dEj = xdync(zk, k, A, dA, C, dC)[:4]
+    Gk, dGk = polmul(Ek, dEk, B, dB)
+    S, dS = Fk, dFk 
+    R, dR = Gk, dGk
+    return S, dS, R, dR
+