1 files changed, 21 insertions, 0 deletions
diff --git a/257/CH7/EX7.19/example_7_19.sce b/257/CH7/EX7.19/example_7_19.sce
new file mode 100644
index 000000000..3c923811e
--- /dev/null
+++ b/257/CH7/EX7.19/example_7_19.sce
@@ -0,0 +1,21 @@
+q=poly([0 20 1],'s','coeff');
+G=400/q //gain FACTOR=k
+H=1
+T=G/(1+G*H) 
+
+omegaN=sqrt(400)
+zeta=20/(2*omegaN)
+disp(omegaN,"omegaN = ")
+disp(zeta,"zeta = ")
+omegaD=omegaN*sqrt(1-zeta^2)
+theta= atan(sqrt(1-zeta^2)/zeta)
+disp(theta,"theta = ")
+
+syms s t
+c=(1-(%e^(-zeta*omegaN*t))/sqrt(1-zeta^2)*sin(omegaD*t+theta))
+disp(c, " c = ")
+
+Kv=limit(s*G*H,s,0)
+disp(Kv, " Kv = ")
+Ess=1/Kv
+disp(Ess, " Ess = ")