initial commit / add all books

29 files changed, 473 insertions, 0 deletions

diff --git a/3137/CH17/EX17.10/Ex17_10.sce b/3137/CH17/EX17.10/Ex17_10.sce
new file mode 100755
index 000000000..7eb700872
--- /dev/null
+++ b/3137/CH17/EX17.10/Ex17_10.sce
@@ -0,0 +1,8 @@
+//Initilization of variables

+speed=90000 //m/h

+P=100*1000 //N

+//Calculations

+Power=P*((speed)/3600) //J/s

+//Result

+clc

+printf('The power developed is %fJ/s',Power) 

diff --git a/3137/CH17/EX17.11/Ex17_11.sce b/3137/CH17/EX17.11/Ex17_11.sce
new file mode 100755
index 000000000..2b9b157ef
--- /dev/null
+++ b/3137/CH17/EX17.11/Ex17_11.sce
@@ -0,0 +1,16 @@
+//Initilization of variables

+d=0.6 //m

+T_t=800 //N

+T_s=180 //N

+w=200 //rpm

+//Calculations

+r=d/2 //m radius

+//Torque

+M=(T_t-T_s)*r //N.m

+//Power

+w_new=(2*%pi*w)/60 //rad/s

+Power=M*(w_new) //W

+//Result

+clc

+printf('The power transmitted is %f W',Power)

+

diff --git a/3137/CH17/EX17.12/Ex17_12.sce b/3137/CH17/EX17.12/Ex17_12.sce
new file mode 100755
index 000000000..25483a403
--- /dev/null
+++ b/3137/CH17/EX17.12/Ex17_12.sce
@@ -0,0 +1,12 @@
+//Initilization of variables

+P=25.6 //lb

+w=600 //rpm

+a=36 //in

+b=12 //in

+//Calculations

+M=P*(((b/2)+a)/12)  //lb-ft

+w_new=(2*%pi*w)/60 //rad/s

+Hp=(M*w_new)/550 //hp

+//Result

+clc

+printf('The power being transmitted is %fhp',Hp)

diff --git a/3137/CH17/EX17.13/Ex17_13.sce b/3137/CH17/EX17.13/Ex17_13.sce
new file mode 100755
index 000000000..1bf032fc0
--- /dev/null
+++ b/3137/CH17/EX17.13/Ex17_13.sce
@@ -0,0 +1,8 @@
+//Initilization of variables

+Pout=3.8 //bhp

+Pin=4.1 //ihp

+//Calculations

+Efficiency=round((Pout/Pin)*100) //Percent

+//Result

+clc

+printf('The efficiency of the engine is %ipercent',Efficiency)

diff --git a/3137/CH17/EX17.15/Ex17_15.sce b/3137/CH17/EX17.15/Ex17_15.sce
new file mode 100755
index 000000000..f587faf35
--- /dev/null
+++ b/3137/CH17/EX17.15/Ex17_15.sce
@@ -0,0 +1,20 @@
+//Initilization Of Variables

+a=3 //Lower Limit oF the Integral

+b=6 //Upper Limit of the Integral

+n=10 //Interval of the integral

+g=9.8 //m/s^2

+w=4/16

+//Calculation

+//Using Trapezoidal Rule for Intergration 

+function[I1]=Trap_Composite1(f,a,b,n)

+    h=(b-a)/n

+    x=linspace(a,b,n+1)

+    I1=(h/2)*((2*sum(f(x)))-f(x(1))-f(x(n+1)))

+endfunction

+deff('[y]=f(x)','y=-3*x^-1')

+an=-Trap_Composite1(f,a,b,n) //ft-lb

+v=sqrt((an*g)/(0.5*w)) //ft/s

+//Result

+clc

+printf('The speed of the disk is %fft/s',v)

+//The answer in the textbook is incorrect

diff --git a/3137/CH17/EX17.16/Ex17_16.sce b/3137/CH17/EX17.16/Ex17_16.sce
new file mode 100755
index 000000000..e627244be
--- /dev/null
+++ b/3137/CH17/EX17.16/Ex17_16.sce
@@ -0,0 +1,15 @@
+//Initilization of variables

+l=2 //m

+m=4 //kg

+w_1=20 //rpm

+w_2=50 //rpm

+rev=10 //no of revolution

+//Calculations

+Io=(1/3)*(m)*l^2 //kg.m^2

+w1=(2*%pi*w_1)/60 //rad/s

