1 files changed, 18 insertions, 0 deletions
diff --git a/1535/CH1/EX1.4/Ch01Ex4.sci b/1535/CH1/EX1.4/Ch01Ex4.sci
new file mode 100755
index 000000000..81515cfc0
--- /dev/null
+++ b/1535/CH1/EX1.4/Ch01Ex4.sci
@@ -0,0 +1,18 @@
+// Scilab Code Ex1.4 :  Page-24 (2010)
+v1 = 16;    // Velocity of particle executing SHM at position 3 cm
+v2 = 12;    // Velocity of particle executing SHM at position 4 cm
+x1 = 3;    // First position of the particle, cm
+x2 = 4;    // Second position of the particle, cm
+// As v = omega*sqrt(A^2 - x^2) so
+// (v1/v2)^2 = (A^2 - x1^2)/(A^2 - x2^2), solving for A
+A = poly(0, 'A');    // Declare variable A
+A = roots((A^2 - x1^2)*v2^2-(A^2 - x2^2)*v1^2);
+printf("\nThe amplitude of SHM = %1d cm", A(1));
+// v = omega*sqrt(A^2 - x^2), solving for omega
+omega = v1/sqrt(A(1)^2 - x1^2);    // Angular speed of the particle, rad per sec
+T = 2*%pi/omega;    // Time period of oscillation, sec
+printf("\nThe time period of oscillation = %5.3f sec", T);
+
+// Result 
+// The amplitude of SHM = 5 cm
+// The time period of oscillation = 1.571 sec