242 files changed, 2559 insertions, 0 deletions
diff --git a/564/CH1/EX1.1/1_01.sce b/564/CH1/EX1.1/1_01.sce
new file mode 100755
index 000000000..7907f07d6
--- /dev/null
+++ b/564/CH1/EX1.1/1_01.sce
@@ -0,0 +1,18 @@
+pathname=get_absolute_file_path('1_01.sce')
+filename=pathname+filesep()+'1_01data.sci'
+exec(filename)
+Sx1= p*d/(4*t);
+Sy= p*d/(2*t);//σy
+printf("\nLongitudnal stress produced by internal pressure: %f N/mm^2",Sx1);
+printf("\nCircumferential stress produced by internal pressure (σy): %f N/mm^2",Sy);
+Sx2= Load/(%pi*d*t);
+printf("\ndirect stress due to the axial load: %f N/mm^2",Sx2);
+Sx= Sx1+Sx2;//σx
+printf("\nσx is: %f N/mm^2",Sx);
+printf("\nσy is: %f N/mm^2",Sy);
+Sn=(Sx)*(cos(theta))^2 + Sy*(sin(theta))^2;//σn
+T= (Sx-Sy)*sin(2*(theta)) /2;//τ
+printf("\nσn is: %f N/mm^2",Sn);
+printf("\nτ is: %f N/mm^2",T);
+Tmax= (Sx-Sy)/2;//τ max
+printf("\nτmax: %f N/mm^2",Tmax);
+\ No newline at end of file
diff --git a/564/CH1/EX1.2/1_02.sce b/564/CH1/EX1.2/1_02.sce
new file mode 100755
index 000000000..a4bf98e69
--- /dev/null
+++ b/564/CH1/EX1.2/1_02.sce
@@ -0,0 +1,18 @@
+pathname=get_absolute_file_path('1_02.sce')
+filename=pathname+filesep()+'1_02data.sci'
+exec(filename)
+printf("\nbending moment due to direct loading in a vertical plane: %f N/mm^2",Load*t);
+Sx1= Load/(%pi*(d^2)/4);//σx (axial load)
+printf("\nσx (axial load) is: %f N/mm^2",Sx1);
+I= %pi*(d^4)/64;// moment of Inertia
+Sx2= t*Load*(d/2)*(1/I);//σx (bending moment)
+printf("\nσx (bending moment) is: %f N/mm^2",Sx2)
+Sx=Sx1+Sx2;//σx 
+J= %pi*(d^4)/32;//torsion constant
+Txy=T*(d/2)*(1/J);//τxy
+printf("\nSince the element is positioned at the bottom of the beam\nτxy:%fN/mm^2",-Txy);
+printf("\nσx: %f N/mm^2",-Sx);
+Sn=(-Sx)*(cos(theta))^2 +(-Txy)*sin(2*theta);//σn
+T= (-Sx-0)*sin(2*theta)/2 -(-Txy)*cos(2*theta);//τ
+printf("\nσn: %f N/mm^2",Sn);
+printf("\n τ: %f N/mm^2)",T);
+\ No newline at end of file
diff --git a/564/CH1/EX1.3/1_03.sce b/564/CH1/EX1.3/1_03.sce
new file mode 100755
index 000000000..3f913c5f6
--- /dev/null
+++ b/564/CH1/EX1.3/1_03.sce
@@ -0,0 +1,27 @@
+pathname=get_absolute_file_path('1_03.sce')
+filename=pathname+filesep()+'1_03data.sci'
+exec(filename)
+Txy=((lLoad-Sx)*(lLoad-Sy))^0.5;
+printf("\nτxy:%f N/mm^2",Txy);
+printf(",%f N/mm^2",-Txy);
+S=poly(0,'S');
+y=(S^2)-(S*(Sx+Sy))+(Sx*Sy)-(Txy^2);
+m=roots(y);
+printf("\nσ1: %f N/mm^2",m(2,:));
+printf("\nσ2:%f N/mm^2",m(1,:));
+Tmax=abs((m(2,:)-m(1,:))/2);
+printf("\nτmax: %f N/mm^2",Tmax);
+//plotting mohr circle
+x=abs((m(2,:)+m(1,:))/2);//Centre of circle
+plot2d(0,0,-1,"031"," ",[x-Tmax-50,-Tmax-50,x+Tmax+100,Tmax+50]);
+xgrid(3);
+xarc(x-Tmax,Tmax,2*Tmax,2*Tmax,0,360*64);
+x1=[m(1,:),m(2,:)],y1=[0,0];
+x2=[Sx,Sy],y2=[Txy,-Txy];
+x3=[Sx,Sx],y3=[Txy,0];
+x4=[Sy,Sy],y4=[-Txy,0];
+x5=[x,x],y5=[Tmax,-Tmax];
+plot(x1,y1,x2,y2,'--',x3,y3,'--',x4,y4,'--',x5,y5,'--');
+datatipToggle();
+printf("\n\nclick on the point to view its coordinate on the plot");
+xtitle( 'Mohr’s circle of stress', ' σ (N/mm^2)', ' τ (N/mm^2)', boxed = 1 );
+\ No newline at end of file
diff --git a/564/CH1/EX1.3/1_03MohrCircle.JPG b/564/CH1/EX1.3/1_03MohrCircle.JPG
new file mode 100755
index 000000000..101cd4dca
--- /dev/null
+++ b/564/CH1/EX1.3/1_03MohrCircle.JPG
diff --git a/564/CH1/EX1.4/1_04.sce b/564/CH1/EX1.4/1_04.sce
new file mode 100755
index 000000000..7100dce78
--- /dev/null
+++ b/564/CH1/EX1.4/1_04.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('1_04.sce')
+filename=pathname+filesep()+'1_04data.sci'
+exec(filename)
+Ex= (1/E)*(Sx-v*Sy);//εx
+Ey= (1/E)*(Sy-v*Sx);//εy
+Ez= (1/E)*(0-v*(Sx+Sy));//εz
+printf("\nεx: %f",Ex);
+printf("\nεy: %f",Ey);
+printf("\nεz: %f",Ez);
+Tmax=(Sx-Sy)/2;//τmax
+Gmax=2*(1+v)*Tmax/E;//γmax
+printf("\nγmax: %f",Gmax);
+printf("\nτmax: %f N/mm^2",Tmax);
+\ No newline at end of file
diff --git a/564/CH1/EX1.5/1_05.sce b/564/CH1/EX1.5/1_05.sce
new file mode 100755
index 000000000..f44bb5b3f
--- /dev/null
+++ b/564/CH1/EX1.5/1_05.sce
@@ -0,0 +1,34 @@
+pathname=get_absolute_file_path('1_05.sce')
+filename=pathname+filesep()+'1_05data.sci'
+exec(filename)
+Ex= (1/E)*(Sx-v*Sy);//εx
+Ey= (1/E)*(Sy-v*Sx);//εy
+G=E/(2*(1+v));//shear modulus
+Gxy=Txy/G;//γxy
+printf("\n εx:  %f",Ex);
+printf("\n εy: %f",Ey);
+printf("\n G: %f",G);
+printf("\nγxy:  %f",Gxy);
+PE1=(Ex+Ey)/2 + 0.5*((Ex-Ey)^2 +Gxy^2)^0.5;//εI
+PE2=(Ex+Ey)/2 - 0.5*((Ex-Ey)^2 +Gxy^2)^0.5;//εII
+theta=0.5*atan(Gxy/(Ex-Ey)) *(180/%pi);//θ in degree
+printf("\n εI:  %f",PE1);
+printf("\nεII: %f",PE2);
+printf("\n  θ:  %f degree",theta);
+printf("\n  θ:  %f degree",theta+90);
+//plotting mohr circle
+x=abs((PE1+PE2)/2);//Centre
+r=(abs(PE1)+abs(PE2))/2;//radius
+plot2d(0,0,-1,"031"," ",[x- 1.5*r,-1.5*r,x+ 1.5*r,1.5*r]);
+xgrid(3);
+xarc(x-r,r,2*r,2*r,0,360*64);
+x1=[PE1,PE2],y1=[0,0];
+x2=[Ex,Ey],y2=0.5*[Gxy,-Gxy];
+theta1=atan(Gxy/(2*Ex)) *(180/%pi);
+plot(x1,y1);
+plot(x2,y2);
+xarc(x-r/4,r/4,r/2,r/2,0,theta1*64);
+xarc(x-r/2,r/2,r,r,180*64,(180+theta1)*64);
+datatipToggle();
+printf("\n\nclick on the point to view its coordinate on the plot");
+xtitle( 'Mohr’s circle of strain', ' ε', ' γ (gamma)', boxed = 1 );
+\ No newline at end of file
diff --git a/564/CH1/EX1.5/1_05mohrcircle.JPG b/564/CH1/EX1.5/1_05mohrcircle.JPG
new file mode 100755
index 000000000..7f24c7de5
--- /dev/null
+++ b/564/CH1/EX1.5/1_05mohrcircle.JPG
diff --git a/564/CH1/EX1.7/1_07.sce b/564/CH1/EX1.7/1_07.sce
new file mode 100755
index 000000000..52a9a57f1
--- /dev/null
+++ b/564/CH1/EX1.7/1_07.sce
@@ -0,0 +1,19 @@
+pathname=get_absolute_file_path('1_07.sce')
+filename=pathname+filesep()+'1_07data.sci'
+exec(filename)
+E1=0.5*(Ea+Ec) + (((Ea-Eb)^2 +(Ec-Eb)^2)/2)^0.5;
+E2=0.5*(Ea+Ec) - (((Ea-Eb)^2 +(Ec-Eb)^2)/2)^0.5;
+disp(E1,"εI=",E2,"εII=");
+S1=E*(E1+ v*E2)/(1-v^2);//σI
+S2=E*(E2+ v*E1)/(1-v^2);//σII
+Sx=S1+S2;//σx
+printf("\nσI: %f N/mm^2",S1);
+printf("\nσII: %f N/mm^2",S2);
+printf("\nσx: %f N/mm^2",Sx);
+P=Sx*(%pi*d^2)/4;//axial tensile load
+Txy=(((S1-S2)^2 -(S1+S2)^2)/4)^0.5;//τxy
+printf("\nτxy: %f N/mm^2",Txy);
+J=(%pi*d^4)/32;//torsion constant
+T=2*Txy*J/d;//Torque
+disp(p/10^3,"P in KN",P,"P(N)=");
+disp(T/10^6,"T in KN.m",T,"T (N.mm)=");
+\ No newline at end of file
diff --git a/564/CH10/EX10.1/10_1.sce b/564/CH10/EX10.1/10_1.sce
new file mode 100755
index 000000000..e97916700
--- /dev/null
+++ b/564/CH10/EX10.1/10_1.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('10_1.sce')
+filename=pathname+filesep()+'10_1data.sci'
+exec(filename)
+lambda=(m*l^3)/(3*48*EI);//λ
+x=poly(0,'x');
+y=21*x^2 -22*x +1;
+m=roots(y);
+omega1=(m(1,:)/lambda)^0.5;//ω1
+omega2=(m(1,:)/lambda)^0.5;//ω2
+f1=omega1/(2*%pi);
+f2=omega2/(2*%pi);
+printf("\nf1:%f",f1);
+printf("\nf2:%f",f2);
+\ No newline at end of file
diff --git a/564/CH10/EX10.2/10_2.sce b/564/CH10/EX10.2/10_2.sce
new file mode 100755
index 000000000..13348740a
--- /dev/null
+++ b/564/CH10/EX10.2/10_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('10_2.sce')
+filename=pathname+filesep()+'10_2data.sci'
+exec(filename)
+lambda=(m*l^3)/(6*EI);//λ
+a=188,b=-44,c=1
+m1=(-b+(b*b -4*a*c)^0.5)/(2*a);
+m2=(-b-(b*b -4*a*c)^0.5)/(2*a);
+omega1=(m1/lambda)^0.5;//ω1
+omega2=(m2/lambda)^0.5;//ω2
+f1=[omega1;omega2];
+f=min(f1)/(2*%pi);
+printf("\nlowest natural frequency is: %f",f);
+\ No newline at end of file
diff --git a/564/CH10/EX10.4/10_4.sce b/564/CH10/EX10.4/10_4.sce
new file mode 100755
index 000000000..667240355
--- /dev/null
+++ b/564/CH10/EX10.4/10_4.sce
@@ -0,0 +1,9 @@
+pathname=get_absolute_file_path('10_4.sce')
+filename=pathname+filesep()+'10_4data.sci'
+exec(filename)
+for i=1:3
+    f(i)=(0.5*%pi*i/L)*(EI/(rho*A))^0.5;
+end
+printf("\nf1: %f",f(1));
+printf("\nf2: %f",f(2));
+printf("\nf3: %f",f(3));
+\ No newline at end of file
diff --git a/564/CH10/EX10.5/10_5.JPG b/564/CH10/EX10.5/10_5.JPG
new file mode 100755
index 000000000..024539a6d
--- /dev/null
+++ b/564/CH10/EX10.5/10_5.JPG
diff --git a/564/CH10/EX10.5/10_5.sce b/564/CH10/EX10.5/10_5.sce
new file mode 100755
index 000000000..d3e33d0e3
--- /dev/null
+++ b/564/CH10/EX10.5/10_5.sce
@@ -0,0 +1,31 @@
+pathname=get_absolute_file_path('10_5.sce')
+filename=pathname+filesep()+'10_5data.sci'
+exec(filename)
+for i=1:3
+    lambda(i)=(i-0.5)*%pi/L;
+    omega(i)=(lambda(i)^2)*(EI/(rho*A))^0.5;
+    fre(i)=omega(i)/(2*%pi);
+    K(i)=(cos(L*lambda(i))+cosh(L*lambda(i)))/(sin(L*lambda(i))+sinh(L*lambda(i)));
+end
+function[V1]=f(z)
+    V1=K(1)*(cosh(z*lambda(1)) -cos(z*lambda(1)) -K(1)*(sinh(z*lambda(1)) -sin(z*lambda(1))));
+endfunction
+function[V2]=f1(z)
+    V2=K(2)*(cosh(z*lambda(2)) -cos(z*lambda(2)) -K(2)*(sinh(z*lambda(2)) -sin(z*lambda(2))));
+endfunction
+function[V3]=f2(z)
+    V3=K(3)*(cosh(z*lambda(3)) -cos(z*lambda(3)) -K(3)*(sinh(z*lambda(3)) -sin(z*lambda(3))));
+endfunction
+z=linspace(0,L,10*L);
+subplot(3,1,1),xgrid(3),xtitle('1st Normal Mode Of Vibration','z','V')
+V1=feval(z,f);
+plot2d(z,V1);
+subplot(3,1,2),xgrid(3),xtitle('2nd Normal Mode Of Vibration','z','V')
+V2=feval(z,f1);
+plot2d(z,V2);
+subplot(3,1,3),xgrid(3),xtitle('3rd Normal Mode Of Vibration','z','V')
+V3=feval(z,f2);
+plot2d(z,V3);
+printf("\nf1: %f ",fre(1));
+printf("\nf2: %f ",fre(2));
+printf("\nf3: %f ",fre(3));
+\ No newline at end of file
diff --git a/564/CH10/EX10.6/10_6.sce b/564/CH10/EX10.6/10_6.sce
new file mode 100755
index 000000000..313569b49
--- /dev/null
+++ b/564/CH10/EX10.6/10_6.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('10_6.sce')
+filename=pathname+filesep()+'10_6data.sci'
+exec(filename)
+omega=1.1584*(EI/(m*l^3));
+printf("\nω1: %f",omega)
+\ No newline at end of file
diff --git a/564/CH12/EX12.1/12_1.sce b/564/CH12/EX12.1/12_1.sce
new file mode 100755
index 000000000..c92e80d1a
--- /dev/null
+++ b/564/CH12/EX12.1/12_1.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('12_1.sce')
+filename=pathname+filesep()+'12_1data.sci'
+exec(filename)
+clear
+Ar=(%pi*d^2)/4;
+b=((2*Ar*Ls)/(t*Lt))+d;
+printf("\nb: %f",b);
+\ No newline at end of file
diff --git a/564/CH12/EX12.2/12_2.sce b/564/CH12/EX12.2/12_2.sce
new file mode 100755
index 000000000..1123d6476
--- /dev/null
+++ b/564/CH12/EX12.2/12_2.sce
@@ -0,0 +1,16 @@
+pathname=get_absolute_file_path('12_2.sce')
+filename=pathname+filesep()+'12_2data.sci'
+exec(filename)
+clear
+Pe=Load*CL;
+V=Load/6;
+r=(DC^2 +AD^2)^0.5;
+sumr=4*r^2 +2*DC^2;//sum(r^2)
+S=(Pe/sumr)*r//shear on rivets A and B
+printf("\nshear force on A & B (S): %f N",S);
+thetaA=3*%pi/4;
+thetaB= %pi/4;
+RA=(S^2 + V^2 +2*S*V*cos(thetaA))^0.5;//resultent force on A
+RB=(S^2 + V^2 +2*S*V*cos(thetaB))^0.5;//resultent force on B
+printf("\nresultent force on A: %f N",RA);
+printf("\nresultent force on B: %f N",RB);
+\ No newline at end of file
diff --git a/564/CH14/EX14.1/14_1.sce b/564/CH14/EX14.1/14_1.sce
new file mode 100755
index 000000000..498d1472a
--- /dev/null
+++ b/564/CH14/EX14.1/14_1.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('14_1.sce')
+filename=pathname+filesep()+'14_1data.sci'
+exec(filename)
+T=((W/g)*a)/cos(10*%pi/180);
+printf("\nT: %f N",T);
+R=W+T*sin(10*%pi/180);
+Ls=(R/2)/cos(20*%pi/180);//Load in each strut
+printf("\nLoad in each strut: %f N",Ls);
+Li=(Wa/g)*a;//inertial load at CG of fuselage aft of AA
+N=T-Li*cos(10*%pi/180)+Wa*sin(10*%pi/180);
+S=Li*sin(10*%pi/180)+Wa*cos(10*%pi/180);
+s=v0^2/(2*a);
+printf("\nN: %f N",N);
+printf("\nS: %f N",S);
+printf("\ns: %f m",s);
+\ No newline at end of file
diff --git a/564/CH14/EX14.2/14_2.sce b/564/CH14/EX14.2/14_2.sce
new file mode 100755
index 000000000..cd8548a55
--- /dev/null
+++ b/564/CH14/EX14.2/14_2.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('14_2.sce')
+filename=pathname+filesep()+'14_2data.sci'
+exec(filename)
+clear
+ax=Rh/(W/g);//horizontal deceleration
+ay=(Rv-W)/(W/g);//vertical deceleration
+Ialpha=Rv*Sh +Rh*Sv;
+alpha=(Ialpha*10^6)/Icg;
+t=v0/ay;
+omega=alpha*t;
+printf("\nhorizontal reaction force: %f kN",W*ax/g);
+printf("\nvertical reaction force: %f kN",W*ay/g);
+printf("\nα: %f rad/s^2",alpha);
+printf("\nt: %f s",t);
+printf("\nω: %f rad/s",omega);
diff --git a/564/CH14/EX14.3/14_3.sce b/564/CH14/EX14.3/14_3.sce
new file mode 100755
index 000000000..6ab7bb516
--- /dev/null
+++ b/564/CH14/EX14.3/14_3.sce
@@ -0,0 +1,16 @@
+pathname=get_absolute_file_path('14_3.sce')
+filename=pathname+filesep()+'14_3data.sci'
+exec(filename)
+clear
+L=n*W;
+Cl1=(L/(0.5*S*rho*v^2));
+l1=4.18*cos((alpha1-2)*%pi/180) +0.31*sin((alpha1-2)*%pi/180);
+Cl2=(L/(0.5*S*rho*v^2))-((c/l1)*Cmcg1);
+l2=4.18*cos((alpha2-2)*%pi/180) +0.31*sin((alpha2-2)*%pi/180);
+Cl3=(L/(0.5*S*rho*v^2))-((c/l2)*Cmcg2);
+//CD from Fig. 14.8(a)
+CD=0.0875;
+Lift=0.5*rho*v*v*S*(Cl1+Cl2+Cl3)/3;
+Drag=0.5*rho*v*v*S*CD;
+printf("\nLift: %f N",Lift);
+printf("\nDrag: %f N",Drag);
+\ No newline at end of file
diff --git a/564/CH15/EX15.1/15_1.sce b/564/CH15/EX15.1/15_1.sce
new file mode 100755
index 000000000..0ae8d2d69
--- /dev/null
+++ b/564/CH15/EX15.1/15_1.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('15_1.sce')
+filename=pathname+filesep()+'15_1data.sci'
+exec(filename)
+clear
+af=(K^2)/(S^2 *alpha^2 *%pi);
+Nf=(2/(C*(n-2)*((S*%pi^0.5)^n)))*((1/ai^((n-2)/2))-(1/af^((n-2)/2)));
+printf("\naf: %f mm",af)
+if(round(Nf)>Nf) then
+    printf("\nNf: %f cycles",round(Nf));
+else
+    printf("\nNf: %f cycles",round(Nf)+1)
+end
diff --git a/564/CH16/EX16.1/16_1.JPG b/564/CH16/EX16.1/16_1.JPG
new file mode 100755
index 000000000..b8b208fbc
--- /dev/null
+++ b/564/CH16/EX16.1/16_1.JPG
diff --git a/564/CH16/EX16.1/16_1.sce b/564/CH16/EX16.1/16_1.sce
new file mode 100755
index 000000000..46ea290a9
--- /dev/null
+++ b/564/CH16/EX16.1/16_1.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('16_1.sce')
+filename=pathname+filesep()+'16_1data.sci'
+exec(filename)
+Ixx=(a*b^3)/12 - ((a-tx)*(b-2*ty)^3)/12;
+deff("[Sz]=f(y)","Sz=M*y/Ixx");
+y=[-b/2:0.05:b/2];
+fplot2d(y,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'Direct stress', ' -y- ', '-σz-');
+printf("\nσz at the top of the beam: %f N/mm^2",f(b/2));
+printf("\nσz at the bottom of the beam: %f N/mm^2",f(-b/2));
+\ No newline at end of file
diff --git a/564/CH16/EX16.10/16_10.JPG b/564/CH16/EX16.10/16_10.JPG
new file mode 100755
index 000000000..dbd630923
--- /dev/null
+++ b/564/CH16/EX16.10/16_10.JPG
diff --git a/564/CH16/EX16.10/16_10.sce b/564/CH16/EX16.10/16_10.sce
new file mode 100755
index 000000000..78d8af05b
--- /dev/null
+++ b/564/CH16/EX16.10/16_10.sce
@@ -0,0 +1,30 @@
+pathname=get_absolute_file_path('16_10.sce')
+filename=pathname+filesep()+'16_10data.sci'
+exec(filename)
+clear
+function[si]=sing(a,b)
+    if(a<b) then
+        si=0;
+    else
+        si=(a-b);
+    end
+endfunction
+deff("[v]=f(z)","v=(1/EI)*(((L*W*z^3)/64) -((W/24)*(sing(z,0.5*L))^4)+((W/24)*(sing(z,0.75*L))^4)-((27*W*z*L^3)/2048))");
+funcprot();
+for i=1:100*L
+    x0=f((i-1)/100),x1=f(i/100),x2=f((i+1)/100);
+    y1=(x0-x1),y2=(x1-x2);
+    if(y1*y2<0) then
+       P=i/100;
+       P1=f(i/100);
+       printf("\nmaximum deflection: %f m",P1);
+       printf("\nat z= %f m",P);
+    end
+end
+
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.11/16_11.JPG b/564/CH16/EX16.11/16_11.JPG
new file mode 100755
index 000000000..15b2f2adc
--- /dev/null
+++ b/564/CH16/EX16.11/16_11.JPG
diff --git a/564/CH16/EX16.11/16_11.sce b/564/CH16/EX16.11/16_11.sce
new file mode 100755
index 000000000..73ee297e6
--- /dev/null
+++ b/564/CH16/EX16.11/16_11.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('16_11.sce')
+filename=pathname+filesep()+'16_11data.sci'
+exec(filename)
+function[si]=sing(a,a1)
+    if(a<a1) then
+        si=0;
+    else
+        si=(a-a1);
+    end
+endfunction
+deff("[v]=f(z)","v=(M0/(EI*6*L))*(-z^3 +(3*L*(sing(z,b))^2)-(2*L^2 -6*L*b +3*b^2)*z)");
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.12/16_12.JPG b/564/CH16/EX16.12/16_12.JPG
new file mode 100755
index 000000000..ed1f629f1
--- /dev/null
+++ b/564/CH16/EX16.12/16_12.JPG
diff --git a/564/CH16/EX16.12/16_12.sce b/564/CH16/EX16.12/16_12.sce
new file mode 100755
index 000000000..631389db1
--- /dev/null
+++ b/564/CH16/EX16.12/16_12.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('16_12.sce')
+filename=pathname+filesep()+'16_12data.sci'
+exec(filename)
+P=(W*Ixy/(E*(Ixx*Iyy -Ixy^2)));
+P1=-(W*Iyy/(E*(Ixx*Iyy -Ixy^2)));
+deff("[u]=f(z)","u=P*(0.5*L*z^2 -(1/6)*z^3)");
+deff("[v]=f1(z)","v=P1*(0.5*L*z^2 -(1/6)*z^3)");
+funcprot();
+z=[0:0.05:L];
+subplot(2,1,1);
+fplot2d(z,f);
+xgrid(3);
+xtitle( 'deflection curve(u)', ' -z- ', '-u-');
+subplot(2,1,2);
+fplot2d(z,f1);
+xgrid(3);
+xtitle( 'deflection curve(v)', ' -z- ', '-v-');
+datatipToggle();
+printf("\nmaximum value of u: %f",f(L));
+printf("\nmaximum value of v: %f",f1(L));
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.13/16_13.sce b/564/CH16/EX16.13/16_13.sce
new file mode 100755
index 000000000..b104a3162
--- /dev/null
+++ b/564/CH16/EX16.13/16_13.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('16_13.sce')
+filename=pathname+filesep()+'16_13data.sci'
+exec(filename)
+dw=d- 2*tf;
+Ixx=2*(((b*tf^3)/12)+ b*tf*((dw+tf)/2)^2) + (tw*dw^3)/12;
+Iyy=(2*tf*b^3)/12 + (dw*dw^3)/12;
+printf("\n Ixx= %f mm^4",Ixx);
+printf("\n Iyy= %f mm^4",Iyy)
+\ No newline at end of file
diff --git a/564/CH16/EX16.14/16_14.JPG b/564/CH16/EX16.14/16_14.JPG
new file mode 100755
index 000000000..6820f4a0d
--- /dev/null
+++ b/564/CH16/EX16.14/16_14.JPG
diff --git a/564/CH16/EX16.14/16_14.sce b/564/CH16/EX16.14/16_14.sce
new file mode 100755
index 000000000..ebb6b5e94
--- /dev/null
+++ b/564/CH16/EX16.14/16_14.sce
@@ -0,0 +1,32 @@
+pathname=get_absolute_file_path('16_14.sce')
+filename=pathname+filesep()+'16_14data.sci'
+exec(filename)
+Ixx=(t*h^3)/3;
+Iyy=(t*h^3)/12;
+Ixy=(t*h^3)/8;
+P=Mx/(Ixx*Iyy -Ixy*Ixy);
+deff("[Sz1]=f(x)","Sz1=P*(Iyy*h*0.5 -Ixy*x)");//distribution of direct stress in top flange
+deff("[Sz2]=f1(y)","Sz2=P*(Iyy*y)");//distribution of direct stress in web
+deff("[Sz3]=f2(x1)","Sz3=P*(-Iyy*h*0.5 -Ixy*x1)");//distribution of direct stress in top flange
+deff("[Sz4]=f3(x)","Sz4=0*x");
+deff("[Sz5]=f4(y)","Sz5=0*y");
+y=[-h/2:0.05:h/2];
+x=[0:0.05:h/2];
+x1=[-h/2:0.05:0];
+funcprot(0);
+subplot(3,1,1);
+fplot2d(x,f);
+fplot2d(x,f3);
+xgrid(3);
+xtitle( 'Direct stress in upper flange', ' -x- ', '-σz-');
+subplot(3,1,2);
+fplot2d(y,f1);
+fplot2d(y,f4);
+xgrid(3);
+xtitle( 'Direct stress in web', ' -y- ', '-σz-');
+subplot(3,1,3);
+fplot2d(x1,f2);
+fplot2d(x1,f3);
+xgrid(3);
+xtitle( 'Direct stress in bottom flange', ' -y- ', '-σz-');
+datatipToggle();
diff --git a/564/CH16/EX16.15/16_15.sce b/564/CH16/EX16.15/16_15.sce
new file mode 100755
index 000000000..53b0fba91
--- /dev/null
+++ b/564/CH16/EX16.15/16_15.sce
@@ -0,0 +1,9 @@
+pathname=get_absolute_file_path('16_15.sce')
+filename=pathname+filesep()+'16_15data.sci'
+exec(filename)
+Nt=4*E*alpha*a*t*T0;
+Mxt=2*E*alpha*t*T0*a^2;
+Myt=-E*alpha*t*T0*a^2;
+printf("\nNt: %f ",Nt);
+printf("\nMxt: %f ",Mxt);
+printf("\nMyt: %f ",Myt);
+\ No newline at end of file
diff --git a/564/CH16/EX16.16/16_16.sce b/564/CH16/EX16.16/16_16.sce
new file mode 100755
index 000000000..cbad2db2d
--- /dev/null
+++ b/564/CH16/EX16.16/16_16.sce
@@ -0,0 +1,9 @@
+pathname=get_absolute_file_path('16_16.sce')
+filename=pathname+filesep()+'16_16data.sci'
+exec(filename)
+Nt=4*E*alpha*a*t*T0;
+Mxt=(8/3)*E*alpha*t*T0*a^2;
+Myt=-E*alpha*t*T0*a^2;
+printf("\nNt: %f ",Nt);
+printf("\nMxt: %f ",Mxt);
+printf("\nMyt: %f ",Myt);
+\ No newline at end of file
diff --git a/564/CH16/EX16.2/16_2.JPG b/564/CH16/EX16.2/16_2.JPG
new file mode 100755
index 000000000..d120afb4b
--- /dev/null
+++ b/564/CH16/EX16.2/16_2.JPG
diff --git a/564/CH16/EX16.2/16_2.sce b/564/CH16/EX16.2/16_2.sce
new file mode 100755
index 000000000..4cd42b9df
--- /dev/null
+++ b/564/CH16/EX16.2/16_2.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('16_2.sce')
+filename=pathname+filesep()+'16_2data.sci'
+exec(filename)
+clear
+Iyy=2*(ty*a^3)/12 +((b- 2*ty)*tx^3)/12;
+deff("[Sz]=f(x)","Sz=M*x/Iyy");
+x=[-a/2:0.05:a/2];
+fplot2d(x,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'Direct stress', ' -x- ', '-σ z-');
+printf("\nσz at the left hand edges of flanges: %f N/mm^2",f(-a/2));
+printf("\nσz at the right hand edges of flanges: %f N/mm^2",f(a/2));
+\ No newline at end of file
diff --git a/564/CH16/EX16.3/16_3.JPG b/564/CH16/EX16.3/16_3.JPG
new file mode 100755
index 000000000..fd647053c
--- /dev/null
+++ b/564/CH16/EX16.3/16_3.JPG
diff --git a/564/CH16/EX16.3/16_3.sce b/564/CH16/EX16.3/16_3.sce
new file mode 100755
index 000000000..3929fb6ae
--- /dev/null
+++ b/564/CH16/EX16.3/16_3.sce
@@ -0,0 +1,35 @@
+pathname=get_absolute_file_path('16_3.sce')
+filename=pathname+filesep()+'16_3data.sci'
+exec(filename)
+clear
+Iyy=2*(ty*a^3)/12 +((b- 2*ty)*tx^3)/12;
+Ixx=(a*b^3)/12 - ((a-tx)*(b-2*ty)^3)/12;
+Mx=M*cos(theta),My=M*sin(theta);
+alpha=(atan((My/Iyy)/(Mx/Ixx)))*180/%pi;
+deff("[Sz1]=f(x)","Sz1=((Mx/Ixx)*(b/2))-((My/Iyy)*x)");
+deff("[Sz2]=f1(x)","Sz2=((Mx/Ixx)*(-b/2))-((My/Iyy)*x)");
+deff("[Sz3]=f2(y)","Sz3=((Mx/Ixx)*y)");
+deff("[Sz4]=f3(x)","Sz4=0*x");
+deff("[Sz5]=f4(y)","Sz5=0*y");
+funcprot();
+y=[-b/2:0.05:b/2];
+x=[-a/2:0.05:a/2];
+funcprot(0);
+subplot(3,1,1);
+fplot2d(x,f);
+fplot2d(x,f3);
+xgrid(3);
+xtitle( 'Direct stress in upper flange', ' -x- ', '-σz-');
+subplot(3,1,3);
+fplot2d(x,f1);
+fplot2d(x,f3);
+xgrid(3);
+xtitle( 'Direct stress in lower flange', ' -x- ', '-σz-');
+subplot(3,1,2);
+fplot2d(y,f2);
+fplot2d(y,f4);
+xgrid(3);
+xtitle( 'Direct stress in web', ' -y- ', '-σz-');
+datatipToggle();
+printf("\nclick on the point on the plot to view its coordinates");
+
diff --git a/564/CH16/EX16.4/16_4.sce b/564/CH16/EX16.4/16_4.sce
new file mode 100755
index 000000000..42ffce120
--- /dev/null
+++ b/564/CH16/EX16.4/16_4.sce
@@ -0,0 +1,29 @@
+pathname=get_absolute_file_path('16_4.sce')
+ilename=pathname+filesep()+'16_4data.sci'
+exec(filename)
+Yc=(b^2 -t^2 +a*t)/(2*(a+b-t));
+Xc=((((a/2)-a1+ 0.5*t)*a) +((b-t)*t/2))/(a+b-t);
+Ixx=(1/3)*((t*((Yc-t)^3 -(Yc-b)^3))+(a*((Yc)^3 -(Yc-t)^3)));
+P=a1- 0.5*t +Xc;
+Iyy=(1/3)*((t*(P^3 -(P-a)^3))+((b-t)*(Xc^3- (Xc-t)^3)));
+Ixy=a*t^2 *(Yc- t*0.5) + (b-t)*t*(Yc+ 0.5*t)*12;
+M1=(1+round(100*Mx*Iyy/(Ixx*Iyy -Ixy^2)))/100,M2=(1+round(100*Mx*Ixy/(Ixx*Iyy -Ixy^2)))/100;
+function[z]=Sz(x,y)
+    z=M1*y -M2*x;
+endfunction
+Load=[Sz(-P,Yc);Sz(a-P,Yc);Sz(-P,Yc-t);Sz(a-P,Yc-t);Sz(-Xc,Yc-b);Sz(-Xc+t,Yc-b)];
+Point1=[-P;a-P;-P;a-P;-Xc;-Xc+t];
+Point2=[Yc;Yc;Yc-t;Yc-t;Yc-b;Yc-b];
+maximum=Load(1);
+for i=2:6
+    if(abs(Load(i))>abs(Load(i-1))) then
+        maximum=abs(Load(i));
+        mm=i;
+    end
+end
+printf("\nσz,max: %f N/mm^2",Load(mm));
+disp("at point");
+printf("\nX: %f ",Point1(mm));
+printf("\nY: %f ",Point2(mm));
+
+
diff --git a/564/CH16/EX16.5/16_5.JPG b/564/CH16/EX16.5/16_5.JPG
new file mode 100755
index 000000000..3a3a24591
--- /dev/null
+++ b/564/CH16/EX16.5/16_5.JPG
diff --git a/564/CH16/EX16.5/16_5.sce b/564/CH16/EX16.5/16_5.sce
new file mode 100755
index 000000000..fd43ab406
--- /dev/null
+++ b/564/CH16/EX16.5/16_5.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('16_5.sce')
+filename=pathname+filesep()+'16_5data.sci'
+exec(filename)
+deff("[v]=f(z)","v=(-W/(6*EI))*(3*L*z^2 -z^3)");
+funcprot();
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\ntip deflection: %f m",f(L));
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.6/16_6.JPG b/564/CH16/EX16.6/16_6.JPG
new file mode 100755
index 000000000..4b98b09c1
--- /dev/null
+++ b/564/CH16/EX16.6/16_6.JPG
diff --git a/564/CH16/EX16.6/16_6.sce b/564/CH16/EX16.6/16_6.sce
new file mode 100755
index 000000000..3c149b27b
--- /dev/null
+++ b/564/CH16/EX16.6/16_6.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('16_6.sce')
+filename=pathname+filesep()+'16_6data.sci'
+exec(filename)
+deff("[v]=f(z)","v=(-W/(24*EI))*(6*(L^2)*z^2 -4*L*z^3 +z^4)");
+funcprot();
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\ntip deflection: %f m",f(L));
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.7/16_7.JPG b/564/CH16/EX16.7/16_7.JPG
new file mode 100755
index 000000000..1c2e96ff6
--- /dev/null
+++ b/564/CH16/EX16.7/16_7.JPG
diff --git a/564/CH16/EX16.7/16_7.sce b/564/CH16/EX16.7/16_7.sce
new file mode 100755
index 000000000..3dd145e1f
--- /dev/null
+++ b/564/CH16/EX16.7/16_7.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('16_7.sce')
+filename=pathname+filesep()+'16_7data.sci'
+exec(filename)
+deff("[v]=f(z)","v=(W/(24*EI))*(2*L*z^3 -z*L^3 -z^4)");
+funcprot();
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\nmaximum deflection: %f m",f(L/2));
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.8/16_8.JPG b/564/CH16/EX16.8/16_8.JPG
new file mode 100755
index 000000000..5e95f34f4
--- /dev/null
+++ b/564/CH16/EX16.8/16_8.JPG
diff --git a/564/CH16/EX16.8/16_8.sce b/564/CH16/EX16.8/16_8.sce
new file mode 100755
index 000000000..84535f919
--- /dev/null
+++ b/564/CH16/EX16.8/16_8.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('16_8.sce')
+filename=pathname+filesep()+'16_8data.sci'
+exec(filename)
+deff("[v]=f(z)","v=(W/(48*EI))*(4*z^3 -3*z*L^2)");
+funcprot();
+z=[0:0.05:L];
+fplot2d(z,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'deflection curve', ' -z- ', '-v-');
+printf("\nmaximum deflection: %f m",f(L/2));
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH16/EX16.9/16_9.JPG b/564/CH16/EX16.9/16_9.JPG
new file mode 100755
index 000000000..355730a8e
--- /dev/null
+++ b/564/CH16/EX16.9/16_9.JPG
diff --git a/564/CH16/EX16.9/16_9.sce b/564/CH16/EX16.9/16_9.sce
new file mode 100755
index 000000000..241f1c205
--- /dev/null
+++ b/564/CH16/EX16.9/16_9.sce
@@ -0,0 +1,35 @@
+pathname=get_absolute_file_path('16_9.sce')
+filename=pathname+filesep()+'16_9data.sci'
+exec(filename)
+clear
+function[si]=sing(a,b)
+    if(a<b) then
+        si=0;
+    else
+        si=(a-b);
+    end
+endfunction
+deff("[v]=f(z)","v=(1/EI)*(((W*z^3)/8) -((W/6)*(sing(z,a))^3)-((W/6)*(sing(z,2*a))^3)+((W/3)*(sing(z,3*a))^3)-((5*W*z*a^2)/8))");
+for i=1:400*a
+    x0=f((i-1)/100),x1=f(i/100),x2=f((i+1)/100);
+    y1=(x0-x1),y2=(x1-x2);
+    if(y1*y2<0) then
+        if(y1<0) then
+            P=i/100;
+            P1=f(i/100);
+        else
+            Q=i/100;
+            Q1=f(i/100);
+        end
+    end
+end
+z=[0:0.05:4*a];
+fplot2d(z,f);
+xgrid(3);
+xtitle( 'deflection curve', ' -z- ', '-v-');
+datatipToggle();
+printf("\nmaximum positive deflection: %f m",P1);
+printf("\nat z= %f m",P);
+printf("\nmaximum negative deflection: %f m",Q1);
+printf("\nat z= %f m",Q);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH17/EX17.1/17_1.JPG b/564/CH17/EX17.1/17_1.JPG
new file mode 100755
index 000000000..43279e3d5
--- /dev/null
+++ b/564/CH17/EX17.1/17_1.JPG
diff --git a/564/CH17/EX17.1/17_1.sce b/564/CH17/EX17.1/17_1.sce
new file mode 100755
index 000000000..04583d5ad
--- /dev/null
+++ b/564/CH17/EX17.1/17_1.sce
@@ -0,0 +1,34 @@
+pathname=get_absolute_file_path('17_1.sce')
+filename=pathname+filesep()+'17_1data.sci'
+exec(filename)
+Ixx=(t*h^3)/3,Iyy=(t*h^3)/12,Ixy=(t*h^3)/8;
+A=0.5*(Sy/(Ixx*Iyy -Ixy^2));
+function[q1]=q12(s1),
+    q1=A*(((h*t*(Iyy-Ixy))*s1)+(Ixy*t*s1^2));
+endfunction
+s1=linspace(0,h/2,5*h);
+q1=feval(s1,q12);
+subplot(3,1,1)
+plot2d(s1,q1);
+xgrid(3);
+xtitle( 'Direct stress in lower flange', ' -s1- ', '-σz-');
+function[q2]=q23(s2),
+    q2=(A*((Iyy*t*h*s2)-(Iyy*t*s2^2)))+q12(h/2);
+endfunction
+s2=linspace(0,h,10*h);
+q2=feval(s2,q23);
+subplot(3,1,2)
+plot2d(s2,q2);
+xgrid(3);
+xtitle( 'Direct stress in the web', ' -s2- ', '-σz-');
+function[q3]=q34(s3),
+    q3=(A*((Ixy*t*s3*s3)-(Iyy*t*h*s3)))+q23(h);
+endfunction
+s3=linspace(0,h/2,5*h);
+q3=feval(s3,q34);
+subplot(3,1,3)
+plot2d(s3,q3);
+xgrid(3);
+xtitle( 'Direct stress in upper flange', ' -s3- ', '-σz-');
+datatipToggle();
+printf("\nclick on the point to view its coordinate on the plot")
+\ No newline at end of file
diff --git a/564/CH17/EX17.2/17_2.sce b/564/CH17/EX17.2/17_2.sce
new file mode 100755
index 000000000..d304277a0
--- /dev/null
+++ b/564/CH17/EX17.2/17_2.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('17_2.sce')
+filename=pathname+filesep()+'17_2data.sci'
+exec(filename)
+Es=(3*b*b)/(h*(1+ 6*b/h));
+printf("\nξs: %f mm",Es);
+\ No newline at end of file
diff --git a/564/CH17/EX17.3/17_3.sce b/564/CH17/EX17.3/17_3.sce
new file mode 100755
index 000000000..1a643be14
--- /dev/null
+++ b/564/CH17/EX17.3/17_3.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('17_3.sce')
+filename=pathname+filesep()+'17_3data.sci'
+exec(filename)
+Es=-3.35*a;
+printf("\nξs: %f mm",Es);
+\ No newline at end of file
diff --git a/564/CH18/EX18.1/18_1.sce b/564/CH18/EX18.1/18_1.sce
new file mode 100755
index 000000000..518fc8e23
--- /dev/null
+++ b/564/CH18/EX18.1/18_1.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('18_1.sce')
+filename=pathname+filesep()+'18_1data.sci'
+exec(filename)
+clear
+A=%pi*d*d/4;
+tmin1=Tmax/(2*A*Smax);
+tmin2= (Tmax*%pi*d*0.5*L)/(4*A*A*G*angle);
+if(tmin1<tmin2) then
+    printf("\nminimum allowable thickness is: %f mm",tmin2);
+else
+    printf("\nminimum allowable thickness is: %f mm",tmin1);
+end
+\ No newline at end of file
diff --git a/564/CH18/EX18.2/18_2.JPG b/564/CH18/EX18.2/18_2.JPG
new file mode 100755
index 000000000..7fb66ad17
--- /dev/null
+++ b/564/CH18/EX18.2/18_2.JPG
diff --git a/564/CH18/EX18.2/18_2.sce b/564/CH18/EX18.2/18_2.sce
new file mode 100755
index 000000000..360a186f9
--- /dev/null
+++ b/564/CH18/EX18.2/18_2.sce
@@ -0,0 +1,37 @@
+pathname=get_absolute_file_path('18_2.sce')
+filename=pathname+filesep()+'18_2data.sci'
+exec(filename)
+clear
+del=2*((b/tb)+(a/ta));//δ
+A=a*b;
+deff("[W41]=f(s1)","W41=(T/(2*A*G))*((s1/tb)+(del*a*s1/(4*A)))");
+deff("[W12]=f1(s2)","W12=(T/(2*A*G))*((s2/ta)+(del*b*s2/(4*A)))");
+deff("[W23]=f2(s1)","W23=-((T/(2*A*G))*((s1/tb)+(del*a*s1/(4*A))))");
+deff("[W34]=f3(s2)","W34=-((T/(2*A*G))*((s2/ta)+(del*b*s2/(4*A))))");
+deff("[Sz]=f4(s1)","Sz=0*s1");
+deff("[Sz1]=f5(s2)","Sz1=0*s2");
+funcprot(0);
+s1=[-b/2:0.05:b/2];
+s2=[-a/2:0.05:a/2];
+subplot(2,2,1)
+fplot2d(s1,f)
+fplot2d(s1,f4)
+xgrid(3);
+xtitle( 'Wraping in 4-1', ' -y- ', '-w-');
+subplot(2,2,2)
+fplot2d(s2,f1)
+fplot2d(s2,f5)
+xgrid(3);
+xtitle( 'Wraping in 1-2', ' -x- ', '-w-');
+subplot(2,2,3)
+fplot2d(s1,f2)
+fplot2d(s1,f4)
+xgrid(3);
+xtitle( 'Wraping in 2-3', ' -y- ', '-w-');
+subplot(2,2,4)
+fplot2d(s2,f3)
+fplot2d(s2,f5)
+xgrid(3);
+xtitle( 'Wraping in 3-4', ' -x- ', '-w-');
+datatipToggle();
+printf("\nclick on the point on the plot to view its coordinates")
diff --git a/564/CH18/EX18.3/18_3.JPG b/564/CH18/EX18.3/18_3.JPG
new file mode 100755
index 000000000..6dd469035
--- /dev/null
+++ b/564/CH18/EX18.3/18_3.JPG
diff --git a/564/CH18/EX18.3/18_3.sce b/564/CH18/EX18.3/18_3.sce
new file mode 100755
index 000000000..bacecbbdb
--- /dev/null
+++ b/564/CH18/EX18.3/18_3.sce
@@ -0,0 +1,37 @@
+pathname=get_absolute_file_path('18_3.sce')
+filename=pathname+filesep()+'18_3data.sci'
+exec(filename)
+clear
+J=(2*a*ty^3 + b*tx^3)/3;
+if(tx>ty) then
+    tmax=tx*T/J;
+else
+    tmax=ty*T/J;
+end
+printf("\nmaximum shear stress: %f N/mm^2",tmax);
+Ixx=a*ty*b*b/2 +(tx*b^3)/12;
+Es=(ty*(a*b)^2)/(4*Ixx);
+deff("[W32]=f(s1)","W32=-2*(T/(J*G))*(0.5*Es*s1)");
+deff("[W21]=f1(s2)","W21=-2*(T/(J*G))*(0.5*Es*a -0.5*a*s2)");
+deff("[W43]=f2(s2)","W43=2*(T/(J*G))*(0.5*Es*a -0.5*a*s2)");
+deff("[Sz]=f3(s1)","Sz=0*s1");
+deff("[Sz1]=f4(s2)","Sz1=0*s2");
+s1=[-b/2:0.05:b/2];
+s2=[0:0.05:a];
+subplot(3,1,1)
+fplot2d(s2,f1)
+fplot2d(s2,f4)
+xgrid(3);
+xtitle( 'Wraping in 2-1', ' -x- ', '-w-');
+subplot(3,1,2)
+fplot2d(s1,f)
+fplot2d(s1,f3)
+xgrid(3);
+xtitle( 'Wraping in 3-2', ' -y- ', '-w-');
+subplot(3,1,3)
+fplot2d(s2,f2)
+fplot2d(s2,f4)
+xgrid(3);
+xtitle( 'Wraping in 4-3', ' -x- ', '-w-');
+datatipToggle();
+printf("\nclick on the point on the plot to view its coordinates")
diff --git a/564/CH19/EX19.1/19_1.JPG b/564/CH19/EX19.1/19_1.JPG
new file mode 100755
index 000000000..3a83a603e
--- /dev/null
+++ b/564/CH19/EX19.1/19_1.JPG
diff --git a/564/CH19/EX19.1/19_1.sce b/564/CH19/EX19.1/19_1.sce
new file mode 100755
index 000000000..058b6e3d4
--- /dev/null
+++ b/564/CH19/EX19.1/19_1.sce
@@ -0,0 +1,54 @@
+pathname=get_absolute_file_path('19_1.sce')
+filename=pathname+filesep()+'19_1data.sci'
+exec(filename)
+clear
+Yc=(2*L45^2 +L12^2)/(4*L12 +4*L45);
+Ixx=((2*L12 +L45)*t*Yc^2) +(L45*t*(L45-Yc)^2)+ (2*t/3)*((2*Yc^3) -((Yc-L45)^3) -((Yc-L12)^3));
+P=-Sy/Ixx;
+function[q1]=q12(s1),
+    q1=P*t*(Yc-L12)*s1 +P*t*0.5*s1^2;
+endfunction
+s1=linspace(0,L12,10*L12);
+q1=feval(s1,q12);
+subplot(2,3,1)
+plot2d(s1,q1)
+xgrid(3);
+xtitle( 'Shear flow in 1-2 and 8-7', ' -s1- ', '-q-');
+function[q2]=q23(s2),
+    q2= P*t*Yc*s2+q12(L12);
+endfunction
+s2=linspace(0,L12,10*L12);
+q2=feval(s2,q23);
+subplot(2,3,2)
+plot2d(s2,q2)
+xgrid(3);
+xtitle( 'Shear flow in 2-3 and 7-6', ' -s2- ', '-q-');
+function[q3]=q03(s3),
+    q3=P*t*Yc*s3;
+endfunction
+s3=linspace(0,L45/2,5*L45);
+q3=feval(s3,q03);
+subplot(2,3,3)
+plot2d(s3,q3)
+xgrid(3);
+xtitle( 'Shear flow in 0-3 and 0-6', ' -s3- ', '-q-');
+function[q4]=q34(s4),
+    q4=(P*t*Yc*s4 -0.5*P*t*s4^2)+q03(L45/2)+q23(L12);
+endfunction
+s4=linspace(0,L45,10*L45);
+q4=feval(s4,q34);
+subplot(2,3,4)
+plot2d(s4,q4)
+xgrid(3);
+xtitle( 'Shear flow in 3-4 and 6-5', ' -s4- ', '-q-');
+function[q5]=q94(s5),
+    q5=P*t*(Yc-L45)*s5;
+endfunction
+s5=linspace(0,L45/2,5*L45);
+q5=feval(s5,q94);
+subplot(2,3,5)
+plot2d(s5,q5)
+xgrid(3);
+xtitle( 'Shear flow in 9-4 and 9-5', ' -s5- ', '-q-');
+datatipToggle();
+printf("\nclick on the point on the plot to view its coordinates");
+\ No newline at end of file
diff --git a/564/CH19/EX19.2/19_2.sce b/564/CH19/EX19.2/19_2.sce
new file mode 100755
index 000000000..1c1976c23
--- /dev/null
+++ b/564/CH19/EX19.2/19_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('19_2.sce')
+filename=pathname+filesep()+'19_2data.sci'
+exec(filename)
+GJcl=((4*A*A*G*t12)/(L12+L34));
+GJo=((L34+L13)*G*t^3)/3
+GJ=GJcl+GJo;
+Dtheta=T/GJ;
+qcl=GJcl*Dtheta/(2*A);
+tmaxcl=qcl/t12;
+tmaxo=G*t13*Dtheta;
+printf("\nτmax,cl: %f N/mm^2",tmaxcl);
+printf("\nτmax,op: %f N/mm^2",tmaxo);
+\ No newline at end of file
diff --git a/564/CH20/EX20.1/20_1.sce b/564/CH20/EX20.1/20_1.sce
new file mode 100755
index 000000000..d978ec44c
--- /dev/null
+++ b/564/CH20/EX20.1/20_1.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('20_1.sce')
+filename=pathname+filesep()+'20_1data.sci'
+exec(filename)
+L25=(L16*L12+L34*L23)/(L12+L23);
+B1=A+ (t16*L16/6) +(t12*L12/6)*(2+(L25/L16));
+B6=B1;
+B2= 2*A + (t12*L12/6)*(2+(L16/L25))+(t25*L25/6) +(t23*L23/6)*(2+(L34/L25));
+B5=B2;
+B3=A + (t23*L23/6)*(2+(L25/L34)) + (t34*L34/6);
+B4=B3;
+printf("\nB1 = B6 = %f mm^2",B1);
+printf("\nB2 = B5 = %f mm^2",B2);
+printf("\nB3 = B4 = %f mm^2",B3)
+\ No newline at end of file
diff --git a/564/CH20/EX20.2/20_2.sce b/564/CH20/EX20.2/20_2.sce
new file mode 100755
index 000000000..32495dff8
--- /dev/null
+++ b/564/CH20/EX20.2/20_2.sce
@@ -0,0 +1,14 @@
+pathname=get_absolute_file_path('20_2.sce')
+filename=pathname+filesep()+'20_2data.sci'
+exec(filename)
+clear
+Yc=(2*(B'*D) -B(1,:)*D(1,:))/((2*sum(B))-(B(1,:)+B(9,:)));
+Y=D-Yc*ones(9,1);
+for i=1:9
+    Ixx(i)=B(i)*(Y(i))^2;
+end
+P=2*sum(Ixx)-Ixx(1)-Ixx(9);
+for i=1:9
+    Sz(i)=(Mx/P)*Y(i);//σz
+end
+printf("\nσz: %f N/mm^2",Sz);
diff --git a/564/CH20/EX20.3/20_3.sce b/564/CH20/EX20.3/20_3.sce
new file mode 100755
index 000000000..4a1112c5a
--- /dev/null
+++ b/564/CH20/EX20.3/20_3.sce
@@ -0,0 +1,14 @@
+pathname=get_absolute_file_path('20_3.sce')
+filename=pathname+filesep()+'20_3data.sci'
+exec(filename)
+clear
+Ixx=4*A*L^2;
+B=[A;A;A;A];
+Y=[L;L;-L;-L];
+q(1)=(-Sy/Ixx)*B(1)*Y(1);
+for i=2:4
+    q(i)=((-Sy/Ixx)*B(i)*Y(i)) +q(i-1);
+end
+printf("\nq12: %f N/mm",q(1));
+printf("\nq23: %f N/mm",q(2));
+printf("\nq34: %f N/mm",q(3));
+\ No newline at end of file
diff --git a/564/CH20/EX20.4/20_4.sce b/564/CH20/EX20.4/20_4.sce
new file mode 100755
index 000000000..4b0de8a93
--- /dev/null
+++ b/564/CH20/EX20.4/20_4.sce
@@ -0,0 +1,22 @@
+pathname=get_absolute_file_path('20_4.sce')
+filename=pathname+filesep()+'20_4data.sci'
+exec(filename)
+clear
+for i=1:8
+    I(i)=B(i)*y(i)*y(i);
+end
+Ixx=sum(I);
+P=(-Sy/Ixx);
+qb1(1)=0,qb2(1)=0;
+for i=2:3
+    qb1(i)=P*B(i+1)*y(i+1) +qb1(i-1);
+    qb2(i)=(P*B(i+5)*y(i+5) +qb2(i-1));
+end
+qb1(4)=qb1(2);
+qb2(4)=qb2(2);
+qb=[qb1;-qb2];
+A= (y(1)+y(2))*x(3) +(y(2)+y(3))*x(2) +(y(3)+y(4))*x(1);
+qs0=(qb(7)*(x(2)+x(3))*(2*y(1)) +2*qb(6)*(x(3)*y(2) +x(3)*(y(2)-y(1))) + 2*qb(1)*x(2)*y(3) -2*qb(2)*x(1)*y(3) -qb(3)*2*y(4)*x(1))/(2*A);
+q=[qb1+qs0*ones(4,1);qb2+qs0*ones(4,1);];
+disp("shear flow(in order q23,q34,q45,q56,q67,q12,q78,q81)");
+printf("\n q: %f N/mm",q)
diff --git a/564/CH21/EX21.1/21_1.JPG b/564/CH21/EX21.1/21_1.JPG
new file mode 100755
index 000000000..5f3d1f416
--- /dev/null
+++ b/564/CH21/EX21.1/21_1.JPG
diff --git a/564/CH21/EX21.1/21_1.sce b/564/CH21/EX21.1/21_1.sce
new file mode 100755
index 000000000..7f66b1481
--- /dev/null
+++ b/564/CH21/EX21.1/21_1.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('21_1.sce')
+filename=pathname+filesep()+'21_1data.sci'
+exec(filename)
+Mx=-Sy*L1;//moment at section AA
+Wa=(We*(L-L1)+ W*L1)/L;//width of section AA
+Ixx=2*A*(Wa/2)^2 + (1/12)*t*Wa^3;
+Sz1=Mx*(Wa/2)/Ixx;
+Pz1=Sz1*A;
+Syz=Sy +Pz1*((W-We)/L);
+deff("[q12]=f(s)","q12=(-Syz/Ixx)*((-s^2 +Wa*s) +(A*Wa/2))");
+s=[0:0.01:Wa];
+fplot2d(s,f);
+xgrid(3)
+xtitle('shear flow','-s-','-q12-');
+datatipToggle();
+\ No newline at end of file
diff --git a/564/CH21/EX21.2/21_2.sce b/564/CH21/EX21.2/21_2.sce
new file mode 100755
index 000000000..e4bd28191
--- /dev/null
+++ b/564/CH21/EX21.2/21_2.sce
@@ -0,0 +1,28 @@
+pathname=get_absolute_file_path('21_2.sce')
+filename=pathname+filesep()+'21_2data.sci'
+exec(filename)
+clear
+Mx=-Sy*(L-Lc);//moment at section
+L3=(L1*(L-Lc)+ L2*Lc)/L;//leangth of CS
+B3=(B1*(L-Lc)+ B2*Lc)/L;//breadth of CS
+Ixx=4*Aco*(B3/2)^2 + 2*Ace*(B3/2)^2;
+B=[Aco;Ace;Aco;Aco;Ace;Aco];
+Y=[B3/2;B3/2;B3/2;-B3/2;-B3/2;-B3/2];
+Dxr=[(L1-L2)/(2*L);0;-(L1-L2)/(2*L);-(L1-L2)/(2*L);0;(L1-L2)/(2*L)];
+Dyr=[-(B1-B2)/(2*L);-(B1-B2)/(2*L);-(B1-B2)/(2*L);(B1-B2)/(2*L);(B1-B2)/(2*L);(B1-B2)/(2*L)];
+for i=1:6
+    Pz(i)=(Mx/Ixx)*Y(i)*B(i);
+    Px(i)=Pz(i)*Dxr(i);
+    Py(i)=Pz(i)*Dyr(i);
+    P(i)=(((Px(i))^2 +(Py(i))^2 +(Pz(i))^2)^0.5)*((abs(Pz(i)))/Pz(i));
+end
+Syw=Sy-sum(Py);
+qb(1)=0
+for i=2:4
+    qb(i)=qb(i-1)+(-Syw/Ixx)*B(i-1)*Y(i-1);
+end
+qb(5)= qb(3);
+qb(6)= qb(2);
+qs0=((-Sy*L3/2)+(qb(2)*L3*B3*0.5 +qb(3)*L3*B3*0.5 +qb(4)*B3*L3*0.5))/(2*L3*B3);
+q=-qb+qs0*ones(6,1);
+disp(q,"shear flow (61,12,23,34,45,56)")
+\ No newline at end of file
diff --git a/564/CH21/EX21.3/21_3.sce b/564/CH21/EX21.3/21_3.sce
new file mode 100755
index 000000000..e9db19df3
--- /dev/null
+++ b/564/CH21/EX21.3/21_3.sce
@@ -0,0 +1,40 @@
+pathname=get_absolute_file_path('21_3.sce')
+filename=pathname+filesep()+'21_3data.sci'
+exec(filename)
+clear
+Lc1=Lc+ 0.1*10^3,Lc2=Lc- 0.1*10^3;
+Mx1=-Sy*(L-Lc1);//moment at section1
+L31=(L1*(L-Lc1)+ L2*Lc1)/L;//leangth of CS1
+B31=(B1*(L-Lc1)+ B2*Lc1)/L;//breadth of CS1
+Mx2=-Sy*(L-Lc2);//moment at section2
+L32=(L1*(L-Lc2)+ L2*Lc2)/L;//leangth of CS2
+B32=(B1*(L-Lc2)+ B2*Lc2)/L;//breadth of CS2
+L3=(L31+L32)/2;
+B3=(B31+B32)/2;
+Ixx1=4*Aco*(B31/2)^2 + 2*Ace*(B31/2)^2;
+Ixx2=4*Aco*(B32/2)^2 + 2*Ace*(B32/2)^2;
+B=[Aco;Ace;Aco;Aco;Ace;Aco];
+Y1=[B31/2;B31/2;B31/2;-B31/2;-B31/2;-B31/2];
+Y2=[B32/2;B32/2;B32/2;-B32/2;-B32/2;-B32/2];
+Dxr=[(L1-L2)/(2*L);0;-(L1-L2)/(2*L);-(L1-L2)/(2*L);0;(L1-L2)/(2*L)];
+Dyr=[-(B1-B2)/(2*L);-(B1-B2)/(2*L);-(B1-B2)/(2*L);(B1-B2)/(2*L);(B1-B2)/(2*L);(B1-B2)/(2*L)];
+for i=1:6
+    Pz1(i)=(Mx1/Ixx1)*Y1(i)*B(i);
+    Px1(i)=Pz1(i)*Dxr(i);
+    Py1(i)=Pz1(i)*Dyr(i);
+    P1(i)=(((Px1(i))^2 +(Py1(i))^2 +(Pz1(i))^2)^0.5)*((abs(Pz1(i)))/Pz1(i));
+    Pz2(i)=(Mx2/Ixx2)*Y2(i)*B(i);
+    Px2(i)=Pz2(i)*Dxr(i);
+    Py2(i)=Pz2(i)*Dyr(i);
+    P2(i)=(((Px2(i))^2 +(Py2(i))^2 +(Pz2(i))^2)^0.5)*((abs(Pz2(i)))/Pz2(i));
+end
+delP=(P1-P2)/200;
+q12=(Sy*L3*0.5 +delP(2)*2*L3*B3*0.25 +(delP(2)+delP(3))*L3*B3*0.5 +delP(6)*L3*B3*0.5)/(2*L3*B3*0.25 +2*L3*B3*0.25 +L3*B3*0.5 +L3*B3*0.5);
+q23=q12-delP(2);
+q34=q12-(delP(2)+delP(3));
+q45=q12-(delP(2)+delP(3)+delP(4));
+q56=q12;
+q61=q12-delP(6);
+q=[q12;q23;q34;q45;q56;q61];
+disp("Shear flow(q12;q23;q34;q45;q56;q61):");
+printf("\n %f N/mm",q);
+\ No newline at end of file
diff --git a/564/CH22/EX22.1/22_1.sce b/564/CH22/EX22.1/22_1.sce
new file mode 100755
index 000000000..c68f9baad
--- /dev/null
+++ b/564/CH22/EX22.1/22_1.sce
@@ -0,0 +1,23 @@
+pathname=get_absolute_file_path('22_1.sce')
+filename=pathname+filesep()+'22_1data.sci'
+exec(filename)
+b=(2*%pi*r)/16;
+for i=1:16
+    angle(i)=(2*%pi/16)*(i-1);
+    y(i)=r*cos(angle(i));
+end
+B1=A+ 2*((t*b)/6)*(2 +y(2)/y(1));
+for i=1:16
+    B(i)=B1;
+    I(i)=B(i)*y(i)*y(i);
+end
+Ixx=sum(I);
+for i=1:16
+    Sz(i)=(Mx*y(i))/Ixx;
+    printf("\nσz%f =",i);
+    printf(" %f N/mm^2",Sz(i))
+end
+
+
+
+
diff --git a/564/CH22/EX22.2/22_2.sce b/564/CH22/EX22.2/22_2.sce
new file mode 100755
index 000000000..2f326bd0b
--- /dev/null
+++ b/564/CH22/EX22.2/22_2.sce
@@ -0,0 +1,33 @@
+pathname=get_absolute_file_path('22_2.sce')
+filename=pathname+filesep()+'22_2data.sci'
+exec(filename)
+b=(2*%pi*r)/16;
+for i=1:16
+    angle(i)=(2*%pi/16)*(i-1);
+    y(i)=r*cos(angle(i));
+end
+B1=A+ 2*((t*b)/6)*(2 +y(2)/y(1));
+for i=1:16
+    B(i)=B1;
+    I(i)=B(i)*y(i)*y(i);
+end
+Ixx=sum(I);
+qs(1)=0;
+qs(8)=qs(1);
+for i=2:5
+    qs(i)=((-Sy*B(i)*y(i))/Ixx)+qs(i-1);
+    qs(9-i)=qs(i);
+end
+qs1(1)=((-Sy*B(1)*y(1))/Ixx);
+qs1(8)=qs1(1);
+for i=2:5
+    qs1(i)=((-Sy*B(18-i)*y(18-i))/Ixx)+qs1(i-1);
+    qs1(9-i)=qs1(i);
+end
+Ac=%pi*r*r/16;
+qs0=-(((Sy*d)-2*Ac*(sum(qs1)-sum(qs)))/(2*Ac*16));
+Q=[qs;-qs1];
+for i=1:16
+    q(i)=Q(i)+qs0;
+    printf("\nq: %f N/mm",q(i))
+end
+\ No newline at end of file
diff --git a/564/CH23/EX23.1/23_1.sce b/564/CH23/EX23.1/23_1.sce
new file mode 100755
index 000000000..a02ead6f3
--- /dev/null
+++ b/564/CH23/EX23.1/23_1.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('23_1.sce')
+filename=pathname+filesep()+'23_1data.sci'
+exec(filename)
+clear
+for i=1:6
+    I(i)=B(i)*y(i)*y(i);
+end
+Ixx=sum(I);
+for i=1:6
+    S(i)=(Mx/Ixx)*y(i);
+end
+printf("\nσz: %f N/mm^2",S);
+\ No newline at end of file
diff --git a/564/CH23/EX23.2/23_2.sce b/564/CH23/EX23.2/23_2.sce
new file mode 100755
index 000000000..29699cc30
--- /dev/null
+++ b/564/CH23/EX23.2/23_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('23_2.sce')
+filename=pathname+filesep()+'23_2data.sci'
+exec(filename)
+clear
+for i=1:8
+    ts(i)=(G(i)*t(i))/Gref;
+    del(i)=L(i)/ts(i);
+end
+P=[(((del(1)+del(2))/A(1))+(del(2)/A(2))) -((del(2)/A(1))+((del(2)+del(3)+del(4)+del(5))/A(2))) (del(5)/A(2));
+   ((del(1)+del(2))/A(1))  ((del(5)/A(3))-(del(2)/A(1)))  (-(del(5)+del(6)+del(7)+del(8))/A(3));
+   2*A(1)  2*A(2)  2*A(3)];
+X=[0;0;T]
+q=inv(P)*X;
+disp("shear flows are");
+printf("\nqI: %f N/mm",q(1,:));
+printf("\nqII: %f N/mm",q(2,:));
+printf("\nqIII: %f N/mm\n",q(3,:));
+disp("shear stress distribution is(in order 12o,12i,13,24,34,35,46,56)");
+X1=[q(1,:)/t(1);(q(1,:)-q(2,:))/t(2);q(2,:)/t(3);q(2,:)/t(4);(q(2,:)-q(3,:))/t(5);
+    q(3,:)/t(6);q(3,:)/t(7);q(3,:)/t(8)]
+printf("\n %f N/mm^2",X1)
diff --git a/564/CH23/EX23.3/23_3.sce b/564/CH23/EX23.3/23_3.sce
new file mode 100755
index 000000000..5adf7c897
--- /dev/null
+++ b/564/CH23/EX23.3/23_3.sce
@@ -0,0 +1,29 @@
+pathname=get_absolute_file_path('23_3.sce')
+filename=pathname+filesep()+'23_3data.sci'
+exec(filename)
+clear
+for i=1:10
+    ts(i)=t(i)*G(i)/Gref;
+    del(i)=L(i)/ts(i);
+end
+for i=1:6
+    I(i)=B(i)*y(i)*y(i);
+end
+Ixx=sum(I);
+for i=1:6
+    q1(i)=(-Sy/Ixx)*B(i)*y(i);
+end
+q=[q1(2);q1(1);0;q1(5);q1(3);q1(4)]
+P=[(((del(4)+del(5)+del(6))/A(1))+(del(6)/A(2))) -(((del(3)+del(6)+del(10)+del(8))/A(2))+(del(6)/A(1))) del(8)/A(2);
+   ((del(4)+del(5)+del(6))/A(1)) -(((-del(8))/A(3))+(del(6)/A(1))) (-((del(2)+del(8)+del(7)+del(1)+del(9))/A(3)));
+   2*A(1) 2*A(2) 2*A(3)]
+X=[-((q(6)*del(5) +q(6)*del(6))/A(1))+((q(4)*del(8) +q(5)*del(6))/A(2));
+   -((q(6)*del(5) +q(6)*del(6))/A(1))+((-q(2)*del(9) +q(1)*del(8)+q(1)*del(7))/A(3));
+   -q(5)*L(10)*L(5)-q(5)*L(10)*L(6)+q(2)*L(1)*L(9)];
+X1=inv(P)*X;
+X2=[X1(1,:);X1(2,:);X1(2,:);X1(3,:);X1(3,:);-q(2)+X1(3,:);q(4)-X1(3,:);q(4)-X1(1,:);q(6)+X1(1,:);-q(5)-X1(1,:)];
+dth=(1/(2*Gref))*((((del(4)+del(5)+del(6))/A(1))*X1(1,:))-(del(6)/A(1))*X1(2,:)+((q(6)*del(5) +q(6)*del(6))/A(1)));
+printf("\ndθ/dz: %f rad/mm\n",dth);
+disp("shear flow distribution is(in order q34,q23,q87,q12,q56,q61,q57,q72,q48,q83)");
+printf("\n %f N/mm",X2)
+
diff --git a/564/CH23/EX23.4/23_4.sce b/564/CH23/EX23.4/23_4.sce
new file mode 100755
index 000000000..52d14df49
--- /dev/null
+++ b/564/CH23/EX23.4/23_4.sce
@@ -0,0 +1,45 @@
+pathname=get_absolute_file_path('23_4.sce')
+filename=pathname+filesep()+'23_4data.sci'
+exec(filename)
+clear
+y1=0.5*[W1;W1;W1;-W1;-W1;-W1];
+y2=0.5*[W2;W2;W2;-W2;-W2;-W2];
+x1=[0;L12;L12+L11;L12+L11;L12;0];
+x2=[0;L22;L22+L21;L22+L21;L22;0];
+Li=[L12;L11;W1;L11;L12;W1;W1];
+for i=1:6;
+    I1(i)=B(i)*y1(i)*y1(i);
+    del(i)=Li(i)/t(3);
+end
+del(7)=Li(7)/t(3);
+Ixx=sum(I1);
+dely=(y1-y2);
+delx=(x1-x2);
+epr=[L12;0;L11;L11;0;L12];
+nr=abs(y1);
+for i=1:6
+    Pz(i)=(Mx/Ixx)*B(i)*y1(i);
+    Py(i)=Pz(i)*dely(i)/L;
+    Px(i)=Pz(i)*delx(i)/L;
+    Pr(i)=((Px(i)^2 +Py(i)^2 +Pz(i)^2)^0.5)*(y1(i)/abs(y1(i)));
+    Pxn(i)=-abs(Px(i)*nr(i));
+    Pyep(i)=Py(i)*epr(i); 
+end
+Pyep(6)=-Pyep(i);
+Sxw=-sum(Px);
+Syw=Sy-sum(Py);
+qb=[0;0;(-Syw/Ixx)*(B(3)*y1(3));0;0;(-Syw/Ixx)*(B(6)*y1(6));(-Syw/Ixx)*B(5)*y1(5)];
+qb1=0;
+for i=1:7
+    qb1=qb(i)*del(i) +qb1;
+end
+A1=L11*W1;
+A2=L12*W1;
+P=[((del(2)+del(3)+del(4)+del(7))/A1)+(del(7))/A2   -((del(7)/A1)+((del(1)+del(5)+del(6)+del(7))/A2));
+   2*A1  2*A2];
+X=[(qb1/(3*A1))+(qb1/(3*A2));-(qb(3)*W1*L12 +qb(6)*W1*L11)-sum(Pxn)+sum(Pyep)];
+qs=inv(P)*X;
+M1=-[qs(2);qs(1);qs(1);qs(1);qs(2);qs(2);qs(1)-qs(2)];
+q=qb+M1;
+disp("shear flow(q12;q23;q34;q45;q56;q61;q52):");
+printf("\n %f N/mm",q)
+\ No newline at end of file
diff --git a/564/CH23/EX23.6/23_6.JPG b/564/CH23/EX23.6/23_6.JPG
new file mode 100755
index 000000000..a151a5d1a
--- /dev/null
+++ b/564/CH23/EX23.6/23_6.JPG
diff --git a/564/CH23/EX23.6/23_6.sce b/564/CH23/EX23.6/23_6.sce
new file mode 100755
index 000000000..7f61e1ef8
--- /dev/null
+++ b/564/CH23/EX23.6/23_6.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('23_6.sce')
+filename=pathname+filesep()+'23_6data.sci'
+exec(filename)
+clear
+A=L1*B1;
+q1=T/(2*A);
+S=T/B1;
+q1=S/L1;
+P=S*(L/2)/L1;
+X=[L -L;A A];
+X1=[P;T];
+q=[q1;inv(X)*X1];
+X2=[0;L*q(2)-L*q(3);-L*q(2)+L*q(3);0];
+Y=[0;1500;3000;4500];
+plot(Y,X2);
+xgrid(3);
+xtitle('Distribution of load in the top flange of the front spar',' ','Load')
+datatipToggle();
+printf("\nq1: %f N/mm",q(1))
+printf("\nq2: %f N/mm",q(2))
+printf("\nq3: %f N/mm",q(3))
+\ No newline at end of file
diff --git a/564/CH24/EX24.1/24_1.JPG b/564/CH24/EX24.1/24_1.JPG
new file mode 100755
index 000000000..5dedf43d1
--- /dev/null
+++ b/564/CH24/EX24.1/24_1.JPG
diff --git a/564/CH24/EX24.1/24_1.sce b/564/CH24/EX24.1/24_1.sce
new file mode 100755
index 000000000..81c58a7c5
--- /dev/null
+++ b/564/CH24/EX24.1/24_1.sce
@@ -0,0 +1,29 @@
+pathname=get_absolute_file_path('24_1.sce')
+filename=pathname+filesep()+'24_1data.sci'
+exec(filename)
+clear
+P=[GH -GH;DK KH],X=[L1*sin(theta);L1*cos(theta)];
+q=inv(P)*X;
+q(3)=L1*cos(theta)/(DK+KH);
+q(4)=(L1*cos(theta)+L2)/(DK+KH);
+PA=GH*q(1) +FG*q(3) +EF*q(4);
+PE=-GH*q(2) -FG*q(3) -EF*q(4);
+X1=[GH*q(1)+FG*q(3)+EF*q(4);GH*q(1)+FG*q(3);GH*q(1);0];
+X2=[-GH*q(2)-FG*q(3)-EF*q(4);-GH*q(2)-FG*q(3);-GH*q(2);0];
+Y=[0;EF;EF+FG;EF+FG+GH];
+subplot(2,1,1)
+plot2d(Y,X1)
+xgrid(3)
+xtitle('Load distribution in ABCD','y','Load')
+subplot(2,1,2)
+plot2d(Y,X2)
+xgrid(3)
+xtitle('Load distribution in EFGH','y','Load')
+datatipToggle();
+printf("\nq1: %f N/mm",q(1));
+printf("\nq2: %f N/mm",q(2));
+printf("\nq3: %f N/mm",q(3));
+printf("\nq4: %f N/mm",q(4));
+printf("\nPA: %f N",PA);
+printf("\nPE: %f N",PE);
+printf("\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH24/EX24.2/24_2.sce b/564/CH24/EX24.2/24_2.sce
new file mode 100755
index 000000000..9282c6a7f
--- /dev/null
+++ b/564/CH24/EX24.2/24_2.sce
@@ -0,0 +1,30 @@
+pathname=get_absolute_file_path('24_2.sce')
+filename=pathname+filesep()+'24_2data.sci'
+exec(filename)
+clear
+P=[2*L11 -2*L11 0;0 -L11 L11;2*A2 2*(A1+A2) 0],X=[L1;L2;440000];
+q=inv(P)*X;
+//actual X is X=[L1;L2;-L2*L11],it leads to wrong answers than book;
+Sy1=-(q(2))*L11;
+Px4= 2*A1*-(q(2))/L11;
+Py2= Px4*tan(theta);
+q1= (Sy1- 2*Py2)/L11;
+P2= (Px4*Px4 +Py2*Py2)^0.5;
+P5= 2*((A1+A22)*(-q(2)) - A21*q(1))/L12;
+q2= ((-q(2))*L11 +(-q(2))*0.5*(L12-L11) -q(1)*0.5*(L12-L11))/L12;
+q31= (q2*L12 +L2)/L12;
+M3= 2*((A1+A2)*(-q(2)) -A2*q(1))+ L2*L11;
+Px1=M3/L11;
+Py1=Px1*tan(theta);
+P1=(Px1*Px1+ Py1*Py1)^0.5;
+q32=((L2+ L11*(-q(2))) -(2*Py1))/L11;
+printf("\nq12: %f N/mm",q(1));
+printf("\nq23: %f N/mm",q(2));
+printf("\nq31: %f N/mm",q(3));
+printf("\nP2= P4= %f N",P2);
+printf("\nP5= P6= %f N",P5);
+printf("\nP1= P3= %f N",P1);
+printf("\nq1: %f N/mm",q1);
+printf("\nq2: %f N/mm",q2);
+printf("\nq3: %f N/mm",q31);
+printf("\nshear flow q3 in the web: %f N/mm",q32)
+\ No newline at end of file
diff --git a/564/CH25/EX25.1/25_1.sce b/564/CH25/EX25.1/25_1.sce
new file mode 100755
index 000000000..d73542279
--- /dev/null
+++ b/564/CH25/EX25.1/25_1.sce
@@ -0,0 +1,18 @@
+pathname=get_absolute_file_path('25_1.sce')
+filename=pathname+filesep()+'25_1data.sci'
+exec(filename)
+clear
+A=(L1*(Bc+Be1+Be2));
+E1= ((Ec*L1*Bc)+(Ee*L1*(Be1+Be2)))/A;
+S1=Load/A;
+e1=S1/E1;
+dell=e1*L;
+v1= ((vc*L1*Bc)+(ve*L1*(Be1+Be2)))/A;
+et=-e1*v1;
+delt=-et*(Bc+Be1+Be2);
+Se=e1*Ee;
+Sc=e1*Ec;
+printf("\nlengthening of the bar: %f mm",dell);
+printf("\nreduction in thickness: %f mm",delt);
+printf("\nσm(epoxy):  %f N/mm^2",Se);
+printf("\nσf(carbon): %f N/mm^2",Sc);
+\ No newline at end of file
diff --git a/564/CH25/EX25.2/25_2.sce b/564/CH25/EX25.2/25_2.sce
new file mode 100755
index 000000000..914107e35
--- /dev/null
+++ b/564/CH25/EX25.2/25_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('25_2.sce')
+filename=pathname+filesep()+'25_2data.sci'
+exec(filename)
+clear
+vtl=vlt*(Et/El);
+e1=(S1/El)-(vtl*S2/Et);
+e2=(S2/Et)-(vlt*S1/El);
+gammalt=S3/Glt;//γlt
+printf("\nνtl:%f ",vtl);
+printf("\nεl:%f ",e1);
+printf("\nεt:%f ",e2);
+printf("\nγlt:%f ",gammalt);
+\ No newline at end of file
diff --git a/564/CH25/EX25.3/25_3.sce b/564/CH25/EX25.3/25_3.sce
new file mode 100755
index 000000000..d857effd9
--- /dev/null
+++ b/564/CH25/EX25.3/25_3.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('25_3.sce')
+filename=pathname+filesep()+'25_3data.sci'
+exec(filename)
+clear
+s11=1/El;
+s22=1/Et;
+s12=-vlt/El;
+s33=1/Glt;
+m=cos(theta),n=sin(theta);
+P=[(s11*m^4)+(s22*n^4)+(2*s12*m*m*n*n)+(s33*m*m*n*n) (m*m*n*n)*(s11+s22-s33)+(m^4 +n^4)*s12 0;
+   (m*m*n*n)*(s11+s22-s33)+(m^4 +n^4)*s12 (s11*n^4)+(s22*m^4)+(2*s12*m*m*n*n)+(s33*m*m*n*n) 0;
+   ((-n*m^3 +m*n^3)*(2*s12 +s33))-(2*s22*m*n^3)+(2*n*s11*m^3) ((n*m^3 -m*n^3)*(2*s12 +s33))+(2*s11*m*n^3)-(2*n*s22*m^3) 0];
+X=[S1;S2;0];
+E=P*X;
+printf("\nεx:  %f ",E(1));
+printf("\nεy: %f ",E(2));
+printf("\nγxy:%f ",E(3));
+\ No newline at end of file
diff --git a/564/CH25/EX25.4/25_4.sce b/564/CH25/EX25.4/25_4.sce
new file mode 100755
index 000000000..9a64d1f9b
--- /dev/null
+++ b/564/CH25/EX25.4/25_4.sce
@@ -0,0 +1,11 @@
+pathname=get_absolute_file_path('25_4.sce')
+filename=pathname+filesep()+'25_4data.sci'
+exec(filename)
+clear
+X1=a*tb*Ez1;
+X2=b*ta*Ez2;
+ez=Load/(2*X1+X2);
+P1=ez*X1;
+P2=ez*X2;
+printf("\nP(flanges) %f N",P1);
+printf("\nP(web) %f N",P2);
+\ No newline at end of file
diff --git a/564/CH25/EX25.5/25_5.sce b/564/CH25/EX25.5/25_5.sce
new file mode 100755
index 000000000..cc3a5d251
--- /dev/null
+++ b/564/CH25/EX25.5/25_5.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('25_5.sce')
+filename=pathname+filesep()+'25_5data.sci'
+exec(filename)
+clear
+Ixx=(2*Ez1*a*tb*(b/2)^2) +(Ez2*ta*b^3)/12;
+Iyy=(Ez1*tb*(2*a)^3)/12;
+Ixy= Ez1*a*tb*a*(b/2) +Ez1*a*tb*(-a)*(-b/2);
+P1=(-Mx*Ixy)/(Ixx*Iyy-Ixy*Ixy);
+P2=(Mx*Iyy)/(Ixx*Iyy-Ixy*Ixy);
+function[S1]=Sz1(x,y),//stress in flanges
+    S1=Ez1*(P1*x+P2*y),
+endfunction
+function[S2]=Sz2(x,y),//stress in web
+    S2=Ez2*(P1*x+P2*y),
+endfunction
+X=[Sz1(a,0.5*b);Sz1(0,0.5*b);Sz1(0,-0.5*b);Sz1(-a,-0.5*b);Sz2(0,0.5*b);Sz2(0,-0.5*b)];
+printf("\nmaximum direct stress in the beam cross-section is: %f N/mm^2",max(X));
+\ No newline at end of file
diff --git a/564/CH25/EX25.6/25_6.JPG b/564/CH25/EX25.6/25_6.JPG
new file mode 100755
index 000000000..b6afbbfc9
--- /dev/null
+++ b/564/CH25/EX25.6/25_6.JPG
diff --git a/564/CH25/EX25.6/25_6.sce b/564/CH25/EX25.6/25_6.sce
new file mode 100755
index 000000000..4db2197fe
--- /dev/null
+++ b/564/CH25/EX25.6/25_6.sce
@@ -0,0 +1,34 @@
+pathname=get_absolute_file_path('25_6.sce')
+filename=pathname+filesep()+'25_6data.sci'
+exec(filename)
+clear
+Ixx=((2*E12*t12*((L23/(2*L12))^2)*L12^3)+(E23*t23*L23^3))/12;
+alpha=asin(L23*0.5/L12);
+function[qb12]=qb12x(s1)
+    qb12= 0.5*(t12*sin(alpha)*E12*Sy/Ixx)*s1^2;
+endfunction
+function[qb23]=qb23x(s2)
+    qb23=((-E23*t23*Sy/Ixx)*(-0.5*L23*s2 +0.5*s2^2))+qb12x(L12);
+endfunction
+funcprot();
+qs0=((Sy*L12)+(L23*L12*qb12x(L12)/3))/(L12*L23);
+function[q12]=q12x(s1)
+    q12=qb12x(s1)-qs0;
+endfunction
+function[q23]=q23x(s2)
+    q23=qb23x(s2)-qs0;
+endfunction
+s1=linspace(0,L12,10*L12);
+s2=linspace(0,L23,10*L23);
+q12=feval(s1,q12x);
+q23=feval(s2,q23x);
+subplot(2,1,1)
+plot2d(s1,q12);
+xgrid(3)
+xtitle('','s1','q12')
+subplot(2,1,2)
+plot2d(s2,q23);
+xgrid(3)
+xtitle('','s2','q23')
+datatipToggle();
+printf("\nclick on the point to view its coordinate on plot")
+\ No newline at end of file
diff --git a/564/CH25/EX25.7/25_7.JPG b/564/CH25/EX25.7/25_7.JPG
new file mode 100755
index 000000000..0f0fad6aa
--- /dev/null
+++ b/564/CH25/EX25.7/25_7.JPG
diff --git a/564/CH25/EX25.7/25_7.sce b/564/CH25/EX25.7/25_7.sce
new file mode 100755
index 000000000..72cd8f699
--- /dev/null
+++ b/564/CH25/EX25.7/25_7.sce
@@ -0,0 +1,23 @@
+pathname=get_absolute_file_path('25_7.sce')
+filename=pathname+filesep()+'25_7data.sci'
+exec(filename)
+clear
+A=l*b;
+q=T/(2*A);
+P=(2*b/(t1*Gc))+(2*l/(t2*Gw));
+w=q*((0.5*l/(Gw*t2))- 0.25*P);
+X1=0.5*[l,-l],Y1=[w,-w];
+X2=0.5*[l,-l],Y2=[-w,w];
+X3=0.5*[b,-b],Y3=[w,-w];
+X4=0.5*[b,-b],Y4=[-w,w];
+subplot(2,2,1)
+plot(X1,Y1),xgrid(3),xtitle('warpinr in 2-3(mm)','y','w')
+subplot(2,2,2)
+plot(X2,Y2),xgrid(3),xtitle('warpinr in 1-4(mm)','y','w')
+subplot(2,2,3)
+plot(X3,Y3),xgrid(3),xtitle('warpinr in 2-1(mm)','x','w')
+subplot(2,2,4)
+plot(X4,Y4),xgrid(3),xtitle('warpinr in 3-4(mm)','x','w')
+datatipToggle();
+printf("\nW1: %f mm",w);
+printf("\nclick on the point to view its coordinate on plot")
diff --git a/564/CH25/EX25.8/25_8.JPG b/564/CH25/EX25.8/25_8.JPG
new file mode 100755
index 000000000..124316bb7
--- /dev/null
+++ b/564/CH25/EX25.8/25_8.JPG
diff --git a/564/CH25/EX25.8/25_8.sce b/564/CH25/EX25.8/25_8.sce
new file mode 100755
index 000000000..4677f1f10
--- /dev/null
+++ b/564/CH25/EX25.8/25_8.sce
@@ -0,0 +1,20 @@
+pathname=get_absolute_file_path('25_8.sce')
+filename=pathname+filesep()+'25_8data.sci'
+exec(filename)
+clear
+GJ=2*G1*a*(t1^3)/3 +G2*b*(t2^3)/3;
+dtheta=T/GJ;
+tmax12=2*G1*(t1/2)*dtheta;
+tmax23=2*G2*(t2/2)*dtheta;
+maximum=tmax23;
+if(tmax12>tmax23) then
+    maximum=tmax12;
+end
+printf("\nmaximum shear stress: %f N/mm^2",maximum);
+w1=-2*(0.5*a*b/2)*dtheta;
+printf("\nw1: %f mm",w1);
+X=[0;a],Y=[0,w1],X1=-X,Y1=-Y;
+subplot(2,1,1),plot(X,Y),xgrid(3),xtitle('warping in upper flange','x','w');
+subplot(2,1,2),plot(X1,Y1),xgrid(3),xtitle('warping in lower flange','x','w');
+datatipToggle();
+printf("\nclick on the point to view its coordinate on plot");
+\ No newline at end of file
diff --git a/564/CH26/EX26.1/26_1.sce b/564/CH26/EX26.1/26_1.sce
new file mode 100755
index 000000000..5975e6471
--- /dev/null
+++ b/564/CH26/EX26.1/26_1.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_1.sce')
+filename=pathname+filesep()+'26_1data.sci'
+exec(filename)
+clear
+L41=L23*cos(asin((L12-L34)/L23));
+t41=t23;
+C=[0;cos(asin((L12-L34)/L23));0;1];
+S=[1;sin(asin((L12-L34)/L23));1;0];
+P=[1 0 0 0 0 -t12*G 0;
+   0 1 0 0 t23*G*cos(asin((L12-L34)/L23)) t23*G*sin(asin((L12-L34)/L23)) -t23*G*L12*cos(asin((L12-L34)/L23));
+   0 0 1 0 0 t34*G -t34*G*L23*cos(asin((L12-L34)/L23));
+   0 0 0 1 -t41*G 0 0;
+   0 1 0 -1 0 0 0;
+   L12 -L23*sin(asin((L12-L34)/L23)) -L34 0 0 0 0;
+   0 L23*L12*cos(asin((L12-L34)/L23)) L34*L23*cos(asin((L12-L34)/L23)) 0 0 0 0];
+X=[0;0;0;0;0;Load;Load*x];
+q=inv(P)*X;
+printf("\nτ12: %f N/mm^2",q(1)/t12);
+printf("\nτ23: %f N/mm^2",q(2)/t23);
+printf("\nτ34: %f N/mm^2",q(3)/t34);
+printf("\nτ41: %f N/mm^2",q(4)/t41);
+\ No newline at end of file
diff --git a/564/CH26/EX26.2/26_2.JPG b/564/CH26/EX26.2/26_2.JPG
new file mode 100755
index 000000000..6d7566002
--- /dev/null
+++ b/564/CH26/EX26.2/26_2.JPG
diff --git a/564/CH26/EX26.2/26_2.sce b/564/CH26/EX26.2/26_2.sce
new file mode 100755
index 000000000..d3955550a
--- /dev/null
+++ b/564/CH26/EX26.2/26_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_2.sce')
+filename=pathname+filesep()+'26_2data.sci'
+exec(filename)
+clear
+du=0;
+P=[1 0 0 G*t*R;
+   0 1 -G*t 2*R*G*t;
+   0 0 1.79 -R
+   0 0 -0.13 R];
+X=[0;0;Load/(2*R*G*t);Load/((-0.5*%pi -6)*G*R*t)];
+q=inv(P)*X;
+printf("\nq12 = q34 = %f N/mm",q(1));
+printf("\nq41= %f N/mm",q(2));
+function[q23]=f(phi)
+    q23=-R*G*t*q(4) -G*t*sin(phi)*q(3);
+endfunction
+phi=linspace(0,%pi,%pi*10);
+q23=feval(phi,f);
+plot2d(phi,q23);
+xgrid(3),xtitle('q23 versus φ','φ','q23');
+datatipToggle();
+\ No newline at end of file
diff --git a/564/CH26/EX26.3/26_3.sce b/564/CH26/EX26.3/26_3.sce
new file mode 100755
index 000000000..7d8f7d373
--- /dev/null
+++ b/564/CH26/EX26.3/26_3.sce
@@ -0,0 +1,19 @@
+pathname=get_absolute_file_path('26_3.sce')
+filename=pathname+filesep()+'26_3data.sci'
+exec(filename)
+clear
+Tr=T*L1/2;
+mu=(8*G*t/(A*E*(b+a)))^0.5;
+L=L1/2;
+k1=((T*(b-a)*10^3)/(8*a*b*G*t));
+k2=1/(mu*cosh(mu*L));
+k3=((4*(b-a))/(a*b*(b+a)));
+k4=(2*T*(10^3)/(a*b*G*t*(b+a)));
+function[th]=f(z)
+    w=(k1*((k2*sinh(mu*z))-z));
+    F=((k1*k3 +k4)*L*L*0.5 -(k1*k2*k3/mu)*cosh(mu*L));
+    th=(k1*k2*k3/mu)*cosh(mu*z) -(k1*k3 +k4)*z*z*0.5 + F;
+endfunction
+funcprot();
+printf("\nangle of twist at mid-span θ: %f rad",f(0));
+
diff --git a/564/CH26/EX26.4/26_4.JPG b/564/CH26/EX26.4/26_4.JPG
new file mode 100755
index 000000000..9babb9f38
--- /dev/null
+++ b/564/CH26/EX26.4/26_4.JPG
diff --git a/564/CH26/EX26.4/26_4.sce b/564/CH26/EX26.4/26_4.sce
new file mode 100755
index 000000000..e0a298ead
--- /dev/null
+++ b/564/CH26/EX26.4/26_4.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_4.sce')
+filename=pathname+filesep()+'26_4data.sci'
+exec(filename)
+B=((a*t2)/6)+((b*t1)/18)*3;
+A=((b*t1)/18)*6;
+L1=Load/4;
+x=L1/a;
+d=b/3;
+E=1,G=GE;
+mu=(((G*t1)/(d*E))*((A+B)/(A*B)))^0.5;
+D=(-x*A)/((A+B)*mu*cosh(L*0.5*mu));
+function[Sb]=f(z)
+     Sb=((D/B)*sinh(mu*z)) -(x/(A+B))*z;
+endfunction
+funcprot()
+z=linspace(0,L,10*L);
+Sb=feval(z,f);
+plot2d(z,Sb);
+xgrid(3),xtitle('σB versus z','z','σB');
+datatipToggle()
+printf("\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH27/EX27.2/27_2.JPG b/564/CH27/EX27.2/27_2.JPG
new file mode 100755
index 000000000..859a67a89
--- /dev/null
+++ b/564/CH27/EX27.2/27_2.JPG
diff --git a/564/CH27/EX27.2/27_2.sce b/564/CH27/EX27.2/27_2.sce
new file mode 100755
index 000000000..94258dbd2
--- /dev/null
+++ b/564/CH27/EX27.2/27_2.sce
@@ -0,0 +1,29 @@
+pathname=get_absolute_file_path('27_2.sce')
+filename=pathname+filesep()+'27_2data.sci'
+exec(filename)
+clear
+TR=((t*a^3 *b^2)*(2*b +a))/(12*(b+ 2*a));
+J=(1/3)*(2*a +b)*t^3;
+mu=((G*J)/(E*TR))^0.5;
+Ar1=(-3/4)*(a*b/4);
+Ar2=(1/4)*(a*b/4);
+MT=P*2*Ar2;
+X=-MT/(E*TR);
+D=X/(mu*cosh(mu*L));
+F=-D*cosh(0)/mu;
+function[theta]=th(z)//θ
+    theta=((D/mu)*cosh(mu*z))+F;
+endfunction
+z=linspace(0,L,10*L);
+theta=feval(z,th);
+plot2d(z,theta);
+xgrid(3)
+xtitle('angle of twist versus Z','Z','θ')
+datatipToggle();
+printf("\nθ(top): %f rad",th(L));
+MT0=-E*TR*X/3;
+A=(2*a +b)*t;
+Sz1=(P/A)+(MT0/TR)*2*Ar1;//σz1
+Sz2=(P/A)+(MT0/TR)*2*Ar2;//σz2
+printf("\nσz1= σz4= %f N/mm^2",Sz1)
+printf("\nσz3= σz2= %f N/mm^2",Sz2)
+\ No newline at end of file
diff --git a/564/CH3/EX3.1/3_1.sce b/564/CH3/EX3.1/3_1.sce
new file mode 100755
index 000000000..79205f5ee
--- /dev/null
+++ b/564/CH3/EX3.1/3_1.sce
@@ -0,0 +1,35 @@
+pathname=get_absolute_file_path('3_1.sce')
+filename=pathname+filesep()+'3_1data.sci'
+exec(filename)
+clear
+J= (%pi*R^4)/2;
+x=poly(0,"x");
+y=poly(0,"y");
+t= T*R/J;// τ
+twistrate= T/(G*J);//rate of twist
+deff("[Tzy]=f(x)","Tzy=T*x/J");//τzy
+deff("[Tzx]=f1(y)","Tzx=-T*y/J");//τzx
+funcprot(0);
+function sh=shear(x,y),
+    sh=(T/J)*(x^2 +y^2)^0.5;
+endfunction
+t=linspace(-5,5,50);
+clf();
+sh=feval(t,t,shear);
+plot3d1(t,t,sh);
+xtitle('STRESS DISTRIBUTION','X','Y','τ');
+x=[-R:0.01:R];
+y=[-R:0.01:R];
+xset('window',1)
+subplot(2,1,1);
+fplot2d(x,f);
+xgrid(3);
+xtitle( 'shear stress', ' -x- ', ' τzy');
+subplot(2,1,2);
+fplot2d(y,f1,[2]);
+xgrid(3);
+xtitle( 'shear stress', ' -y-', ' τzx');
+datatipToggle();
+printf("\nrate of twist: %f",twistrate);
+printf("\nτ= %f n/mm^2",t);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH3/EX3.1/3_1a.JPG b/564/CH3/EX3.1/3_1a.JPG
new file mode 100755
index 000000000..d5fa46275
--- /dev/null
+++ b/564/CH3/EX3.1/3_1a.JPG
diff --git a/564/CH3/EX3.1/3_1b.JPG b/564/CH3/EX3.1/3_1b.JPG
new file mode 100755
index 000000000..ebb0132ef
--- /dev/null
+++ b/564/CH3/EX3.1/3_1b.JPG
diff --git a/564/CH3/EX3.2/3_2.sce b/564/CH3/EX3.2/3_2.sce
new file mode 100755
index 000000000..957727b2f
--- /dev/null
+++ b/564/CH3/EX3.2/3_2.sce
@@ -0,0 +1,32 @@
+pathname=get_absolute_file_path('3_2.sce')
+filename=pathname+filesep()+'3_2data.sci'
+exec(filename)
+clear
+J= (%pi*a^3 *b^3)/(a^2 +b^2);
+twistrate= T/(G*J);//rate of twist
+x=poly(0,"x");
+y=poly(0,"y");
+deff("[Tzy]=f(x)","Tzy=2*T*x/(%pi*b*a^3)");//τzy
+deff("[Tzx]=f1(y)","Tzx=-2*T*y/(%pi*a*b^3)");//τzx
+function w=wrap(x,y),
+    w=((T*(b^2 -a^2))/(%pi*a^3 *b^3 *G))*x*y*1000;//warping
+endfunction
+t=linspace(-a,a,10*a);
+clf();
+w=feval(t,t,wrap);
+plot3d1(t,t,w);
+xtitle('wraping of cross section','x','y','w * 10^3');
+x=[-a:0.01:a];
+y=[-b:0.01:b];
+xset('window',1)
+subplot(2,1,1);
+fplot2d(x,f);
+xgrid(3);
+xtitle( 'shear stress', ' -x- ', ' τzy');
+subplot(2,1,2);
+fplot2d(y,f1,[2]);
+xgrid(3);
+xtitle( 'shear stress', ' -y-', ' τzx');
+datatipToggle();
+printf("\nrate of twist: %f",twistrate);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH3/EX3.2/3_2a.JPG b/564/CH3/EX3.2/3_2a.JPG
new file mode 100755
index 000000000..2460dc1e9
--- /dev/null
+++ b/564/CH3/EX3.2/3_2a.JPG
diff --git a/564/CH3/EX3.2/3_2b.JPG b/564/CH3/EX3.2/3_2b.JPG
new file mode 100755
index 000000000..7a2845bdd
--- /dev/null
+++ b/564/CH3/EX3.2/3_2b.JPG
diff --git a/564/CH4/EX4.1/4_1.sce b/564/CH4/EX4.1/4_1.sce
new file mode 100755
index 000000000..0b602cf42
--- /dev/null
+++ b/564/CH4/EX4.1/4_1.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('4_1.sce')
+filename=pathname+filesep()+'4_1data.sci'
+exec(filename)
+Ra= W*(L-a)/L;
+Rb= W*a/L;
+printf("\nRA: %f N",Ra);
+printf("\nRB: %f N",Rb);
+\ No newline at end of file
diff --git a/564/CH4/EX4.2/4_2.JPG b/564/CH4/EX4.2/4_2.JPG
new file mode 100755
index 000000000..ae479fb88
--- /dev/null
+++ b/564/CH4/EX4.2/4_2.JPG
diff --git a/564/CH4/EX4.2/4_2.sce b/564/CH4/EX4.2/4_2.sce
new file mode 100755
index 000000000..bf9549940
--- /dev/null
+++ b/564/CH4/EX4.2/4_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('4_2.sce')
+filename=pathname+filesep()+'4_2data.sci'
+exec(filename)
+M=(W*a*(L-a))/(L);
+deff("[y]=f(x)","y=(W*x*(L-x))/(L)")//manding moment
+x=[0:0.05:L];
+fplot2d(x,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'Banding Moment versus a', ' -a- ', '-M-');
+printf("\nMB: %f N.m",M);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH4/EX4.3/4_3.sce b/564/CH4/EX4.3/4_3.sce
new file mode 100755
index 000000000..f95995f20
--- /dev/null
+++ b/564/CH4/EX4.3/4_3.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('4_3.sce')
+filename=pathname+filesep()+'4_3data.sci'
+exec(filename)
+Fba=Loadc*(CD/BD);//FBA
+printf("\nFBA: %f KN",Fba);
+\ No newline at end of file
diff --git a/564/CH4/EX4.4/4_4.JPG b/564/CH4/EX4.4/4_4.JPG
new file mode 100755
index 000000000..88c5468c7
--- /dev/null
+++ b/564/CH4/EX4.4/4_4.JPG
diff --git a/564/CH4/EX4.4/4_4.sce b/564/CH4/EX4.4/4_4.sce
new file mode 100755
index 000000000..5dfbe5043
--- /dev/null
+++ b/564/CH4/EX4.4/4_4.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('4_4.sce')
+filename=pathname+filesep()+'4_4data.sci'
+exec(filename)
+clear
+Vb=(w*L^4)/(8*E*I);//υB
+deff("[y]=f(x)","y= -((w*x^2)/(24*E*I))*(6*L^2 -4*L*x +x^2)");//deflection
+x=[0:0.05:L];
+fplot2d(x,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'vertical deflection versus x', ' -x- ', 'vertical deflection');
+printf("\nυB: %f m//",Vb);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH4/EX4.5/4_5.sce b/564/CH4/EX4.5/4_5.sce
new file mode 100755
index 000000000..50ce46195
--- /dev/null
+++ b/564/CH4/EX4.5/4_5.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('4_5.sce')
+filename=pathname+filesep()+'4_5data.sci'
+exec(filename)
+clear
+theta=(W*L^2)/(16*E*I);//θA
+printf("\nθA: %f",theta);
+\ No newline at end of file
diff --git a/564/CH4/EX4.6/4_6.sce b/564/CH4/EX4.6/4_6.sce
new file mode 100755
index 000000000..a44767450
--- /dev/null
+++ b/564/CH4/EX4.6/4_6.sce
@@ -0,0 +1,37 @@
+pathname=get_absolute_file_path('4_6.sce')
+filename=pathname+filesep()+'4_6data.sci'
+exec(filename)
+L=[(h^2 +AB^2)^0.5;BC;(h^2 +CD^2)^0.5;CD;BC;AB;h;(h^2 +BC^2)^0.5;h];//column 2
+A=atan(h/AB),B=atan(h/BC),D=atan(h/CD);
+lFC=Load2;
+P=[
+   0 0 sin(D) 0 0 0 0 sin(B);
+   cos(A) 0 0 0 0 1 0 0;
+   cos(A) -1 0 0 0 0 0 0;
+   sin(A) 0 0 0 0 0 1 0;
+   0 0 0 -1 1 0 0 0;
+   0 0 0 0 1 -1 0 cos(B);
+   0 0 0 0 0 0 1 -sin(B);
+   0 1 -cos(D) 0 0 0 0 cos(B);];
+X=[-Load2;0;0;-Load1;0;0;Load1;0];
+D1=[0;0;0;1;1;1;0;0;0];//F1,D
+FC=0;
+P1=[cos(A) 0 0 0 0 1 0 0;
+    0 0 cos(D) 1 0 0 0 0;
+    cos(A) -1 0 0 0 0 0 0;
+    sin(A) 0 0 0 0 0 1 0;
+    0 0 0 -1 1 0 0 0;
+    0 0 0 0 1 -1 0 cos(B);
+    0 0 0 0 0 0 1 sin(B);
+    0 0 1 0 0 0 0 1]
+X1=[0;0;0;0;0;0;1;0]
+FA=[inv(P)*X;lFC];//FA
+B1=[inv(P1)*X1;FC];//F1,B
+for i=1:9
+    K1(i)=FA(i)*B1(i)*L(i);//FA.F1,B.L
+    K2(i)=FA(i)*D1(i)*L(i);//FA.F1,D.L
+end
+delB=(sum(K1)*10^3)/(E*CS);
+delD=(sum(K2)*10^3)/(E*CS);
+printf("\nδB,v: %f mm",delB);
+printf("\nδD,h: %f mm",delD);
diff --git a/564/CH5/EX5.1/5_1.sce b/564/CH5/EX5.1/5_1.sce
new file mode 100755
index 000000000..0dcb9281c
--- /dev/null
+++ b/564/CH5/EX5.1/5_1.sce
@@ -0,0 +1,37 @@
+pathname=get_absolute_file_path('5_1.sce')
+filename=pathname+filesep()+'5_1data.sci'
+exec(filename)
+L=[(h^2 +AB^2)^0.5;BC;(h^2 +CD^2)^0.5;CD;BC;AB;h;(h^2 +BC^2)^0.5;h];
+A=atan(h/AB),B=atan(h/BC),D=atan(h/CD);
+lFC=Load2;
+P=[
+   0 0 sin(D) 0 0 0 0 sin(B);
+   cos(A) 0 0 0 0 1 0 0;
+   cos(A) -1 0 0 0 0 0 0;
+   sin(A) 0 0 0 0 0 1 0;
+   0 0 0 -1 1 0 0 0;
+   0 0 0 0 1 -1 0 cos(B);
+   0 0 0 0 0 0 1 -sin(B);
+   0 1 -cos(D) 0 0 0 0 cos(B);];
+X=[-Load2;0;0;-Load1;0;0;Load1;0];
+D1=[0;0;0;1;1;1;0;0;0];
+FC=0;
+P1=[cos(A) 0 0 0 0 1 0 0;
+    0 0 cos(D) 1 0 0 0 0;
+    cos(A) -1 0 0 0 0 0 0;
+    sin(A) 0 0 0 0 0 1 0;
+    0 0 0 -1 1 0 0 0;
+    0 0 0 0 1 -1 0 cos(B);
+    0 0 0 0 0 0 1 sin(B);
+    0 0 1 0 0 0 0 1]
+X1=[0;0;0;0;0;0;1;0]
+FA=[inv(P)*X;lFC];
+B1=[inv(P1)*X1;FC];
+for i=1:9
+    K1(i)=FA(i)*B1(i)*L(i);
+    K2(i)=FA(i)*D1(i)*L(i);
+end
+delB=(sum(K1)*10^3)/(E*CS);
+delD=(sum(K2)*10^3)/(E*CS);
+printf("\nδB,v: %f mm",delB);
+printf("\nδD,h: %f mm",delD);
diff --git a/564/CH5/EX5.10/5_10.sce b/564/CH5/EX5.10/5_10.sce
new file mode 100755
index 000000000..7e5beb393
--- /dev/null
+++ b/564/CH5/EX5.10/5_10.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('5_10.sce')
+filename=pathname+filesep()+'5_10data.sci'
+exec(filename)
+c=((a*a +b*b)^0.5);
+L=[a;b;a;b;(a*a +b*b)^0.5;(a*a +b*b)^0.5];
+F1=[a/b;1;a/b;1;-c/b;-c/b]
+for i=1:6
+    X(i)=F1(i)*F1(i)*L(i);
+end
+a11=sum(X)/(A*E);
+X1=-alpha*b*T/a11;
+Fa=[X1*a/b;X1;X1*a/b;X1;-X1*c/b;-X1*c/b]
+printf("\nX1: %f N",X1);
+disp("Forces in the Mambrane are (AB,BC,CD,DA,AC,DB):");
+printf("\n %f KN",Fa);
+\ No newline at end of file
diff --git a/564/CH5/EX5.11/5_11.sce b/564/CH5/EX5.11/5_11.sce
new file mode 100755
index 000000000..741274fa1
--- /dev/null
+++ b/564/CH5/EX5.11/5_11.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('5_11.sce')
+filename=pathname+filesep()+'5_11data.sci'
+exec(filename)
+VB=(W*L^3)/(48*EI);
+printf("\nVB: %f m",VB);
+\ No newline at end of file
diff --git a/564/CH5/EX5.12/5_12.sce b/564/CH5/EX5.12/5_12.sce
new file mode 100755
index 000000000..daad7d644
--- /dev/null
+++ b/564/CH5/EX5.12/5_12.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('5_12.sce')
+filename=pathname+filesep()+'5_12data.sci'
+exec(filename)
+v1= (L1/L0)*V(3);
+v2= (L2/L0)*V(4);
+v3= v1+v2;
+theta=atan(v3/(D(9)-D(6)));
+printf("\nθB: %f rad",theta);
+\ No newline at end of file
diff --git a/564/CH5/EX5.2/5_2.sce b/564/CH5/EX5.2/5_2.sce
new file mode 100755
index 000000000..7304a2e5e
--- /dev/null
+++ b/564/CH5/EX5.2/5_2.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('5_2.sce')
+filename=pathname+filesep()+'5_2data.sci'
+exec(filename)
+delB=(119*w*L^4)/(24576*EI);
+delC=(5*w*L^4)/(384*EI);
+printf("\ndelB: %f m",delB);
+printf("\ndelC: %f m",delC);
+\ No newline at end of file
diff --git a/564/CH5/EX5.3/5_3.sce b/564/CH5/EX5.3/5_3.sce
new file mode 100755
index 000000000..2efbdd8cf
--- /dev/null
+++ b/564/CH5/EX5.3/5_3.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('5_3.sce')
+filename=pathname+filesep()+'5_3data.sci'
+exec(filename)
+clear
+L1=[L*tan(theta);L;L*tan(theta);L;L];
+A=[AC;A;AC;A;A];
+F1=[Load*cos(theta);Load*sin(theta);0;0;0];
+F2=[-sin(theta);cos(theta);-sin(theta);1;cos(theta)];
+for i=1:5
+    X(i)=L1(i)/A(i);
+    X1(i)=F1(i)*F2(i)*X(i);
+    X2(i)=F2(i)*F2(i)*X(i);
+end
+R=-sum(X1)/sum(X2);
+F=F1+ R*F2;
+disp("Force (AC,CB,BD,CD,AD):");
+printf("\n %f N",F);
+\ No newline at end of file
diff --git a/564/CH5/EX5.4/5_4.sce b/564/CH5/EX5.4/5_4.sce
new file mode 100755
index 000000000..6121bb9c3
--- /dev/null
+++ b/564/CH5/EX5.4/5_4.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('5_4.sce')
+filename=pathname+filesep()+'5_4data.sci'
+exec(filename)
+clear
+AC=(AB^2 +BC^2)^0.5;
+delBC=BC*delT*alpha;
+F1=[AB/BC;1;AB/BC;1;-AC/BC;-AC/BC];
+L=[AB;BC;AB;BC;AC;AC];
+for i=1:6
+    X(i)=L(i)*F1(i)*F1(i);
+end
+R=(-delBC*A*E)/sum(X);
+F=R*F1;
+disp("Force (AB,BC,CD,DA,AC,DB):");
+printf("\n %f N",F);
+\ No newline at end of file
diff --git a/564/CH5/EX5.5/5_5.sce b/564/CH5/EX5.5/5_5.sce
new file mode 100755
index 000000000..74e4e4efa
--- /dev/null
+++ b/564/CH5/EX5.5/5_5.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('5_5.sce')
+filename=pathname+filesep()+'5_5data.sci'
+exec(filename)
+clear
+L1=0.5*L*ones(7,1);
+A1=[Ab;Ab;A;A;A;A;A];
+R=(11*3^0.5 *P*A*Ab*L^2)/(48*(L*L*A*Ab + 4*I*(A+ 10*Ab)));
+delF=[-cos(theta);-cos(theta);1;1;-1;-1;1];
+F=R*delF;
+disp(R,"R:");
+disp("Force (AB,BC,CD,DE,BD,EB,AE):");
+printf("\n %f N",F)
diff --git a/564/CH5/EX5.6/5_6.sce b/564/CH5/EX5.6/5_6.sce
new file mode 100755
index 000000000..a8a220fcf
--- /dev/null
+++ b/564/CH5/EX5.6/5_6.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('5_6.sce')
+filename=pathname+filesep()+'5_6data.sci'
+exec(filename)
+SA=0.187*M0/r;
+SB=0.44*M0/r;
+SC=0.373*M0/r;
+printf("\nSA: %f N ",SA)
+printf("\nSB: %f N",SB)
+printf("\nSC:%f  N",SC)
+printf("\nM1: %f N.mm",SA*r)
+printf("\nM2: %f N.mm",SB*r)
+printf("\nM3: %f N.mm",SC*r)
+\ No newline at end of file
diff --git a/564/CH5/EX5.7/5_7.sce b/564/CH5/EX5.7/5_7.sce
new file mode 100755
index 000000000..f00563639
--- /dev/null
+++ b/564/CH5/EX5.7/5_7.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('5_7.sce')
+filename=pathname+filesep()+'5_7data.sci'
+exec(filename)
+clear
+delY=(W*L^4)*((11/(24*EI))+ 1/(2*GJ));
+delZ=(W*L^4)*((1/(6*EI))+ 1/(2*GJ));
+printf("\ndel Y %f mm",delY);
+printf("\ndel Z %f mm",delZ);
+\ No newline at end of file
diff --git a/564/CH5/EX5.8/5_8.sce b/564/CH5/EX5.8/5_8.sce
new file mode 100755
index 000000000..ff68c2399
--- /dev/null
+++ b/564/CH5/EX5.8/5_8.sce
@@ -0,0 +1,14 @@
+pathname=get_absolute_file_path('5_8.sce')
+filename=pathname+filesep()+'5_8data.sci'
+exec(filename)
+L1=[L;L;L;1.41*L;1.41*L;L];
+F0=[0 0 -P 0 1.41*P 0];
+F1=[-0.71;-0.71;-0.71;1;1;-0.71];
+for i=1:6
+    X(i)=F0(i)*F1(i)*L1(i);
+    X1(i)=F1(i)*F1(i)*L1(i);
+end
+XBD= -sum(X)/sum(X1);
+Fa=[-XBD/1.41;-XBD/1.41;((XBD/1.41)+1);XBD;-(((XBD/1.41)+1))*1.41;-XBD/1.41];
+disp("Forces in the Mambrane are:");
+printf("\n %f N",Fa)
+\ No newline at end of file
diff --git a/564/CH5/EX5.9/5_9.sce b/564/CH5/EX5.9/5_9.sce
new file mode 100755
index 000000000..3cb00572e
--- /dev/null
+++ b/564/CH5/EX5.9/5_9.sce
@@ -0,0 +1,24 @@
+pathname=get_absolute_file_path('5_9.sce')
+filename=pathname+filesep()+'5_9data.sci'
+exec(filename)
+M=2^0.5;
+L1=[L;M*L;L;L;M*L;M*L;L];
+F0=[P;0;0;0;0;-M*P;0];
+F1=[-0.71;0;0;-0.71;1;1;-0.71];
+F1R2=[-2;-M;1;1;0;M;0];
+for i=1:7
+    X(i)=F0(i)*F1(i)*L1(i);
+    X1(i)=F0(i)*F1R2(i)*L1(i);
+    X2(i)=F1(i)*F1(i)*L1(i);
+    X3(i)=F1R2(i)*F1R2(i)*L1(i);
+    X4(i)=F1(i)*F1R2(i)*L1(i);
+end
+X5=[sum(X2) sum(X4);
+    sum(X4) sum(X3)];
+X6=-[sum(X);sum(X1)]
+X7=inv(X5)*X6;
+printf("\nX1: %f KN",X7(1,:));
+printf("\nR2: %f KN",X7(2,:))
+Fa=[X7(1,:)*M-X7(2,:)*(1+(1/M));-X7(2,:)*M;X7(2,:);X7(2,:)-X7(1,:)/M;X7(1,:);-X7(2,:)*(1+(1/M));-X7(1,:)/M]
+disp("Forces in the Mambrane are:");
+printf("\n %f KN",Fa);
+\ No newline at end of file
diff --git a/564/CH6/EX6.1/6_1.JPG b/564/CH6/EX6.1/6_1.JPG
new file mode 100755
index 000000000..3d3c44983
--- /dev/null
+++ b/564/CH6/EX6.1/6_1.JPG
diff --git a/564/CH6/EX6.1/6_1.sce b/564/CH6/EX6.1/6_1.sce
new file mode 100755
index 000000000..29664b24c
--- /dev/null
+++ b/564/CH6/EX6.1/6_1.sce
@@ -0,0 +1,41 @@
+pathname=get_absolute_file_path('6_1.sce')
+filename=pathname+filesep()+'6_1data.sci'
+exec(filename)
+//theta=θ,lambda=λ,mu=μ
+theta12=0,lambda12=cos(theta12),mu12=sin(theta12);
+theta13=%pi/2,lambda13=cos(theta13),mu13=sin(theta13);
+theta23=%pi/2 + atan(L13/L12);lambda23=cos(theta23),mu23=sin(theta23);
+L23=(L12^2 +L13^2)^0.5;
+LL12=lambda12^2,LM12=lambda12*mu12,MM12=mu12^2;
+LL13=lambda13^2,LM13=lambda13*mu13,MM13=mu13^2;
+LL23=lambda23^2,LM23=lambda23*mu23,MM23=mu23^2;
+K220=[LL12 LM12;LM12 MM12];
+K222=(L12/L23)*[LL23 LM23;LM23 MM23];
+K210=[-LL12 -LM12;-LM12 -MM12];
+K232=(L12/L23)*[-LL23 -LM23;-LM23 -MM23];
+K=K220+K222;
+F=[FX2;FY2];
+V=(L12/AE)*inv(K)*F;
+u2=V(1,:),v2=V(2,:);
+u1=0;v1=0;u3=0;v3=0;
+F1=[K210';K232']*inv(K)*F;
+S12=(AE/L12)*[lambda12 mu12]*[u2-u1;v2-v1];
+S13=(AE/L13)*[lambda13 mu13]*[u3-u1;v3-v1];
+S23=(AE/L23)*[lambda23 mu23]*[u3-u2;v3-v2];
+x1=[-L12,-L12],y1=[0,L13];
+x2=[-L12,u2],y2=[0,v2];
+x3=[-L12,u2],y3=[L13,v2];
+plot(x1,y1,x2,y2,x3,y3);
+xgrid(3);
+datatipToggle();
+xtitle( 'framework after deformation', ' -x- ', ' -y- ',boxed = 1);
+printf("\n\nu2: %f m",u2);
+printf("\nv2: %f m",v2);
+printf("\nFx,1: %f N",F1(1,:));
+printf("\nFy,1: %f N",F1(2,:));
+printf("\nFx,3: %f N",F1(3,:));
+printf("\nFy,3: %f N",F1(4,:));
+printf("\nS12: %f N",S12);
+printf("\nS13: %f N",S13);
+printf("\nS23: %f N",S23);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH6/EX6.2/6_2.JPG b/564/CH6/EX6.2/6_2.JPG
new file mode 100755
index 000000000..8ef98acb3
--- /dev/null
+++ b/564/CH6/EX6.2/6_2.JPG
diff --git a/564/CH6/EX6.2/6_2.sce b/564/CH6/EX6.2/6_2.sce
new file mode 100755
index 000000000..aa9952125
--- /dev/null
+++ b/564/CH6/EX6.2/6_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('6_2.sce')
+filename=pathname+filesep()+'6_2data.sci'
+exec(filename)
+F=[-W;M;0;0];
+P=EI*[(27/(2*L^3)) (9/(2*L^2)) (6/L^2) -(3/(2*L^2));
+   (9/(2*L^2)) (6/L) (2/L) (1/L);
+   (6/L^2) (2/L) (4/L) 0;
+   -(3/(2*L^2)) (1/L) 0 (2/L)];
+V=inv(P)*F;
+Sy12= EI*(-(6*V(3)/(L^2))-(12*V(1)/(L^3)) -(6*V(2)/(L*L)));//Sy,12
+deff("[M12]=f(x)","M12=EI*(((-6*x/L*L)+(4/L))*V(3) + ((-12*x/L^3)+(6/L^2))*V(1) +((-6*x/L*L)+(2/L))*V(2))");
+x=[0:0.05:L];
+fplot2d(x,f);
+xgrid(3);
+xtitle( 'M12 versus x', ' -x- ', 'M12');
+datatipToggle();
+printf("\nυ2: %f m",V(1));//v2
+printf("\nθ2: %f ",V(2));//θ2
+printf("\nθ2: %f ",V(3));//θ1
+printf("\nθ2: %f ",V(4));//θ3
+printf("\nSy,12: %f",Sy12);//Sy,12
+\ No newline at end of file
diff --git a/564/CH6/EX6.4/6_4.sce b/564/CH6/EX6.4/6_4.sce
new file mode 100755
index 000000000..b40aa1f33
--- /dev/null
+++ b/564/CH6/EX6.4/6_4.sce
@@ -0,0 +1,25 @@
+pathname=get_absolute_file_path('6_4.sce')
+filename=pathname+filesep()+'6_4data.sci'
+exec(filename)
+P=[1 P1(1) P1(2) P1(1)*P1(2);
+   1 P2(1) P2(2) P2(1)*P2(2);
+   1 P3(1) P3(2) P3(1)*P3(2);
+   1 P4(1) P4(2) P4(1)*P4(2)];
+alpha1=inv(P)*u;
+alpha2=inv(P)*v;
+alpha=[alpha1;alpha2];
+deff("[Ex]=f(y)","Ex=alpha(2)+ y*alpha(4)");//εx
+deff("[Ey]=f1(x)","Ey=alpha(7)+ x*alpha(8)");//εy
+function[G]=Gxy(x,y)//γxy
+    G=x*alpha(4) +y*alpha(8) +alpha(3)+alpha(6);
+endfunction
+//at the centre
+Pc(1)=(P1(1)+P3(1))/2;
+Pc(2)=(P1(2)+P3(2))/2;
+Sy=(E/(1-V^2))*(f1(Pc(2)) +V*f(Pc(1)));
+Sx=(E/(1-V^2))*(f(Pc(1)) +V*f1(Pc(2)));
+gxy=Gxy(0,0);
+txy=(E/(1-V^2))*0.5*(1-V)*gxy;//τxy
+printf("\nσx: %f N/mm^2",Sx);
+printf("\nσy: %f N/mm^2",Sy);
+printf("\nτxy: %f N/mm^2",txy);
+\ No newline at end of file
diff --git a/564/CH7/EX7.1/7_1.JPG b/564/CH7/EX7.1/7_1.JPG
new file mode 100755
index 000000000..74a74c7bc
--- /dev/null
+++ b/564/CH7/EX7.1/7_1.JPG
diff --git a/564/CH7/EX7.1/7_1.sce b/564/CH7/EX7.1/7_1.sce
new file mode 100755
index 000000000..2e018a1e6
--- /dev/null
+++ b/564/CH7/EX7.1/7_1.sce
@@ -0,0 +1,28 @@
+pathname=get_absolute_file_path('7_1.sce')
+filename=pathname+filesep()+'7_1data.sci'
+exec(filename)
+D=(E*t^3)/(12*(1-v^2));
+X=(16*q0/(D*%pi^6))/(((1/a^2)+(1/b^2))^2);
+X1=(16*q0/%pi^4)*((1/a^2)+(v/b^2))/(((1/a^2)+(1/b^2))^2);
+X2=(16*q0/%pi^4)*((v/a^2)+(1/b^2))/(((1/a^2)+(1/b^2))^2);
+function[w]=f(x,y),//taking first term only
+    w=X*(sin(%pi*x/a))*sin(%pi*y/b);
+endfunction
+x=linspace(0,a,10*a);
+y=linspace(0,b,10*b);
+w=feval(x,y,f);
+plot3d1(x,y,w);
+Wmax=(0.0443*q0*a^4)/(E*t^3);//from summation of first 4 terms
+printf("\nWmax: %f mm",Wmax);
+function[Mx,My]=f1(x,y)//taking first term only
+    Mx=X1*(sin(%pi*x/a))*sin(%pi*y/b);
+    My=X2*(sin(%pi*x/a))*sin(%pi*y/b);
+endfunction
+function[Sx,Sy]=f2(x,y,z)//taking first term only
+    Sx=12*X1*(sin(%pi*x/a))*sin(%pi*y/b)*z/t^3;
+    Sy=12*X2*(sin(%pi*x/a))*sin(%pi*y/b)*z/t^3;
+endfunction
+Mxmax=0.0479*q0*a^2;//from summation of first 5 terms
+printf("\nMx,max= My,max= %f N.mm",Mxmzx);
+Sxmax=(0.287*q0*a^2)/t^2;
+printf("\nσx,max= σy,max= %f N.mm",Sxmax);
diff --git a/564/CH7/EX7.3/7_3.JPG b/564/CH7/EX7.3/7_3.JPG
new file mode 100755
index 000000000..2f45a12a3
--- /dev/null
+++ b/564/CH7/EX7.3/7_3.JPG
diff --git a/564/CH7/EX7.3/7_3.sce b/564/CH7/EX7.3/7_3.sce
new file mode 100755
index 000000000..c2646e723
--- /dev/null
+++ b/564/CH7/EX7.3/7_3.sce
@@ -0,0 +1,14 @@
+pathname=get_absolute_file_path('7_3.sce')
+filename=pathname+filesep()+'7_3data.sci'
+exec(filename)
+D=(E*t^3)/(12*(1-v^2));
+X=(16*q0/(D*%pi^6))/((((1/a^2)+(1/b^2))^2)+(Nx/(D*a*a*%pi^2)));
+function[w]=f(x,y),//taking first term only
+    w=X*(sin(%pi*x/a))*sin(%pi*y/b);
+endfunction
+x=linspace(0,a,10*a);
+y=linspace(0,b,10*b);
+w=feval(x,y,f);
+plot3d1(x,y,w);
+datatipToggle();
+printf("\nclick on the point to view its coordinate on the plot")
+\ No newline at end of file
diff --git a/564/CH7/EX7.4/7_4.JPG b/564/CH7/EX7.4/7_4.JPG
new file mode 100755
index 000000000..c966c9020
--- /dev/null
+++ b/564/CH7/EX7.4/7_4.JPG
diff --git a/564/CH7/EX7.4/7_4.sce b/564/CH7/EX7.4/7_4.sce
new file mode 100755
index 000000000..4bcec2022
--- /dev/null
+++ b/564/CH7/EX7.4/7_4.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('7_4.sce')
+filename=pathname+filesep()+'7_4data.sci'
+exec(filename)
+clear
+D=(E*t^3)/(12*(1-v^3));
+A11= (16*q0*(a^4)*(b^4))/((%pi^6)*D*((a*a +b*b)^2));
+function[w]=f(x,y)
+    w=A11*(sin(%pi*x/a))*(sin(%pi*y/b));
+endfunction
+funcprot()
+x=linspace(0,a,10*a);
+y=linspace(0,b,10*b);
+clf();
+w=feval(x,y,f);
+plot3d1(x,y,w);
+datatipToggle()
+printf("\nA11= Wmax: %f mm",f(a/2,b/2));
+\ No newline at end of file
diff --git a/564/CH8/EX8.2/8_2.JPG b/564/CH8/EX8.2/8_2.JPG
new file mode 100755
index 000000000..58bb74203
--- /dev/null
+++ b/564/CH8/EX8.2/8_2.JPG
diff --git a/564/CH8/EX8.2/8_2.sce b/564/CH8/EX8.2/8_2.sce
new file mode 100755
index 000000000..27034e601
--- /dev/null
+++ b/564/CH8/EX8.2/8_2.sce
@@ -0,0 +1,19 @@
+pathname=get_absolute_file_path('8_2.sce')
+filename=pathname+filesep()+'8_2data.sci'
+exec(filename)
+mu=(P/(E*I))^0.5;//μ
+deff("[v]=f(z)","v=e*((cos(mu*(z-L/2))/cos(mu*L/2))-1)");//v
+deff("[m]=f1(z)","m=P*(e+e*((cos(mu*(z-L/2))/cos(mu*L/2))-1))");//M
+funcprot(0);
+z=[0:0.05:L];
+subplot(2,1,1);
+fplot2d(z,f);
+xgrid(3);
+xtitle( 'deflection', ' -z- ', '-v-');
+subplot(2,1,2);
+fplot2d(z,f1);
+xgrid(3);
+xtitle( 'bending moment', ' -z- ', ' -M-');
+datatipToggle();
+printf("\nmaximum deflection: %f",f(L/2));
+printf("\nM(max)= %f N.m",f1(L/2));
+\ No newline at end of file
diff --git a/564/CH8/EX8.3/8_3.sce b/564/CH8/EX8.3/8_3.sce
new file mode 100755
index 000000000..a8e224b3f
--- /dev/null
+++ b/564/CH8/EX8.3/8_3.sce
@@ -0,0 +1,20 @@
+pathname=get_absolute_file_path('8_3.sce')
+filename=pathname+filesep()+'8_3data.sci'
+exec(filename)
+A=t*(2*a +b);
+Ixx=2*a*t*(b/2)^2 +(t*b^3)/12;
+Iyy=(2*t*a^3)/12;
+I0=Ixx+Iyy;
+J=(2*a*t^3)/3 +(b*t^3)/3;
+tau=(t*a^3 *b^2)/24;
+Px=(%pi^2 *E*Ixx/L^2);//PCR(xx)
+Py=(%pi^2 *E*Iyy/L^2);//PCR(yy)
+Pth=(A/I0)*(G*J +(%pi^2 *E*tau)/L^2);//PCR(θ)
+P=poly(0,'P');
+y=(I0/A)*(P^2 -(Px +Pth)*P +Px*Pth);
+m=roots(y);
+Ptb1=m(1,:)
+Ptb2=m(2,:)
+Load=[Px;Py;Pth];
+printf("\nbuckling Load: %f N",min(Load));
+printf("\nflexural-torsional buckling Load: %f N",roots(y));
+\ No newline at end of file
diff --git a/564/CH8/EX8.4/8_4.sce b/564/CH8/EX8.4/8_4.sce
new file mode 100755
index 000000000..b1d2c9477
--- /dev/null
+++ b/564/CH8/EX8.4/8_4.sce
@@ -0,0 +1,31 @@
+pathname=get_absolute_file_path('8_4.sce')
+filename=pathname+filesep()+'8_4data.sci'
+exec(filename)
+clear
+Xc=(a2^2)/(b+ 2*a2);//x bar
+Xs=-(Xc+((3*a2^2)/(b*(1+ 6*a2/b))));
+A= 2*a2*t + b*t;
+Ixx=2*a2*t*(b/2)^2 +(t*b^3)/12;
+Iyy=b*t*Xc^2 +(2*t/3)*((a2-Xc)^3 +Xc^3);
+I0=Ixx+Iyy+ A*Xs^2;
+J=(2*a2*t^3)/3 +(b*t^3)/3;
+tau= 0.1244*t*a2^3 *b^2;
+Px=((%pi^2)*E*Ixx/(L^2));//PCR(xx)
+Py=((%pi^2)*E*Iyy/(L^2));//PCR(yy)
+Pth=(A/I0)*(G*J +((%pi^2)*E*tau)/(L^2));//PCR(θ)
+a=(1-(A*Xs^2)/I0),b=-(Px +Pth),c=Px*Pth;
+P1=(-b +(b*b -4*a*c)^0.5)/(2*a);
+P2=(-b -(b*b -4*a*c)^0.5)/(2*a);
+Load=[Px;Py;Pth;P1;P2];
+minimum=Load(1);
+for i=2:5
+    if(Load(i)<Load(i-1)) then
+        minimum=Load(i);
+        a1=i;
+    end
+end
+if(a1>3) then
+    printf("\nflexural-torsional buckling will happen.\nand buckling load is: %f N",minimum);
+else
+    printf("\nuncoupled buckling will happen\nand buckling Load is %f N",minimum);
+end
diff --git a/564/CH9/EX9.1/9_1.sce b/564/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..1ddde273f
--- /dev/null
+++ b/564/CH9/EX9.1/9_1.sce
@@ -0,0 +1,25 @@
+pathname=get_absolute_file_path('9_1.sce')
+filename=pathname+filesep()+'9_1data.sci'
+exec(filename)
+alpha= atan(((1+ t*d/(2*Af))/(1+ t*b/As))^0.25);//α
+Ft=(W*z/d)+(W/(2*tan(alpha)));//FT
+printf("\nFT: %f N",Ft);
+Mmax= (W*tan(alpha)*b^2)/(12*d); 
+Smax= (Mmax/ESM)+(Ft/Af);//σ max
+printf("\nMaximum Stress in Flange: %f N/mm^2",Smax);
+P=(W*b*tan(alpha))/d;
+if(b<1.5*d) then//le
+    Le=d/((4-(2*b/d))^0.5);
+else
+    Le=d;
+end
+Pcr= ((%pi^2)*E*I)/(Le^2);
+printf("\nPcr: %f KN",Pcr);
+printf("\nP: %f KN",P);
+if(P<Pcr) then
+    printf("\nstiffener will not buckle")
+else
+    printf("\nstiffener will buckle")
+end
+
+
diff --git a/564/DEPENDENCIES/10_1data.sci b/564/DEPENDENCIES/10_1data.sci
new file mode 100755
index 000000000..2d7364887
--- /dev/null
+++ b/564/DEPENDENCIES/10_1data.sci
@@ -0,0 +1,3 @@
+m=1;//in Kg
+EI=70000*2000;//given in N.mm^2
+l=2000;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/10_2data.sci b/564/DEPENDENCIES/10_2data.sci
new file mode 100755
index 000000000..2d7364887
--- /dev/null
+++ b/564/DEPENDENCIES/10_2data.sci
@@ -0,0 +1,3 @@
+m=1;//in Kg
+EI=70000*2000;//given in N.mm^2
+l=2000;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/10_4data.sci b/564/DEPENDENCIES/10_4data.sci
new file mode 100755
index 000000000..46cd644cf
--- /dev/null
+++ b/564/DEPENDENCIES/10_4data.sci
@@ -0,0 +1,5 @@
+m=1;//in Kg
+EI=70000*2000;//given in N.mm^2
+L=2000;//in mm
+rho=1;//in Kg/mm^3
+A=1;//in mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/10_5data.sci b/564/DEPENDENCIES/10_5data.sci
new file mode 100755
index 000000000..e04893eb7
--- /dev/null
+++ b/564/DEPENDENCIES/10_5data.sci
@@ -0,0 +1,4 @@
+EI=70000*2000;//given in N.mm^2
+L=2000;//in mm
+rho=1;//in Kg/mm^3
+A=1;//in mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/10_6data.sci b/564/DEPENDENCIES/10_6data.sci
new file mode 100755
index 000000000..2d7364887
--- /dev/null
+++ b/564/DEPENDENCIES/10_6data.sci
@@ -0,0 +1,3 @@
+m=1;//in Kg
+EI=70000*2000;//given in N.mm^2
+l=2000;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/12_1data.sci b/564/DEPENDENCIES/12_1data.sci
new file mode 100755
index 000000000..28faf0910
--- /dev/null
+++ b/564/DEPENDENCIES/12_1data.sci
@@ -0,0 +1,5 @@
+t=2.5;//skin thickness,given in mm
+St=1.2;//straps thickness,given in mm
+d=4;//rivets diameter,given in mm
+Lt=125;//limit tensile stress
+Ls=120;//limit shear stress
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/12_2data.sci b/564/DEPENDENCIES/12_2data.sci
new file mode 100755
index 000000000..b735096a8
--- /dev/null
+++ b/564/DEPENDENCIES/12_2data.sci
@@ -0,0 +1,6 @@
+AD=25;//distence between point A and D,given in mm
+DC=20;//distence between point D and centroid,given in mm
+DG=25;//distence between point G and D,given in mm
+CF=25;//distence between point F and centroid,given in mm
+Load=5000;//load,given(5kN in N)
+CL=75;//distence between the centroid and load
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/14_1data.sci b/564/DEPENDENCIES/14_1data.sci
new file mode 100755
index 000000000..e5f225428
--- /dev/null
+++ b/564/DEPENDENCIES/14_1data.sci
@@ -0,0 +1,5 @@
+W=45;//given in kN
+Wa=4.5;//wight of aircraft aft,given in kN
+v0=25;//given in m/s
+g=9.81;//in m/s^2;
+a=3*g;//given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/14_2data.sci b/564/DEPENDENCIES/14_2data.sci
new file mode 100755
index 000000000..5b89652ee
--- /dev/null
+++ b/564/DEPENDENCIES/14_2data.sci
@@ -0,0 +1,8 @@
+W=250;//given in kN
+Rh=400;//horizontal reaction,given in kN
+Rv=1200;//vertical reaction,given in kN
+Icg=5.65*10^8;//given in N.s^2
+Sh=1;//horizontal distence of CG from main wheels,given in m
+Sv=2.5;//vertical distence of CG from main wheels,given in m
+v0=3.7;//initial vertical velocity,given in m/s
+g=9.81;//in m/s^2;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/14_3data.sci b/564/DEPENDENCIES/14_3data.sci
new file mode 100755
index 000000000..13b0f4a67
--- /dev/null
+++ b/564/DEPENDENCIES/14_3data.sci
@@ -0,0 +1,10 @@
+W=8000;//given in N
+S=14.5;//wing area,given in m^2
+c=1.35;//mean chord,given in m
+n=4.5;
+v=60;//speed,given in m/s
+rho=1.223;//air density,given in kg/m^3
+alpha1=13.75;//from fig. 14.8(a),in degree
+Cmcg1=0.075;//Cm,cg  from fig. 14.8(a)
+alpha2=13.3;//given in degree
+Cmcg2=0.073;//given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/15_1data.sci b/564/DEPENDENCIES/15_1data.sci
new file mode 100755
index 000000000..b4fe95124
--- /dev/null
+++ b/564/DEPENDENCIES/15_1data.sci
@@ -0,0 +1,6 @@
+K=1708;//fracture toughness, given in N/mm^2
+S=175;//given in N/mm^2
+alpha=1;///given
+C=40/10^15;//rete of crack growth, given in mm/cycle
+n=4;//given
+ai=0.2;//initial leangth of crack,given in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_10data.sci b/564/DEPENDENCIES/16_10data.sci
new file mode 100755
index 000000000..df8cf2454
--- /dev/null
+++ b/564/DEPENDENCIES/16_10data.sci
@@ -0,0 +1,3 @@
+W=10;//load in N/m
+L=20;//L in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_11data.sci b/564/DEPENDENCIES/16_11data.sci
new file mode 100755
index 000000000..d8d2d92eb
--- /dev/null
+++ b/564/DEPENDENCIES/16_11data.sci
@@ -0,0 +1,4 @@
+M0=10;//load in N/m
+L=20;//L in m
+b=15;//b in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_12data.sci b/564/DEPENDENCIES/16_12data.sci
new file mode 100755
index 000000000..7c031249e
--- /dev/null
+++ b/564/DEPENDENCIES/16_12data.sci
@@ -0,0 +1,6 @@
+Ixx=1.937D+08;
+Iyy=27005208;
+Ixy=562584521;
+E=20000;
+L=20;
+W=100000
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_13data.sci b/564/DEPENDENCIES/16_13data.sci
new file mode 100755
index 000000000..cf98a373a
--- /dev/null
+++ b/564/DEPENDENCIES/16_13data.sci
@@ -0,0 +1,4 @@
+b=100;//in mm
+tf=2;//in mm
+tw=2.5;//in mm
+d=150;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_14data.sci b/564/DEPENDENCIES/16_14data.sci
new file mode 100755
index 000000000..97cd862f4
--- /dev/null
+++ b/564/DEPENDENCIES/16_14data.sci
@@ -0,0 +1,3 @@
+h=100;//in mm
+t=2;//in mm
+Mx=1000;//N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_15data.sci b/564/DEPENDENCIES/16_15data.sci
new file mode 100755
index 000000000..df2d8eee4
--- /dev/null
+++ b/564/DEPENDENCIES/16_15data.sci
@@ -0,0 +1,5 @@
+a=100;//in mm
+t=2;//in mm
+E=70000;//N/mm^2
+T0=10;//in K
+alpha=0.01;//mm/k
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_16data.sci b/564/DEPENDENCIES/16_16data.sci
new file mode 100755
index 000000000..df2d8eee4
--- /dev/null
+++ b/564/DEPENDENCIES/16_16data.sci
@@ -0,0 +1,5 @@
+a=100;//in mm
+t=2;//in mm
+E=70000;//N/mm^2
+T0=10;//in K
+alpha=0.01;//mm/k
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_1data.sci b/564/DEPENDENCIES/16_1data.sci
new file mode 100755
index 000000000..4f3da5f30
--- /dev/null
+++ b/564/DEPENDENCIES/16_1data.sci
@@ -0,0 +1,5 @@
+a=200;//leangth of cross section along x-axis,given in mm
+b=300;//leangth of cross section along y-axis,given in mm
+tx=25;//thickness along x-axis,given in mm
+ty=20;//thickness along y-axis,given in mm
+M=-100*10^6;//bending moment,given in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_2data.sci b/564/DEPENDENCIES/16_2data.sci
new file mode 100755
index 000000000..97adbac5b
--- /dev/null
+++ b/564/DEPENDENCIES/16_2data.sci
@@ -0,0 +1,5 @@
+a=200;//leangth of cross section along x-axis,given in mm
+b=300;//leangth of cross section along y-axis,given in mm
+tx=25;//thickness along x-axis,given in mm
+ty=20;//thickness along y-axis,given in mm
+M=100*10^6;//bending moment,given in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_3data.sci b/564/DEPENDENCIES/16_3data.sci
new file mode 100755
index 000000000..9edce537c
--- /dev/null
+++ b/564/DEPENDENCIES/16_3data.sci
@@ -0,0 +1,7 @@
+a=200;//leangth of cross section along x-axis,given in mm
+b=300;//leangth of cross section along y-axis,given in mm
+tx=25;//thickness along x-axis,given in mm
+ty=20;//thickness along y-axis,given in mm
+M=100*10^6;//bending moment,given in N.mm
+th=30;//inclination of plane of bending moment,given in degrees
+theta=th*%pi/180;//th in radians
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_4data.sci b/564/DEPENDENCIES/16_4data.sci
new file mode 100755
index 000000000..5e355ca73
--- /dev/null
+++ b/564/DEPENDENCIES/16_4data.sci
@@ -0,0 +1,5 @@
+a=120;//leangth of cross section along x-axis,given in mm
+a1=40;//leangth AE
+b=88;//leangth of cross section along y-axis,given in mm
+t=8;//thickness,given in mm
+Mx=1500*10^3;//given in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_5data.sci b/564/DEPENDENCIES/16_5data.sci
new file mode 100755
index 000000000..c641d272f
--- /dev/null
+++ b/564/DEPENDENCIES/16_5data.sci
@@ -0,0 +1,3 @@
+W=10;//load in N
+L=20;//leanth in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_6data.sci b/564/DEPENDENCIES/16_6data.sci
new file mode 100755
index 000000000..c234798b0
--- /dev/null
+++ b/564/DEPENDENCIES/16_6data.sci
@@ -0,0 +1,3 @@
+W=5;//load in N/m
+L=20;//leanth in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_7data.sci b/564/DEPENDENCIES/16_7data.sci
new file mode 100755
index 000000000..c234798b0
--- /dev/null
+++ b/564/DEPENDENCIES/16_7data.sci
@@ -0,0 +1,3 @@
+W=5;//load in N/m
+L=20;//leanth in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_8data.sci b/564/DEPENDENCIES/16_8data.sci
new file mode 100755
index 000000000..c234798b0
--- /dev/null
+++ b/564/DEPENDENCIES/16_8data.sci
@@ -0,0 +1,3 @@
+W=5;//load in N/m
+L=20;//leanth in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/16_9data.sci b/564/DEPENDENCIES/16_9data.sci
new file mode 100755
index 000000000..8958bf7f7
--- /dev/null
+++ b/564/DEPENDENCIES/16_9data.sci
@@ -0,0 +1,3 @@
+W=10;//load in N/m
+a=5;//L/4, in m
+EI=70000;//flexural rigidity,in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/17_1data.sci b/564/DEPENDENCIES/17_1data.sci
new file mode 100755
index 000000000..ac2423937
--- /dev/null
+++ b/564/DEPENDENCIES/17_1data.sci
@@ -0,0 +1,3 @@
+h=100;//in mm
+t=2;//in mm
+Sy=200;//in N/mm^2
diff --git a/564/DEPENDENCIES/17_2data.sci b/564/DEPENDENCIES/17_2data.sci
new file mode 100755
index 000000000..0a19372e3
--- /dev/null
+++ b/564/DEPENDENCIES/17_2data.sci
@@ -0,0 +1,3 @@
+h=50;//in mm
+t=1.5;//in mm
+b=25;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/17_3data.sci b/564/DEPENDENCIES/17_3data.sci
new file mode 100755
index 000000000..7f2e0b89f
--- /dev/null
+++ b/564/DEPENDENCIES/17_3data.sci
@@ -0,0 +1 @@
+a=10;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/18_1data.sci b/564/DEPENDENCIES/18_1data.sci
new file mode 100755
index 000000000..9372b94a9
--- /dev/null
+++ b/564/DEPENDENCIES/18_1data.sci
@@ -0,0 +1,7 @@
+d=200;//given, in mm
+L=2000;//given,2m in mm
+T=30*10^6;//given, in N.mm
+Tmax=15*10^6;//given
+Smax=200;//maximum allowable shear stress,given,in N/mm^2
+angle=2*%pi/180;//maximum angle of twist,given,in degrees
+G=25000;//given, in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/18_2data.sci b/564/DEPENDENCIES/18_2data.sci
new file mode 100755
index 000000000..c7d898878
--- /dev/null
+++ b/564/DEPENDENCIES/18_2data.sci
@@ -0,0 +1,6 @@
+a=200;//in mm
+b=150;//in mm
+ta=2.5;//in mm
+tb=2;//in mm
+T=1000;//in N.mm
+G=25000;//given in N/mm^2
diff --git a/564/DEPENDENCIES/18_3data.sci b/564/DEPENDENCIES/18_3data.sci
new file mode 100755
index 000000000..4328ca669
--- /dev/null
+++ b/564/DEPENDENCIES/18_3data.sci
@@ -0,0 +1,6 @@
+a=25;//leangth along x-axis,given in mm
+b=50;//leangth along y-axis,given in mm
+tx=2.5;//in mm
+ty=1.5;//in mm
+T=10*10^3;//in N.mm
+G=25000;//given in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/19_1data.sci b/564/DEPENDENCIES/19_1data.sci
new file mode 100755
index 000000000..2a3d8aa24
--- /dev/null
+++ b/564/DEPENDENCIES/19_1data.sci
@@ -0,0 +1,5 @@
+L12=100;//leangth of 1-2,given in mm
+L45=200;//leangth of 4-5(3-4-5-6 is a square),given in mm
+t=2;//thickness,given in mm
+Sy=100*10^3;//in N
+G=25000;//given in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/19_2data.sci b/564/DEPENDENCIES/19_2data.sci
new file mode 100755
index 000000000..f8d398311
--- /dev/null
+++ b/564/DEPENDENCIES/19_2data.sci
@@ -0,0 +1,8 @@
+L12=900;//leangth of section 1-2(outer),given in mm
+L34=300;//leangth of section 3-4,given in mm
+L13=600;//leangth of section 1-3,given in mm
+t13=2;//thickness of section 1-3,given in mm
+t12=1.5;//thickness of section 1-2,given in mm
+T=10*10^6;//given in N.mm
+A=20000;//given in N.mm^2
+G=25000;//given in N.mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_01data.sci b/564/DEPENDENCIES/1_01data.sci
new file mode 100755
index 000000000..4e9b14913
--- /dev/null
+++ b/564/DEPENDENCIES/1_01data.sci
@@ -0,0 +1,6 @@
+d=2000;//given(2m in mm)
+t=20;//given(in mm)
+p=1.5;//pressure inside the vessel(N/mm^2)
+Load=2500*10^3;//axial tensile load(N)
+angle= %pi/3;//angle of plane with axis(radian)
+theta=%pi/2 -angle;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_02data.sci b/564/DEPENDENCIES/1_02data.sci
new file mode 100755
index 000000000..33e6a6333
--- /dev/null
+++ b/564/DEPENDENCIES/1_02data.sci
@@ -0,0 +1,6 @@
+t=1.5;//given(in mm)
+d=60;//given(in mm)
+T=1200*10^3;//in N.mm, given(in N.m)
+Load=50*10^3;//Compressive load(N)
+angle= %pi/3;//angle of plane with axis(radian)
+theta=%pi/2 - angle;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_03data.sci b/564/DEPENDENCIES/1_03data.sci
new file mode 100755
index 000000000..18d0be463
--- /dev/null
+++ b/564/DEPENDENCIES/1_03data.sci
@@ -0,0 +1,3 @@
+Sx=160;//σx, tension stress, given in N/mm^2
+Sy=-120;//σy, compression stress, given in N/mm^2
+lLoad=200;//limit load, given in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_04data.sci b/564/DEPENDENCIES/1_04data.sci
new file mode 100755
index 000000000..e349e669f
--- /dev/null
+++ b/564/DEPENDENCIES/1_04data.sci
@@ -0,0 +1,4 @@
+Sx=83;//σx, given (in N/mm^2)
+Sy=65;//σy, given (in N/mm^2)
+E=200000;// given (in N/mm^2)
+v=0.3;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_05data.sci b/564/DEPENDENCIES/1_05data.sci
new file mode 100755
index 000000000..b846e6421
--- /dev/null
+++ b/564/DEPENDENCIES/1_05data.sci
@@ -0,0 +1,5 @@
+Sx=60;//σx, given(in N/mm^2)
+Sy=-40;//-σy, given(in N/mm^2)
+Txy=50;//τxy,given(in N/mm^2)
+E=200000;//given(in N/mm^2)
+v=0.3;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/1_07data.sci b/564/DEPENDENCIES/1_07data.sci
new file mode 100755
index 000000000..cce15096a
--- /dev/null
+++ b/564/DEPENDENCIES/1_07data.sci
@@ -0,0 +1,6 @@
+d=50;//given (in mm)
+Ea=1000/10^6;
+Eb=-200/10^6;
+Ec=-300/10^6;
+v=0.3;
+E=70000;//given (in N/mm^2)
diff --git a/564/DEPENDENCIES/20_1data.sci b/564/DEPENDENCIES/20_1data.sci
new file mode 100755
index 000000000..4400a90b0
--- /dev/null
+++ b/564/DEPENDENCIES/20_1data.sci
@@ -0,0 +1,10 @@
+L16=400;//given in mm
+L34=200;//given in mm
+L12=600;//given in mm
+L23=600;//given in mm
+t12=2;//given in mm
+t23=1.5;//given in mm
+t34=2;//given in mm
+t25=2.5;//given in mm
+t16=3;//given in mm
+A=300;//given in mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/20_2data.sci b/564/DEPENDENCIES/20_2data.sci
new file mode 100755
index 000000000..b8f93bc61
--- /dev/null
+++ b/564/DEPENDENCIES/20_2data.sci
@@ -0,0 +1,5 @@
+Mx=100*10^6;
+B=[640;600;600;600;620;640;640;850;640];//matrix having the values of B
+D=[1200;1140;960;768;565;336;144;38;0];//matrix having the distence of all booms from boom 9
+
+
diff --git a/564/DEPENDENCIES/20_3data.sci b/564/DEPENDENCIES/20_3data.sci
new file mode 100755
index 000000000..a338bb644
--- /dev/null
+++ b/564/DEPENDENCIES/20_3data.sci
@@ -0,0 +1,3 @@
+Sy=4.8*10^3;//shear load,given,in N
+A=300;//Boom area,given,in mm^2
+L=200;//leangth,given,in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/20_4data.sci b/564/DEPENDENCIES/20_4data.sci
new file mode 100755
index 000000000..f08424b3a
--- /dev/null
+++ b/564/DEPENDENCIES/20_4data.sci
@@ -0,0 +1,4 @@
+Sy=10*10^3;//vertical load,given,in N
+B=[200;250;400;100;100;400;250;200];//boom area,given
+y=[30;100;100;50;-50;-100;-100;-30];//x-coordinate of each boom
+x=[120;240;240];//distences along x-axis between booms
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/20_5data.sci b/564/DEPENDENCIES/20_5data.sci
new file mode 100755
index 000000000..71b03545f
--- /dev/null
+++ b/564/DEPENDENCIES/20_5data.sci
@@ -0,0 +1,7 @@
+L=2000;//langth,given,in mm
+a=200;//channal dimension,given,in mm
+t=1;//thickness,given,in mm
+E=70000;//given,in N/mm^2
+G=30000;//given,in N/mm^2
+A=300;//boom area,given,in mm^2
+Sy=4.8*10^3;//given,in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/21_1data.sci b/564/DEPENDENCIES/21_1data.sci
new file mode 100755
index 000000000..16c0b06e7
--- /dev/null
+++ b/564/DEPENDENCIES/21_1data.sci
@@ -0,0 +1,7 @@
+Sy=-20*10^3;//given in N
+A=400;//given in mm^2
+t=2;//thickness,given in mm
+L=2*10^3;//leangth ,given in mm
+L1=1*10^3;//distence of section AA from the end
+We=200;//width at the end,given in mm
+W=400;//width at the joint,given in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/21_2data.sci b/564/DEPENDENCIES/21_2data.sci
new file mode 100755
index 000000000..b500fac15
--- /dev/null
+++ b/564/DEPENDENCIES/21_2data.sci
@@ -0,0 +1,11 @@
+Sy=100*10^3;//given in N
+Aco=900;//area of corner booms,given in mm^2
+Ace=1200;//area of central booms,given in mm^2
+tl=2;//thickness,given in mm
+tb=3;//thickness,given in mm
+L1=1.6*10^3;//leangth of built in end,given in mm
+L2=0.8*10^3;//leangth of open end,given in mm
+B1=0.8*10^3;//breadth of built in end,given in mm
+B2=0.4*10^3;//breadth of open end,given in mm
+Lc=2*10^3;//distence of the section from the built in end,given in mm
+L=4*10^3;//leangth of Beam,given in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/21_3data.sci b/564/DEPENDENCIES/21_3data.sci
new file mode 100755
index 000000000..b500fac15
--- /dev/null
+++ b/564/DEPENDENCIES/21_3data.sci
@@ -0,0 +1,11 @@
+Sy=100*10^3;//given in N
+Aco=900;//area of corner booms,given in mm^2
+Ace=1200;//area of central booms,given in mm^2
+tl=2;//thickness,given in mm
+tb=3;//thickness,given in mm
+L1=1.6*10^3;//leangth of built in end,given in mm
+L2=0.8*10^3;//leangth of open end,given in mm
+B1=0.8*10^3;//breadth of built in end,given in mm
+B2=0.4*10^3;//breadth of open end,given in mm
+Lc=2*10^3;//distence of the section from the built in end,given in mm
+L=4*10^3;//leangth of Beam,given in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/22_1data.sci b/564/DEPENDENCIES/22_1data.sci
new file mode 100755
index 000000000..996db5d0d
--- /dev/null
+++ b/564/DEPENDENCIES/22_1data.sci
@@ -0,0 +1,4 @@
+A=100;//cross sectional area of a stringer
+r=381;//radius of the fuselage,given in mm
+t=0.8;//thickness,given in mm
+Mx=200*10^6;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/22_2data.sci b/564/DEPENDENCIES/22_2data.sci
new file mode 100755
index 000000000..9103ea8b9
--- /dev/null
+++ b/564/DEPENDENCIES/22_2data.sci
@@ -0,0 +1,5 @@
+A=100;//cross sectional area of a stringer
+r=381;//radius of the fuselage,given in mm
+t=0.8;//thickness,given in mm
+Sy=100*10^3;//given in N
+d=150;//distence of action of load from center,given in mm
diff --git a/564/DEPENDENCIES/23_1data.sci b/564/DEPENDENCIES/23_1data.sci
new file mode 100755
index 000000000..69435ae1b
--- /dev/null
+++ b/564/DEPENDENCIES/23_1data.sci
@@ -0,0 +1,3 @@
+B=[2580;3880;3230;3230;3880;2580];//values of boom area,given in mm^2
+y=[165;230;200;-200;-230;-165];//values of y,given in mm
+Mx=300*10^6;//given in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/23_2data.sci b/564/DEPENDENCIES/23_2data.sci
new file mode 100755
index 000000000..5c96bd674
--- /dev/null
+++ b/564/DEPENDENCIES/23_2data.sci
@@ -0,0 +1,7 @@
+//matrices are in the order[120;12i;13;24;34;35;46;56]
+L=[1650;508;775;775;380;508;508;254];//length,given in mm
+t=[1.22;2.03;1.22;1.22;1.63;0.92;0.92;0.92];//thickness,given in mm
+G=[24200;27600;24200;24200;27600;20700;20700;20700];//given in N.mm^2
+Gref=27600;//given in N^mm^2
+A=[258000;355000;161000];//cell area,given in mm^2
+T=11.3*10^6;//torque,given in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/23_3data.sci b/564/DEPENDENCIES/23_3data.sci
new file mode 100755
index 000000000..49586b929
--- /dev/null
+++ b/564/DEPENDENCIES/23_3data.sci
@@ -0,0 +1,9 @@
+//in order 12,56,26,34,48,83,57,72,61,78
+L=[1023;1023;1274;2200;250;150;280;180;330;1270];//length
+t=[1.22;1.22;1.63;2.03;2.64;2.64;2.64;2.64;1.63;1.22];//thickness
+A=[265000;213000;413000];//cell area
+Gref=27600;
+G=[27600*ones(9,1);3*27600];
+B=[2580;3880;3230;3230;3880;2580];//boom areas
+y=[165;230;200;-200;-230;-165];
+Sy=86.8*10^3;//shear load
diff --git a/564/DEPENDENCIES/23_4data.sci b/564/DEPENDENCIES/23_4data.sci
new file mode 100755
index 000000000..1aac85227
--- /dev/null
+++ b/564/DEPENDENCIES/23_4data.sci
@@ -0,0 +1,12 @@
+B=[600;900;600;600;900;600];//matrix having the values of B
+t=[0.8;0.8;1;0.8;0.8;1;1];//thickness,in order-12,23,34,45,56,61,25;
+L11=200;//leangth 1 of built in end,given in mm
+L12=400;//leangth 2 of built in end,given in mm
+W1=180;//breadth of built in end,given in mm
+L21=150;//leangth 1 of open end,given in mm
+L22=300;//leangth 2 of open end,given in mm
+W2=80;//breadth of open end,given in mm
+L=1.2*10^3;//leangth of Beam,given in mm
+Mx=1.65*10^6;//given,in N.mm
+Sy=10*10^3;//given,in N
+G=27600;//in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/23_6data.sci b/564/DEPENDENCIES/23_6data.sci
new file mode 100755
index 000000000..e930cd1fd
--- /dev/null
+++ b/564/DEPENDENCIES/23_6data.sci
@@ -0,0 +1,4 @@
+T=10*10^6;//given, in N.mm
+L1=200;//Cross section dimention,given,in mm
+B1=800;//Cross section dimention,given,in mm
+L=1500;//in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/24_1data.sci b/564/DEPENDENCIES/24_1data.sci
new file mode 100755
index 000000000..e489dc8b7
--- /dev/null
+++ b/564/DEPENDENCIES/24_1data.sci
@@ -0,0 +1,8 @@
+GH=250;//given in mm
+DK=200;//given in mm
+KH=100;//given in mm
+FG=250;//given in mm
+EF=250;//given in mm
+L1=4000;//given in N
+L2=5000;//given in N
+theta=60*(%pi/180);//given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/24_2data.sci b/564/DEPENDENCIES/24_2data.sci
new file mode 100755
index 000000000..e9555e1e9
--- /dev/null
+++ b/564/DEPENDENCIES/24_2data.sci
@@ -0,0 +1,9 @@
+L11=300;
+L12=320;
+L1=12000;
+L2=15000;
+A1=50000;
+A2=95000;
+A21=49000;
+A22=A2-A21;
+theta=15*(%pi/180);
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_1data.sci b/564/DEPENDENCIES/25_1data.sci
new file mode 100755
index 000000000..2730d4fd4
--- /dev/null
+++ b/564/DEPENDENCIES/25_1data.sci
@@ -0,0 +1,10 @@
+Ee=5000;//modulus of epoxy,given,in N/mm^2
+Ec=200000;//modulus of carbon,given,in N/mm^2
+ve=0.2;//poission's ratio of epoxy
+vc=0.3;//poission's ratio of carbon
+L=500;//langth of bar,given,in mm
+L1=80;//langth of cross section ,given,in mm
+Be1=20;//width of first section of epoxy,given,in mm
+Be2=20;//width of other section of epoxy,given,in mm
+Bc=10;//width of section of carbon,given,in mm
+Load=100*10^3;//axial tensile load,given,in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_2data.sci b/564/DEPENDENCIES/25_2data.sci
new file mode 100755
index 000000000..6d78834a8
--- /dev/null
+++ b/564/DEPENDENCIES/25_2data.sci
@@ -0,0 +1,7 @@
+S1=50;//longitudnal direct stress,given,in N/mm^2
+S2=25;//transverse direct stress,given,in N/mm^2
+S3=40;//shear stress,given,in N/mm^2
+El=120000;//given,in N/mm^2
+Et=80000;//given,in N/mm^2
+Glt=5000;//given,in N/mm^2
+vlt=0.3;//given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_3data.sci b/564/DEPENDENCIES/25_3data.sci
new file mode 100755
index 000000000..a16697c14
--- /dev/null
+++ b/564/DEPENDENCIES/25_3data.sci
@@ -0,0 +1,7 @@
+S1=60;//stress parallel to the x reference axis,given,in N/mm^2
+S2=40;//stress perpendicular to the x reference axis,given,in N/mm^2
+theta=45*(%pi/180);//given
+El=150000;//given,in N/mm^2
+Et=90000;//given,in N/mm^2
+Glt=5000;//given,in N/mm^2
+vlt=0.3;//given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_4data.sci b/564/DEPENDENCIES/25_4data.sci
new file mode 100755
index 000000000..cbb2e1107
--- /dev/null
+++ b/564/DEPENDENCIES/25_4data.sci
@@ -0,0 +1,7 @@
+a=100;
+b=150;
+ta=1;
+tb=2;
+Ez1=60000;
+Ez2=20000;
+Load=40*10^3;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_5data.sci b/564/DEPENDENCIES/25_5data.sci
new file mode 100755
index 000000000..0a0dad43e
--- /dev/null
+++ b/564/DEPENDENCIES/25_5data.sci
@@ -0,0 +1,7 @@
+a=50;//cross section dimension,given,in mm
+b=100;//cross section dimension,given,in mm
+ta=1;//thickness,given,in mm
+tb=2;//thickness,given,in mm
+Ez1=50000;//Young’s modulus for the flange,given,in N/mm^2
+Ez2=15000;//Young’s modulus for the web,given,in N/mm^2
+Mx=1*10^6;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_6data.sci b/564/DEPENDENCIES/25_6data.sci
new file mode 100755
index 000000000..8bf002b48
--- /dev/null
+++ b/564/DEPENDENCIES/25_6data.sci
@@ -0,0 +1,7 @@
+E23=20000;//laminate Young’s modulus of 23,given,in N/mm^2
+E12=45000;//laminate Young’s modulus of 12,given,in N/mm^2
+L12=250;//langth of 12,given,in mm
+L23=300;//langth of 23,given,in mm
+t12=2;//thickness of 12,given,in mm
+t23=1.5;//thickness of 23,given,in mm
+Sy=2*10^3;//load
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_7data.sci b/564/DEPENDENCIES/25_7data.sci
new file mode 100755
index 000000000..9f80d9851
--- /dev/null
+++ b/564/DEPENDENCIES/25_7data.sci
@@ -0,0 +1,7 @@
+T=10*10^6;//Torque,given,in N.mm
+l=100;//given,in mm
+b=200;//given,in mm
+t1=2;//given,in mm
+t2=1;//given,in mm
+Gw=35000;//given,in N/mm^2
+Gc=20000;//given,in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/25_8data.sci b/564/DEPENDENCIES/25_8data.sci
new file mode 100755
index 000000000..b319f607a
--- /dev/null
+++ b/564/DEPENDENCIES/25_8data.sci
@@ -0,0 +1,7 @@
+T=10*10^3;//given,in N.mm
+G1=20000;//laminate shear modulus of flanges,given,in N/mm^2
+G2=15000;//laminate shear modulus of web,given,in N/mm^2
+a=25;
+b=50;
+t1=1.5;
+t2=2.5;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/26_1data.sci b/564/DEPENDENCIES/26_1data.sci
new file mode 100755
index 000000000..5cb35cd0d
--- /dev/null
+++ b/564/DEPENDENCIES/26_1data.sci
@@ -0,0 +1,9 @@
+L12=375;//given,in mm
+L23=500;//given,in mm
+L34=125;//given,in mm
+t12=1.6;//given,in mm
+t23=1;//given,in mm
+t34=1.2;//given,in mm
+x=100;//given,in mm
+Load=22000;//given,in mm
+G=70000;//in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/26_2data.sci b/564/DEPENDENCIES/26_2data.sci
new file mode 100755
index 000000000..dc24cd5c3
--- /dev/null
+++ b/564/DEPENDENCIES/26_2data.sci
@@ -0,0 +1,4 @@
+R=200;//given,in mm
+Load=20000;//given,in N
+G=70000;//in N/mm^2
+t=2;//thickness in mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/26_3data.sci b/564/DEPENDENCIES/26_3data.sci
new file mode 100755
index 000000000..90f42987e
--- /dev/null
+++ b/564/DEPENDENCIES/26_3data.sci
@@ -0,0 +1,8 @@
+L1=5000;//given,in mm
+A=800;//booms area,given,in mm^2
+G=20000;//given,in N/mm^2
+E=G/0.36;//given
+T=20;//torque loading,given,in N.m/mm
+t=1;//wall thickness,given,in mm
+a=500;
+b=200;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/26_4data.sci b/564/DEPENDENCIES/26_4data.sci
new file mode 100755
index 000000000..ae97a0d1c
--- /dev/null
+++ b/564/DEPENDENCIES/26_4data.sci
@@ -0,0 +1,7 @@
+t1=2;
+t2=3;
+a=100;
+b=600;
+L=2000;//length
+Load=20*10^3;//given in N
+GE=0.36;//  G/E
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/27_2data.sci b/564/DEPENDENCIES/27_2data.sci
new file mode 100755
index 000000000..40a20eb75
--- /dev/null
+++ b/564/DEPENDENCIES/27_2data.sci
@@ -0,0 +1,7 @@
+a=100;//cross section dimention,given,in mm
+b=200;//cross section dimention,given,in mm
+t=5;//thickness,given,in mm
+L=3000;//length,given,in mm
+E=200000;//given,in N/mm^2
+G=0.36*200000;//given,in N/mm^2
+P=-100*10^3;//load,given,in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/3_1data.sci b/564/DEPENDENCIES/3_1data.sci
new file mode 100755
index 000000000..4ac8460f8
--- /dev/null
+++ b/564/DEPENDENCIES/3_1data.sci
@@ -0,0 +1,3 @@
+R=16;//radius of bar,in mm
+T=100000;//applied torque,in N.mm
+G=76923;//shear modulus,in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/3_2data.sci b/564/DEPENDENCIES/3_2data.sci
new file mode 100755
index 000000000..0e2bb9a3d
--- /dev/null
+++ b/564/DEPENDENCIES/3_2data.sci
@@ -0,0 +1,4 @@
+a=4;//major axis,in mm
+b=3;//minor axis of bar,in mm
+T=100000;//applied torque,in N.mm
+G=76923;//shear modulus,in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_1data.sci b/564/DEPENDENCIES/4_1data.sci
new file mode 100755
index 000000000..fbfdc14ce
--- /dev/null
+++ b/564/DEPENDENCIES/4_1data.sci
@@ -0,0 +1,3 @@
+W=100;//in N
+L=20;//in m
+a=14;//in m
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_2.sce b/564/DEPENDENCIES/4_2.sce
new file mode 100755
index 000000000..dbcce4e4f
--- /dev/null
+++ b/564/DEPENDENCIES/4_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('4.2.sce')
+filename=pathname+filesep()+'4.2data.sci'
+exec(filename)
+M=(W*a*(L-a))/(L);
+deff("[y]=f(x)","y=(W*x*(L-x))/(L)")//manding moment
+x=[0:0.05:L];
+fplot2d(x,f);
+xgrid(3);
+datatipToggle();
+xtitle( 'Banding Moment versus a', ' -a- ', '-M-');
+printf("\nMB: %f N.m",M);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_2data.sci b/564/DEPENDENCIES/4_2data.sci
new file mode 100755
index 000000000..bd83d6b6b
--- /dev/null
+++ b/564/DEPENDENCIES/4_2data.sci
@@ -0,0 +1,3 @@
+W=90;//in Newton
+a=10;//in m
+L=15;//(a+b),in m
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_3data.sci b/564/DEPENDENCIES/4_3data.sci
new file mode 100755
index 000000000..d298ea740
--- /dev/null
+++ b/564/DEPENDENCIES/4_3data.sci
@@ -0,0 +1,3 @@
+CD=4;//given in m
+BD=3;//given in m
+Loadc=30;//given in KN
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_4data.sci b/564/DEPENDENCIES/4_4data.sci
new file mode 100755
index 000000000..3d1ff1584
--- /dev/null
+++ b/564/DEPENDENCIES/4_4data.sci
@@ -0,0 +1,4 @@
+w=10;//intensity of distributed load, in N/m
+L=10;//leangth in m
+E=200000;//in N/mm^2
+I=0.5;// moment of Inertia of cross section, in m^4
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_5data.sci b/564/DEPENDENCIES/4_5data.sci
new file mode 100755
index 000000000..c7b2ebf64
--- /dev/null
+++ b/564/DEPENDENCIES/4_5data.sci
@@ -0,0 +1,4 @@
+W=100;//in N
+L=20;//in m
+E=200000;//in N/mm^2
+I=0.5// moment of Inertia of cross section, in m^4
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/4_6data.sci b/564/DEPENDENCIES/4_6data.sci
new file mode 100755
index 000000000..da5750c4a
--- /dev/null
+++ b/564/DEPENDENCIES/4_6data.sci
@@ -0,0 +1,8 @@
+CS=1800;//cross-section area,given in mm^2
+E=200000;//given in N/mm^2
+h=4;//height, given in m
+AB=4;//given in m
+BC=4;//given in m
+CD=4;//given in m
+Load1=40*10^3;//load at point E,given in N
+Load2=100*10^3;//load at point C,given in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_10data.sci b/564/DEPENDENCIES/5_10data.sci
new file mode 100755
index 000000000..55806b80e
--- /dev/null
+++ b/564/DEPENDENCIES/5_10data.sci
@@ -0,0 +1,6 @@
+a=4000;//langth,given,in mm
+b=3000;//breadth,given,in mm
+T=30;//change in temprature,given,in degrees
+A=200;//given,in mm^2
+alpha=7/10^6;//given,in 1/C;
+E=200000;//given,in N/mm^2
diff --git a/564/DEPENDENCIES/5_11data.sci b/564/DEPENDENCIES/5_11data.sci
new file mode 100755
index 000000000..cc7f33623
--- /dev/null
+++ b/564/DEPENDENCIES/5_11data.sci
@@ -0,0 +1,3 @@
+EI=70000;//in N.mm^2
+W=10;//in N
+L=20;//in m
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_12data.sci b/564/DEPENDENCIES/5_12data.sci
new file mode 100755
index 000000000..acc2070cc
--- /dev/null
+++ b/564/DEPENDENCIES/5_12data.sci
@@ -0,0 +1,5 @@
+D=[0;100;200;300;400;500;600;700;800];
+V=[0;-0.3;-1.4;-2.5;-1.9;0;2.3;4.8;10.6];
+L0=40;
+L1=30;
+L2=10;
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_1data.sci b/564/DEPENDENCIES/5_1data.sci
new file mode 100755
index 000000000..da5750c4a
--- /dev/null
+++ b/564/DEPENDENCIES/5_1data.sci
@@ -0,0 +1,8 @@
+CS=1800;//cross-section area,given in mm^2
+E=200000;//given in N/mm^2
+h=4;//height, given in m
+AB=4;//given in m
+BC=4;//given in m
+CD=4;//given in m
+Load1=40*10^3;//load at point E,given in N
+Load2=100*10^3;//load at point C,given in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_2data.sci b/564/DEPENDENCIES/5_2data.sci
new file mode 100755
index 000000000..9c011b944
--- /dev/null
+++ b/564/DEPENDENCIES/5_2data.sci
@@ -0,0 +1,4 @@
+clear all
+w=5;//in N/m
+L=20;//in m
+EI=70000;//in N.m^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_3data.sci b/564/DEPENDENCIES/5_3data.sci
new file mode 100755
index 000000000..058ae4be0
--- /dev/null
+++ b/564/DEPENDENCIES/5_3data.sci
@@ -0,0 +1,6 @@
+AC=30;//C.S. area of AC,given,in mm^2
+A=20;//C.S. areaa of other elements,given,in mm^2
+L=800;//langth,given in mm
+E=200000;//given in N/mm^2
+theta=60*(%pi/180);
+Load=100;//given,in N
diff --git a/564/DEPENDENCIES/5_4data.sci b/564/DEPENDENCIES/5_4data.sci
new file mode 100755
index 000000000..58db9bfb4
--- /dev/null
+++ b/564/DEPENDENCIES/5_4data.sci
@@ -0,0 +1,7 @@
+clear all
+AB=4000;//langth AB,given,in mm
+BC=3000;//langth BC,given,in mm
+A=200;//area of bar,given,in mm^2
+delT=30;//change in temprature,given,in celcius;
+alpha=7/10^6;//coefficient of linear expansion,given,in 1/celcius
+E=200000;//given in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_5data.sci b/564/DEPENDENCIES/5_5data.sci
new file mode 100755
index 000000000..0d4569d42
--- /dev/null
+++ b/564/DEPENDENCIES/5_5data.sci
@@ -0,0 +1,7 @@
+L=2000;//langth.in mm
+Ab=200;//in mm^2
+A=180;//area of bar,in mm^2
+P=1000;//in N
+theta=60*(%pi/180);
+E=200000;//in N/mm^2
+I=100000;//in mm^4
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_6data.sci b/564/DEPENDENCIES/5_6data.sci
new file mode 100755
index 000000000..16e25e8ad
--- /dev/null
+++ b/564/DEPENDENCIES/5_6data.sci
@@ -0,0 +1,2 @@
+r=1000;//in mm
+M0=1000;//in N.mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_7data.sci b/564/DEPENDENCIES/5_7data.sci
new file mode 100755
index 000000000..0ce9fcd0e
--- /dev/null
+++ b/564/DEPENDENCIES/5_7data.sci
@@ -0,0 +1,4 @@
+L=2000;//langth.in mm
+EI=2*10^10;
+GJ=5*10^10;
+W=0.1;//in N/mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_8data.sci b/564/DEPENDENCIES/5_8data.sci
new file mode 100755
index 000000000..92cf4da60
--- /dev/null
+++ b/564/DEPENDENCIES/5_8data.sci
@@ -0,0 +1,2 @@
+L=1000;//langth in mm
+P=1;//Load in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/5_9data.sci b/564/DEPENDENCIES/5_9data.sci
new file mode 100755
index 000000000..1956f0af7
--- /dev/null
+++ b/564/DEPENDENCIES/5_9data.sci
@@ -0,0 +1,2 @@
+L=1;//000;//langth,given,in mm
+P=10;//*10^3;//load,given,in N
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/6_1data.sci b/564/DEPENDENCIES/6_1data.sci
new file mode 100755
index 000000000..c1a4bd1f0
--- /dev/null
+++ b/564/DEPENDENCIES/6_1data.sci
@@ -0,0 +1,5 @@
+L12=10;//in m
+L13=10;//in m
+AE=200000;//in N
+FX2=0;//FORCE ON POINT 2 IN X DIRECTION
+FY2=-10;//-W,FORCE ON POINT 2 IN Y DIRECTION
diff --git a/564/DEPENDENCIES/6_2data.sci b/564/DEPENDENCIES/6_2data.sci
new file mode 100755
index 000000000..63b8d16e1
--- /dev/null
+++ b/564/DEPENDENCIES/6_2data.sci
@@ -0,0 +1,4 @@
+L=10;//in m
+W=1000;//FORCE ON POINT 2
+M=100;//MOMENT ON POINT 2 
+EI=20000;
diff --git a/564/DEPENDENCIES/6_4data.sci b/564/DEPENDENCIES/6_4data.sci
new file mode 100755
index 000000000..80fa904be
--- /dev/null
+++ b/564/DEPENDENCIES/6_4data.sci
@@ -0,0 +1,8 @@
+P1=[-2;-1];
+P2=[2;-1];
+P3=[2;1];
+P4=[-2;1];
+u=[0.001;0.003;-0.003;0];//given in m
+v=[-0.004;-0.002;0.001;0.001];//given in m
+E=200000;//given in N/mm^2
+V=0.3;//poission's ratio,given
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/7_1data.sci b/564/DEPENDENCIES/7_1data.sci
new file mode 100755
index 000000000..e18b04467
--- /dev/null
+++ b/564/DEPENDENCIES/7_1data.sci
@@ -0,0 +1,6 @@
+t=2;//in mm
+E=70000;//in N/mm^2
+v=0.3;
+a=100;//in mm
+b=100;//in mm
+q0=10;//in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/7_3data.sci b/564/DEPENDENCIES/7_3data.sci
new file mode 100755
index 000000000..0c6581500
--- /dev/null
+++ b/564/DEPENDENCIES/7_3data.sci
@@ -0,0 +1,7 @@
+t=2;//in mm
+E=70000;//in N/mm^2
+v=0.3;
+a=100;//in mm
+b=100;//in mm
+q0=10;//in N/mm^2
+Nx=10;//in n/mm
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/7_4data.sci b/564/DEPENDENCIES/7_4data.sci
new file mode 100755
index 000000000..b55fdc564
--- /dev/null
+++ b/564/DEPENDENCIES/7_4data.sci
@@ -0,0 +1,6 @@
+a=100;//in mm
+b=80;//in mm
+t=2;//in mm
+E=70000;//in N/mm^2
+v=0.3;
+q0=10;//in N/mm^2
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/8_2data.sci b/564/DEPENDENCIES/8_2data.sci
new file mode 100755
index 000000000..d98d24bbb
--- /dev/null
+++ b/564/DEPENDENCIES/8_2data.sci
@@ -0,0 +1,5 @@
+L=20;//in mm
+P=100;//in N
+e=0.1;//in m
+E=200000;//in N/^2
+I=0.5// moment of Inertia of cross section, in m^4
+\ No newline at end of file
diff --git a/564/DEPENDENCIES/8_3data.sci b/564/DEPENDENCIES/8_3data.sci
new file mode 100755
index 000000000..05cfe4c3e
--- /dev/null
+++ b/564/DEPENDENCIES/8_3data.sci
@@ -0,0 +1,6 @@
+L=2000;//given (2m in mm)
+a=37.5;//width of cross section(along x axis),given in mm
+b=75;//leanght of cross section(along y axis),given in mm
+t=2.5;//thickness
+E=75000;//given in N/mm^2
+G=21000;//given in N/mm^2
diff --git a/564/DEPENDENCIES/8_4data.sci b/564/DEPENDENCIES/8_4data.sci
new file mode 100755
index 000000000..5a4f5d572
--- /dev/null
+++ b/564/DEPENDENCIES/8_4data.sci
@@ -0,0 +1,6 @@
+L=1000;//given (1m in mm)
+a2=100;//width of cross section(along x axis),given in mm
+b=100;//leanght of cross section(along y axis),given in mm
+t=2;//thickness
+E=70000;//given in N/mm^2
+G=30000;//given in N/mm^2
diff --git a/564/DEPENDENCIES/9_1data.sci b/564/DEPENDENCIES/9_1data.sci
new file mode 100755
index 000000000..45b5c06a6
--- /dev/null
+++ b/564/DEPENDENCIES/9_1data.sci
@@ -0,0 +1,11 @@
+Af=350;//cross sectional area of flanges,given in mm^2
+As=300;//Cross sectional area of stiffners,given in mm^2
+ESM=750;//elastic section modulus,given in mm^3
+t=2;//thickness of web,given in mm
+MA=2000;//2nd moment of area of a stiffner,given in mm^4
+E=70000;//given in N/mm^2
+I=2000;//moment of inertia
+d=400;//given in mm
+b=300;//given in mm
+W=5000;//given (5 KN in N)
+z=1200;//given in mm
+\ No newline at end of file