+w2=(2*%pi*w_2)/60 //rad/s

+theta=2*%pi*rev //rad

+M=(0.5*Io*(w2^2-w1^2))/theta //N.m

+//result

+clc

+printf('The constant moment required is %fN.m',M)

diff --git a/3137/CH17/EX17.18/Ex17_18.sce b/3137/CH17/EX17.18/Ex17_18.sce
new file mode 100755
index 000000000..323132044
--- /dev/null
+++ b/3137/CH17/EX17.18/Ex17_18.sce
@@ -0,0 +1,13 @@
+//Initilization of variables

+W=1000 //lb

+w_w=200 //lb weight of the individual wheel

+d_w=2.5 //ft diameter of the wheel

+v=22 //ft/s

+t=2 //minutes

+//Calculations

+//T1=Initial Kinetic Energy and T2=Final Kinetic Energy

+F=(-0.5*W*32.2^-1*v^2-4*0.5*w_w*32.2^-1*(v^2+0.5*v^2))/(10560) //lb

+//Negative sign in the answer tells it oposses the motion 

+//Result

+clc

+printf('The rolling resistance is %flb',F) 

diff --git a/3137/CH17/EX17.19/Ex17_19.sce b/3137/CH17/EX17.19/Ex17_19.sce
new file mode 100755
index 000000000..c2dfa8a6e
--- /dev/null
+++ b/3137/CH17/EX17.19/Ex17_19.sce
@@ -0,0 +1,22 @@
+//Initilization of variables

+W=100 //lb

+lo=4 //ft

+theta=45 //degrees

+g=32.2 //ft/s^2

+l=8/3 //ft

+//Calculations

+//Taking moment about point O and equating it to zero

+alpha=(W*(lo*0.5)*cosd(theta))/((W/g)*(l)*2) //rad/s^2

+//Summing forces in the t direction

+Ot=(W*cosd(theta))-((W/g)*lo*0.5*alpha) //lb

+//Work Done

+Work=W*(lo*0.5*cosd(theta)) //ft/lb

+//Moment of inertia

+Io=(1/3)*(W/g)*(lo^2) //kg-ft^2

+//Using the concept for work done=chane in K.E

+w=sqrt(Work/(0.5*Io)) //rad/s

+//Summing forces along the bar

+On=-(-((W/g)*lo*0.5*w^2)-(W*cosd(theta))) //lb

+//Result

+clc

+printf('The bearing reaction at O on the rod is %flb',On)

diff --git a/3137/CH17/EX17.21/Ex17_21.sce b/3137/CH17/EX17.21/Ex17_21.sce
new file mode 100755
index 000000000..54ec01987
--- /dev/null
+++ b/3137/CH17/EX17.21/Ex17_21.sce
@@ -0,0 +1,10 @@
+//Initilization of variables

+vo=9 //m/s

+theta=30 //degrees

+g=9.8 //m/s^2

+//Calculations

+x=((7/10)*vo^2)/(g*sind(theta)) //m

+//Result

+clc

+printf('The ball will roll %f m up the plane',x)

+//The textbook wrongly mentions the unit of displacement as in 

diff --git a/3137/CH17/EX17.22/Ex17_22.sce b/3137/CH17/EX17.22/Ex17_22.sce
new file mode 100755
index 000000000..d025c6222
--- /dev/null
+++ b/3137/CH17/EX17.22/Ex17_22.sce
@@ -0,0 +1,18 @@
+//Initilization of variables

+W=322 //lb

+F=12 //lb

+a=0 //lower limit (where the cyliner starts rolling)

+b=%pi/2 //Upper Limit (where the cyliner stops rolling)

+d=3.2 //ft

+g=32.2 //ft/s^2

+//Calculations

+dR=1.6 //Differential Radius

+d_U=2*dR*F //differential work done

+//Integration Calculations

+//As it is a simple integration we can resort to this

+U=d_U*(b-a) //ft-lb

+//Determination of K.E

+w=sqrt(U/((0.5*(W/g)*(1/(d/2)^2))+((0.5*0.5)*(W/g)*(d/2)^2))) //rad/s

+//Result  

+clc

+printf('the angular velocity of the is %f rad/s',w)

diff --git a/3137/CH17/EX17.23/Ex17_23.sce b/3137/CH17/EX17.23/Ex17_23.sce
new file mode 100644
index 000000000..5769fad12
--- /dev/null
+++ b/3137/CH17/EX17.23/Ex17_23.sce
@@ -0,0 +1,19 @@
+//Initilization of variables

