From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 2048/CH11/EX11.14/gmv_visc.sce | 35 +++++++++++++++++++++++++++++++++++ 1 file changed, 35 insertions(+) create mode 100755 2048/CH11/EX11.14/gmv_visc.sce (limited to '2048/CH11/EX11.14') 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); + -- cgit