updated the code

1 files changed, 11 insertions, 10 deletions

diff --git a/22/CH4/EX4.27/ch4ex27.sce b/22/CH4/EX4.27/ch4ex27.sce
index 0389caa79..85993f0a8 100755
--- a/22/CH4/EX4.27/ch4ex27.sce
+++ b/22/CH4/EX4.27/ch4ex27.sce
@@ -1,10 +1,11 @@
-omega_0=2*%pi*60; theta = [60 80 87]*(%pi/180);

-omega = (0:0.5:1000)'; mag = zeros(3,length(omega));

-s=poly(0,'s')

-for m =1:length(theta)

-    H=syslin('c',((s^2+omega_0^2)/(s^2+2*omega_0*cos(theta(m))*s +omega_0^2)));

-    bode(H,10,100);

-end

-f=omega/((2*%pi))plot(f,mag(1,:),'k-',f mag(2,:),'k--',f,mag(3,:),'k-.');

-xlabel('f[hz]'); ylabel('|H(j2/pi f)|');

-legend('\theta=60^\circ','\theta = 80^\circ','\theta = 87^\circ',0)
\ No newline at end of file
+omega_0=2*%pi*60; theta = [60 80 87]*(%pi/180);
+omega = (0:0.5:1000)'; mag = zeros(3,length(omega));
+s=poly(0,'s')
+for m =1:length(theta)
+    H=syslin('c',((s^2+omega_0^2)/(s^2+2*omega_0*cos(theta(m))*s +omega_0^2)));
+    bode(H,10,100);
+end
+f=omega/((2*%pi));
+plot(f,mag(1,:)','k-',f,mag(2,:)','k--',f,mag(3,:)','k-.');
+xlabel('f[hz]'); ylabel('|H(j2/pi f)|');
+legend('\theta=60^\circ','\theta = 80^\circ','\theta = 87^\circ',1)
\ No newline at end of file