+m=90 //kg

+k=450 //N/m

+lo=0.6 //m

+r=0.15 //m

+x=0.9 //m

+y=0.4 //m

+//Calculations

+//Initial KE=0

+I=0.5*m*r^2 //kg.m^2

+s1=sqrt((lo^2)+(x^2)) //m

+s2=sqrt((lo^2)+(y^2)) //m

+V1=0.5*k*(s1-lo)^2 //N.m

+V2=0.5*k*(s2-lo)^2 //N.m

+//Applying Conservation of Energy 

+w=sqrt((V1-V2)/(0.5*m*r^2+0.5*I)) //rad/s

+//Result

+clc

+printf('The value of  angular speed is: %f rad/s',w)

diff --git a/3137/CH17/EX17.24/Ex17_24.sce b/3137/CH17/EX17.24/Ex17_24.sce
new file mode 100755
index 000000000..03cb1c69b
--- /dev/null
+++ b/3137/CH17/EX17.24/Ex17_24.sce
@@ -0,0 +1,21 @@
+//Initilization of variables

+Wa=161 //lb

+Wb=193.2 //lb

+Wc=322 //lb

+v1=5 //ft/s

+lc=6 //in

+k=6 //lb/ft

+l=4 //ft

+u=0.2 //coefficient of friction

+g=32.2 //ft/s^2

+//Calculations

+Ib=(1/2)*(Wb/g)*(1/2)^2 //Moment of inertia

+w1=v1/0.5 //rad/s

+T1=(0.5*(Wc/g)*v1^2)+(0.5*Ib*w1^2)+(0.5*(Wa/g)*v1^2) //ft-lb

+//Work Done on the system

+//The textbook is ambigious on the calculations hence the result is dispalyed directly

+U=26.4 //ft-lb

+//Velocity Calculations

+v=sqrt((T1+U)/9) //ft/s

+//Result

+printf('The velocity of the block is %f',v)

diff --git a/3137/CH17/EX17.25/Ex17_25.sce b/3137/CH17/EX17.25/Ex17_25.sce
new file mode 100755
index 000000000..01f9df76f
--- /dev/null
+++ b/3137/CH17/EX17.25/Ex17_25.sce
@@ -0,0 +1,14 @@
+//Initilization of variables

+Mm=70 //kg

+Mc=45 //kg

+R=0.6 //m

+g=9.8 //m/s^2

+l=5 //m

+theta=50 //degrees

+//Calculations

+//T2 calculations except for v term in it as it cannot be declared as a number

+T2=68.7 //without the v term in it

+v=sqrt((g*Mm*l-g*Mc*l*sind(theta))/T2) //m/s

+//Result

+clc

+printf('The speed is %fm/s',v)

diff --git a/3137/CH17/EX17.26/Ex17_26.sce b/3137/CH17/EX17.26/Ex17_26.sce
new file mode 100755
index 000000000..9ed79b372
--- /dev/null
+++ b/3137/CH17/EX17.26/Ex17_26.sce
@@ -0,0 +1,17 @@
+//The textbook has a typo in printing the question number

+//Initilization of variables

+W1=96.6 //lb

+W2=128.8 //lb

+v=8 //ft/s

+g=32.2 //ft/s^2

+theta=30 //degrees

+//Calculations

+//Initial KE of the system is T1=0

+T2=(0.5*(W1/g)*v^2)+(0.5*(W2/g)*(v/2)^2) //ft-lb

+//Work Done without s term

+U=-(W1*sind(theta))+W2*0.5

+//S calculations

+s=T2/U //ft

+//Result

+clc

+printf('The block attains a speed of 8ft/s in %f ft',s)

diff --git a/3137/CH17/EX17.28/Ex17_28.sce b/3137/CH17/EX17.28/Ex17_28.sce
new file mode 100755
index 000000000..1a1513b2c
--- /dev/null
+++ b/3137/CH17/EX17.28/Ex17_28.sce
@@ -0,0 +1,20 @@
+//Initilization Of Variables

+a=0 //Lower Limit oF the Integral

+b=6 //Upper Limit of the Integral

+n=10 //Interval of the integral

+m=50 //kg

+l=6 //m

+g=9.8 //m/s^2

+//Calculation

