initial commit / add all books

1 files changed, 35 insertions, 0 deletions

diff --git a/2048/CH11/EX11.14/gmv_visc.sce b/2048/CH11/EX11.14/gmv_visc.sce
new file mode 100755
index 000000000..e8adeb79f
--- /dev/null
+++ b/2048/CH11/EX11.14/gmv_visc.sce
@@ -0,0 +1,35 @@
+// GMVC design of viscosity problem, as described in Example 11.10 on page 423.

+// 11.14

+

+// MacGregor's Viscosity control problem by gmv

+

+exec('gmv.sci',-1);

+exec('cl.sci',-1);

+exec('polmul.sci',-1);

+exec('polsize.sci',-1);

+exec('xdync.sci',-1);

+exec('rowjoin.sci',-1);

+exec('left_prm.sci',-1);

+exec('t1calc.sci',-1);

+exec('indep.sci',-1);

+exec('seshft.sci',-1);

+exec('makezero.sci',-1);

+exec('move_sci.sci',-1);

+exec('colsplit.sci',-1);

+exec('clcoef.sci',-1);

+exec('cindep.sci',-1);

+exec('poladd.sci',-1);

+exec('zpowk.sci',-1);

+exec('tfvar.sci',-1);

+exec('l2r.sci',-1);

+exec('transp.sci',-1);

+exec('tf.sci',-1);

+exec('covar_m.sci',-1);

+

+A = [1 -0.44]; dA = 1; B = [0.51 1.21]; dB = 1;

+C = [1 -0.44]; dC = 1; k = 1; int1 = 1;

+rho = 1;

+[Sc,dSc,R1,dR1] = gmv(A,dA,B,dB,C,dC,k,rho,int1);

+[Nu,dNu,Du,dDu,Ny,dNy,Dy,dDy,yvar,uvar] = ...

+        cl(A,dA,B,dB,C,dC,k,Sc,dSc,R1,dR1,int1);

+