1 files changed, 12 insertions, 0 deletions
diff --git a/2048/CH13/EX13.1/spec_ex.sce b/2048/CH13/EX13.1/spec_ex.sce
new file mode 100755
index 000000000..a120413fa
--- /dev/null
+++ b/2048/CH13/EX13.1/spec_ex.sce
@@ -0,0 +1,12 @@
+// Spectral factorization, as discussed in Example 13.3 on page 467.
+// 13.1
+
+exec('spec1.sci',-1);
+exec('flip.sci',-1);
+exec('polmul.sci',-1);
+exec('polsize.sci',-1);
+exec('poladd.sci',-1);
+
+A = convol([-0.5 1],[-0.9 1]); dA = 2;
+B = 0.5*[-0.9 1]; dB = 1; rho = 1;
+[r,beta1,sigma] = spec1(A,dA,B,dB,rho)