1 files changed, 33 insertions, 0 deletions
diff --git a/1205/CH3/EX3.9/S_3_9.sce b/1205/CH3/EX3.9/S_3_9.sce
new file mode 100644
index 000000000..51a34d218
--- /dev/null
+++ b/1205/CH3/EX3.9/S_3_9.sce
@@ -0,0 +1,33 @@
+clc;
+// Resolving forces 125 N and 90  N
+f1=125;//N
+f2=90;//N
+f3=200;//N
+theta=45;// degree, angle made by f1 with x axis
+theta=theta*%pi/180;// Conversion of angle into radian
+alpha=30;// degree, angle made by force f2 with y axis
+alpha=alpha*%pi/180;// Conversion of angle into radian
+f1x=f1*cos(theta);//N, X component of 125 N.
+f1y=-f1*sin(theta);//N, Y component of 125 N
+
+f2x=-f2*sin(alpha);//N, X component of 90N
+f2y=-f2*cos(alpha);//N, Y component of 90N
+
+// Equivalent force at A
+Rx=f1x+f2x;//N
+Ry=f1y+f2y-f3;//N
+R=[Rx,Ry];//N
+R=norm(R);//N, magnitude of resultant
+theta=atan(Ry/Rx);//radian, angle of resultant with X axis
+theta=theta*180/%pi;
+
+printf("Equivalent force at A is R=%.2fN with angle %.2f Degrees \n ",R,theta);
+
+//equivalent couple at A
+//Clockwise moments are negative
+Meq=-550*f2*sin(35*%pi/180)-800*f3*sin(65*%pi/180)-1200*f1*sin(65*%pi/180);//N.mm , sum of all moments
+Meq=Meq/1000;// N.m conversion into N.m
+
+
+printf("Equivalent couple at A is Meq= %.2f N.m \n",Meq);
+