1 files changed, 46 insertions, 0 deletions
diff --git a/1151/CH3/EX3.1/test.sce b/1151/CH3/EX3.1/test.sce
new file mode 100755
index 000000000..a3b21be32
--- /dev/null
+++ b/1151/CH3/EX3.1/test.sce
@@ -0,0 +1,46 @@
+s=%s
+p= poly([0 0 50],'s','coeff');
+q= poly([0.1 2 10],'s','coeff');
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+     printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+     Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R
+disp (Ess ," f o r u n i t s t e p input Ess=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H
+disp (Ess ," for unit ramp input Ess=")
+ else printf (" for unit ramp input Es s=infinite\n")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" f o r u n i t p a r a b o l i c input Es s=0 \n")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H
+ disp (Ess ," for unit parabolic input Ess=")
+ else printf (" for unit parabolic  input Ess=infinite \n")
+ end
+ end