2 files changed, 23 insertions, 0 deletions
diff --git a/3875/CH1/EX1.3/1_3.txt b/3875/CH1/EX1.3/1_3.txt
new file mode 100644
index 000000000..c564134cb
--- /dev/null
+++ b/3875/CH1/EX1.3/1_3.txt
@@ -0,0 +1,2 @@
+Amplitude of oscillation = 1.26e-02 m
+Phase relative to the applied force is = 161.6 degree
+\ No newline at end of file
diff --git a/3875/CH1/EX1.3/Ex1_3.sce b/3875/CH1/EX1.3/Ex1_3.sce
new file mode 100644
index 000000000..55cb9d52e
--- /dev/null
+++ b/3875/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,21 @@
+clc ;
+clear ;
+m=0.1 // mass in kg
+K=100 //spring constant in N/m
+c=1 //resistive force in Nsm^-1
+F0=2 //force in N
+omega=50 //frequency in rad/s
+
+//calculation
+
+omega_n=sqrt(K/m) //in rad/s
+r=omega/omega_n 
+delta_st=F0/K //in m
+damp_ratio=c/(2*m*omega_n)
+A=delta_st/(sqrt((1-r^2)^2+(2*r*damp_ratio)^2))
+tan_phi=(2*r*damp_ratio)/(1-r^2) //in degree
+phi=180+atand(tan_phi) //converting degree to postive form
+
+mprintf("Amplitude of oscillation = %1.2e m\n",A)
+mprintf("Phase relative to the applied force is = %1.1f degree",phi) 
+//The answers vary due to round of errors