1 files changed, 30 insertions, 0 deletions
diff --git a/1205/CH18/EX18.6/S_18_6.sce b/1205/CH18/EX18.6/S_18_6.sce
new file mode 100644
index 000000000..ec68e9cd4
--- /dev/null
+++ b/1205/CH18/EX18.6/S_18_6.sce
@@ -0,0 +1,30 @@
+clc;
+m=300;//g, mass of each rod
+mom=1/1000;//given mo/m
+vo=2000;//m/s, relative velocity
+wo=60;//rpm
+wo=wo*2*%pi/60;//rad/s, conversion
+a=800;//mm, radius of disk
+a=a/1000;//m, conversion into meter
+//By theorytical work we get
+wx=-4/5*mom*vo/a;//rad/s , x component of angular velocity 
+wy=0;//rad/s , y component of angular velocity 
+wz=wo;//rad/s , x component of angular velocity 
+w=norm([wx,wy,wz]);//rad/s
+w=w/2/%pi*60;//rpm, conversion into rpm
+gama=atan(-wx/wz);//rad, 
+
+
+//precession axis
+
+theta=atan(2*mom*vo/a/wo);//rad, angle forme by precession axis and z axis
+
+//by law of sines
+phi=w*sin(gama)/sin(theta);// rpm, rates of precession
+psi=w*sin(theta-gama)/sin(theta);// rpm, rates of spin
+
+gama=gama*180/%pi;//degree, conversion into degree
+theta=theta*180/%pi;//degree, conversion into degree
+printf("The angular velocity after impact is w= %.1f rpm with angle gamma= %.1f degree\n",w,gama);
+printf("Angle formed by precession axis and z axis is theta= %.1f degree\n",theta);
+printf("Rate of precession is psi= %.1f rpm \n rate of spin is psi=%.1f rpm\n",phi,psi)