summaryrefslogtreecommitdiff
path: root/3792/CH6/EX6.14/Ex6_14.sce
diff options
context:
space:
mode:
Diffstat (limited to '3792/CH6/EX6.14/Ex6_14.sce')
-rw-r--r--3792/CH6/EX6.14/Ex6_14.sce22
1 files changed, 22 insertions, 0 deletions
diff --git a/3792/CH6/EX6.14/Ex6_14.sce b/3792/CH6/EX6.14/Ex6_14.sce
new file mode 100644
index 000000000..fdcb83c0a
--- /dev/null
+++ b/3792/CH6/EX6.14/Ex6_14.sce
@@ -0,0 +1,22 @@
+// SAMPLE PROBLEM 6/14
+clc;funcprot(0);
+// Given data
+// P=1.5*t;
+r_i=9/12;// ft
+r_o=18/12;// ft
+t_1=0;// s
+t_2=10;// s
+k=10/12;// ft
+W=120;// lb
+g=32.2;// The acceleration due to gravity in ft/sec^2
+v_1=-3;// ft/sec
+
+// Calculation
+function[X]=velocity(y)
+ X(1)=(((W/g)*v_1)+integrate('((1.5*t)-y(2))','t',t_1,t_2))-((W/g)*(r_o*y(1)));
+ X(2)=(((W/g)*(k)^2*(v_1/r_o))+integrate('((r_o*y(2))-(r_i*(1.5*t)))','t',t_1,t_2))-((W/g)*(k^2*y(1)));
+endfunction;
+y=[1 10];
+z=fsolve(y,velocity);
+omega_2=z(1);// rad/sec clockwise
+printf("\nThe angular velocity of the wheel,omega_2=%1.2f rad/sec",omega_2);