1 files changed, 18 insertions, 0 deletions
diff --git a/1205/CH10/EX10.4/S_10_4.sce b/1205/CH10/EX10.4/S_10_4.sce
new file mode 100644
index 000000000..b6e80d1e5
--- /dev/null
+++ b/1205/CH10/EX10.4/S_10_4.sce
@@ -0,0 +1,18 @@
+clc;
+m=10;//kg, mass of rim
+r=300;//mm, radius of disk
+a=0.08;//m
+b=0.3;//m
+k=4;//kN/m
+g=9.81;//m/s^2  gravity
+//From theory we get
+
+//sin(theta)=k*a^2/m/g/b*theta
+dif=1;
+for theta=0:0.001:1
+    dif=sin(theta)-k*a^2/m/g/b*theta;
+    if dif<=0.001 then printf("theta= %.3f rad or %.1f degrees\n",theta,theta/%pi*180);
+    end
+end
+
+