1 files changed, 34 insertions, 0 deletions
diff --git a/2048/DEPENDENCIES/polmul.sci b/2048/DEPENDENCIES/polmul.sci
new file mode 100755
index 000000000..67e4adb9a
--- /dev/null
+++ b/2048/DEPENDENCIES/polmul.sci
@@ -0,0 +1,34 @@
+// polmul
+// The command
+//    [C,degA] = polmul(A,degA,B,degB)
+// produces the polynomial matrix C that equals the product A*B of the
+// polynomial matrices A and B.
+//
+// H. Kwakernaak, July, 1990
+
+
+function [C,degC] = polmul(A,degA,B,degB)
+[rA,cA] = polsize(A,degA);
+[rB,cB] = polsize(B,degB);
+if cA ~= rB
+   error('polmul: Inconsistent dimensions of input matrices');
+end
+
+degC = degA+degB;
+C = [];
+for k = 0:degA+degB
+    mi = 0;
+    if k-degB > mi
+       mi = k-degB;
+    end
+    ma = degA;
+    if k < ma
+       ma = k;
+    end
+    Ck = zeros(rA,cB);
+    for i = mi:ma
+        Ck = Ck + A(:,i*cA+1:(i+1)*cA)*B(:,(k-i)*cB+1:(k-i+1)*cB);
+    end
+    C = [C Ck];
+end
+endfunction