+//Gravatational Force is 

+Fg=g*(m/l) //dx

+//Using Trapezoidal Rule for Intergration 

+function[I1]=Trap_Composite1(f,a,b,n)

+    h=(b-a)/n

+    t=linspace(a,b,n+1)

+    I1=(h/2)*((2*sum(f(t)))-f(t(1))-f(t(n+1)))

+endfunction

+deff('[y]=f(t)','y=Fg*(6-t)') 

+//Result

+clc

+printf('The Work done is %f N.m',Trap_Composite1(f,a,b,n))

diff --git a/3137/CH17/EX17.31/Ex17_31.sci b/3137/CH17/EX17.31/Ex17_31.sci
new file mode 100755
index 000000000..77d798fc9
--- /dev/null
+++ b/3137/CH17/EX17.31/Ex17_31.sci
@@ -0,0 +1,18 @@
+//Initilization Of Variables

+x1=150 //mm

+x2=450 //mm

+a=0 //Lower Limit oF the Integral

+b=(x2-x1) //Upper Limit of the Integral

+n=10 //Interval of the integral

+k=0.044 //N/m

+//Calculation

+//Using Trapezoidal Rule for Intergration 

+function[I1]=Trap_Composite1(f,a,b,n)

+    h=(b-a)/n

+    t=linspace(a,b,n+1)

+    I1=(h/2)*((2*sum(f(t)))-f(t(1))-f(t(n+1)))

+endfunction

+deff('[y]=f(t)','y=k*t') 

+//Result

+clc

+printf('The Work done is %f N.m',Trap_Composite1(f,a,b,n)/1000)

diff --git a/3137/CH17/EX17.32/Ex17_32.sce b/3137/CH17/EX17.32/Ex17_32.sce
new file mode 100755
index 000000000..f1f2ab3c2
--- /dev/null
+++ b/3137/CH17/EX17.32/Ex17_32.sce
@@ -0,0 +1,15 @@
+//Initilization of variables

+m=10 //kg

+d=1.2 //m

+g=9.8 //m/s^2

+//Calculations

+//Initilial KE is zero

+//Final KE is(without v^2 term in it)

+KE2=(3/4)*10

+//Work Done

+U=m*g*d //N.m

+//Velocity calculations

+v=sqrt(U/KE2) //m/s

+//Result

+clc

+printf('The velocity is %fm/s',v)

diff --git a/3137/CH17/EX17.33/Ex17_33.sce b/3137/CH17/EX17.33/Ex17_33.sce
new file mode 100755
index 000000000..78cc1b20b
--- /dev/null
+++ b/3137/CH17/EX17.33/Ex17_33.sce
@@ -0,0 +1,27 @@
+//Initilization of variables

+W=161 //lb

+wa=150 //lb

+wb=100 //lb

+la=2 //ft

+lb=4 //ft

+//Calculations

+//Work Done

+T1=wb*lb-wa*la //ft-lb

+//Final KE=zero

+T2=0 //ft-lb

+//Work Done on the system=T2-T1

+//Hence the equation becomes

+//50x-50x^2+100=0

+//where

+a=-50

+b=50

+c=100

+//Solution

+d=sqrt(b^2-4*a*c) 

+x1=(-b+d)/(2*a) //ft

+x2=(-b-d)/(2*a) //ft

+//Result

+clc

+printf('The stretch of the spring is %f',x2)

+//Here even x1 could have been the solution,but the stretch in the string is elongation not compression hence x2 is the valid answer

+

diff --git a/3137/CH17/EX17.34/Ex17_34.sce b/3137/CH17/EX17.34/Ex17_34.sce
new file mode 100755
index 000000000..d4837a8ff
--- /dev/null
+++ b/3137/CH17/EX17.34/Ex17_34.sce
@@ -0,0 +1,13 @@
+//Initilization of variables

+I=100 //slug-ft^2

+w=4 //rad/s

+theta=6 //rad

+Mc=64.4 //lb

+g=32.2 //ft/s^2

+//Calculations

+vb=2*w //ft/s

+vc=0.5*w //ft/s

+Mb=(0.5*I*w^2+0.5*(Mc/g)*vc^2+0.5*Mc*theta)/(2*theta-(0.5*vb^2*(1/g))) //lb

+//Result

+clc

+printf('The weight of the block B is %f lb',Mb)

