From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3792/CH3/EX3.30/Ex3_30.sce | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100644 3792/CH3/EX3.30/Ex3_30.sce

(limited to '3792/CH3/EX3.30/Ex3_30.sce')

diff --git a/3792/CH3/EX3.30/Ex3_30.sce b/3792/CH3/EX3.30/Ex3_30.sce
new file mode 100644
index 000000000..67df7f584
--- /dev/null
+++ b/3792/CH3/EX3.30/Ex3_30.sce
@@ -0,0 +1,32 @@
+// SAMPLE PROBLEM 3/30
+clc;funcprot(0);
+// Given data
+v_1=6;// m/s
+v_2=0;// m/s
+e=0.6;// The coefficient-of-restitution
+theta=30;// degree
+
+// Calculation
+// Assume a for '
+v_1n=v_1*cosd(theta);// m/s
+v_1t=v_1*sind(theta);// m/s
+v_2n=0;// m/s
+v_2t=v_2n;// m/s
+function[X]=velocity(y)
+    X(1)=(v_1n+v_2n)-(y(1)+y(2));
+    X(2)=(e*(v_1n+v_2n))-(y(2)-y(1));
+endfunction
+y=[1,1];
+z=fsolve(y,velocity);
+v_1an=z(1);// m/s
+v_2an=z(2);// m/s
+v_1at=v_1t;// m/s
+v_2at=v_2t;// m/s
+v_1a=sqrt((v_1an)^2+(v_1at)^2);// m/s
+v_2a=sqrt((v_2an)^2+(v_2at)^2);// m/s
+thetaa=atand(v_1an/v_1at);// m/s
+// The kinetic energies just before and just after impact, with m=m1=m2,are
+T=18;// m
+T_a=13.68;// m
+E_l=((T-T_a)/T)*100;// The percentage energy loss(%)
+printf("\nThe final speeds of the particles v_1a=%1.2f m/s ,v_2a=%1.2f m/s \nThe angle which v_1a makes with the t-direction,theta=%2.2f degree \nThe percentage energy loss is %2.0f percentage.",v_1a,v_2a,thetaa,E_l);
-- 
cgit