1 files changed, 50 insertions, 0 deletions
diff --git a/104/CH7/EX7.4/7_4.sce b/104/CH7/EX7.4/7_4.sce
new file mode 100755
index 000000000..a50229140
--- /dev/null
+++ b/104/CH7/EX7.4/7_4.sce
@@ -0,0 +1,50 @@
+//steady state errors from closed loop transfer functions
+s=%s
+p=poly([5 5 0],'s','coeff');
+q=poly([5 5 60 17 1],'s','coeff');
+M=p/q
+disp(M,"M(s)=")
+H=1;
+R=1;
+b=coeff(p)
+a=coeff(q)
+
+//step input
+if (a(1,1)==b(1,1)) then
+    printf("for unit step input Ess=0 \n" )
+else 
+    Ess=1/H*(1-(b(1,1)*H/a(1,1)))*R
+    disp(Ess,"for unit step 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("for unit ramp input Ess=0 \n")
+    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 Ess=inf \n")
+     end 
+ end
+
+//parabolic 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("for unit parabolic input Ess=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=inf \n")
+     end 
+ end
+\ No newline at end of file