diff --git a/3137/CH17/EX17.35/Ex17_35.sce b/3137/CH17/EX17.35/Ex17_35.sce
new file mode 100755
index 000000000..9fe3257a8
--- /dev/null
+++ b/3137/CH17/EX17.35/Ex17_35.sce
@@ -0,0 +1,21 @@
+//Initilization of variables

+Wa=96.6 //lb

+Wb=128.8 //lb

+g=32.2 //ft/s^2

+I=12 //slug-ft^2

+v=16 //ft/s

+ratio=1/3 //ratio of Sb/Sa

+r=3//ft

+va=6 //ft/s

+vb=2 //ft/s

+//Calculations

+//Work Done without S in it

+W=Wa-(ratio*Wb)

+//System has zero KE initially and final KE is given by

+w=va/r //rad/s

+T2=(0.5*(Wa/g)*va^2+0.5*I*w^2+0.5*(Wb/g)*vb^2) //ft-lb

+//Distance Calculations

+S=T2/W //ft

+//Result

+clc

+printf('The distance through which A falls is %f ft',S)

diff --git a/3137/CH17/EX17.36/Ex17_36.sce b/3137/CH17/EX17.36/Ex17_36.sce
new file mode 100755
index 000000000..2e6bc13c4
--- /dev/null
+++ b/3137/CH17/EX17.36/Ex17_36.sce
@@ -0,0 +1,21 @@
+//initilization of variables

+u=0.25 //coefficient of friction

+k=2800 //N/m

+x=0.075 //m

+g=9.8 //m/s^2

+m=7 //kg

+theta=30 //degrees

+//Calculations

+//Normal Reaction

+N=g*m*cosd(theta) //N

+//Frictional Force

+Fr=u*N //N

+//Component of force along the plane

+F=g*m*sind(theta) //N

+//Spring work is

+W=0.5*k*x*x //N.m

+s=(W+Fr*x-F*x)/(F-Fr) //m

+S=round(s*1000) //mm

+//Result

+clc

+printf('The value of S is %i mm',S) 

diff --git a/3137/CH17/EX17.37/Ex17_37.sce b/3137/CH17/EX17.37/Ex17_37.sce
new file mode 100755
index 000000000..22ffc4465
--- /dev/null
+++ b/3137/CH17/EX17.37/Ex17_37.sce
@@ -0,0 +1,19 @@
+//Initilization of variables

+m=5 //kg

+l=2 //m

+k=10000 //N/m

+x=0.1 //m

+g=9.8 //m/s^2

+//Calculations

+drop=l+x //m mass drop length

+//Work Done by Gravity

+Wg=g*m*drop //N.m

+//Work Done by Spring

+Ws=0.5*k*x^2 //N.m

+//Increase in KE is without v^2

+KE=0.5*m //kg

+//Velocity Calculations

+v=sqrt((Wg-Ws)/KE) //m/s

+//Result

+clc

+printf('The speed is %f m/s',v)

diff --git a/3137/CH17/EX17.38/Ex17_38.sce b/3137/CH17/EX17.38/Ex17_38.sce
new file mode 100755
index 000000000..44844b862
--- /dev/null
+++ b/3137/CH17/EX17.38/Ex17_38.sce
@@ -0,0 +1,14 @@
+//Initilization of variables

+l=6 //ft

+k=20 //lb/in

+x=8 //in

+//Calculations

+//Work Done by Gravity

+Wg=(l*12+x) //in without W

+//Work Done by Spring

+Ws=0.5*k*x^2 //in-lb

+//Change in the kinetic energy is zero

+W=Ws/Wg //lb

+//Result

+clc

+printf('The weight is %i lb',W) 

diff --git a/3137/CH17/EX17.4/Ex17_4.sci b/3137/CH17/EX17.4/Ex17_4.sci
new file mode 100755
index 000000000..20cabc3fa
--- /dev/null
+++ b/3137/CH17/EX17.4/Ex17_4.sci
@@ -0,0 +1,16 @@
+//Initilization Of Variables

+s1=2 //Lower Limit oF the Integral

+s2=5 //Upper Limit of the Integral

+n=10 //Interval of the integral

+k=20 //lb/in

+//Calculation

+//Using Trapezoidal Rule for Intergration 

+function[I1]=Trap_Composite1(f,s1,s2,n)

+    h=(s2-s1)/n

+    s=linspace(s1,s2,n+1)

+    I1=(h/2)*((2*sum(f(s)))-f(s(1))-f(s(n+1)))

+endfunction

+deff('[y]=f(s)','y=k*s')

+//Result

+clc

+printf('The work done is %f in-lb',Trap_Composite1(f,s1,s2,n) )

diff --git a/3137/CH17/EX17.40/Ex17_40.sce b/3137/CH17/EX17.40/Ex17_40.sce
new file mode 100755
index 000000000..852096309
--- /dev/null
+++ b/3137/CH17/EX17.40/Ex17_40.sce
@@ -0,0 +1,12 @@
+//Initilization of variables

+W=8 //lb

+//Calculations

+//work done by the spring woithout k

+Ws=0.5*((9/12)^2-(1/12)^2) 

+//Work done by gravity

+Wg=W*(10.5/12) //ft-lb

+//Change in KE is zero

+k=Wg/Ws //lb/ft

+//Result

+clc

+printf('The value of k is %f lb/ft',k)

diff --git a/3137/CH17/EX17.41/Ex17_41.sce b/3137/CH17/EX17.41/Ex17_41.sce
new file mode 100755
index 000000000..bb104a5c2
--- /dev/null
+++ b/3137/CH17/EX17.41/Ex17_41.sce
@@ -0,0 +1,15 @@
+//Initilization of variables

+Wc=100 //lb

+r= 1 //ft

+F=80 //lb

+k=50 //lb/ft

+s=6 //in

+g=32.2 //ft/s^2

+//Calculations

+//Work done on the system

+U=-0.5*k*(1)+F*(s/12) //ft-lb

+//Initial KE is zero

+Vo=sqrt(U/(0.5*(Wc/g+0.5*(Wc/g)*r))) //ft/s

+//Result

+clc

+printf('The initial speed is %f ft/s',Vo)

diff --git a/3137/CH17/EX17.6/Ex17_6.sci b/3137/CH17/EX17.6/Ex17_6.sci
new file mode 100755
index 000000000..f35d9b4aa
--- /dev/null
+++ b/3137/CH17/EX17.6/Ex17_6.sci
@@ -0,0 +1,21 @@
+//Initilization of variables

+m=5 //kg

+d=6 //m

+theta1=30 //degrees

+theta2=10 //degrees

+u=0.2 //coefficient of friction

+g=9.8 //m/s^2

+F=70 //N

+//Calculations

+//Using free body diagram

+Na=m*g*cosd(theta1)-(F*sind(theta2)) //N

+//work done by each force

+W=[F*cosd(theta2),-m*g*sind(theta1),0,-u*Na*d] //N.m

+//Total Work Done

+W_tot=W(1)+W(2)+W(3)+W(4) //N.m

+//Using resultant

+R=F*cosd(theta2)-(u*Na)-(m*g*sind(theta1)) //N

+W_d=R*d //N.m (Work Done)

+//Result

+clc

+printf('The work done is %f N.m',W_d)

diff --git a/3137/CH17/EX17.7/Ex17_7.sce b/3137/CH17/EX17.7/Ex17_7.sce
new file mode 100755
index 000000000..d60d79c69
--- /dev/null
+++ b/3137/CH17/EX17.7/Ex17_7.sce
@@ -0,0 +1,12 @@
+//Initilization of variables

+m=20 //kg

+d=1.5 //m

+theta=30 //degrees

+u=0.25 //coefficient of friction

+g=9.8 //m/s^2

+F=130 //N

+//Calculations

+W=F*d-(m*g*sind(theta)*d) //N.m

+//Result

+clc

+printf('The work done is %i N.m',W)

diff --git a/3137/CH17/EX17.9/Ex17_9.sce b/3137/CH17/EX17.9/Ex17_9.sce
new file mode 100755
index 000000000..1988e6bcd
--- /dev/null
+++ b/3137/CH17/EX17.9/Ex17_9.sce
@@ -0,0 +1,16 @@
+//Initilization of variables

+d=6/12 //ft

+l=8/12 //ft

+l_c=3.2 //in

+y=1.82 //in^2

+//Calculations

+V=(1/4)*%pi*d^2*l //ft^3

+//One horizontal inch 

+h_i=V/l_c //ft^3

+//One vertical inch

+v_i=100*144 //lb/ft^2

+//Then 1.82 in^2 represents

+x=y*v_i*h_i //ft-lb

+//Result

+clc

+printf('The work capacity is %f ft-lb',x)