initial commit / add all books

author: priyanka 2015-06-24 15:03:17 +0530
committer: priyanka 2015-06-24 15:03:17 +0530
commit: b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
tree: ab291cffc65280e58ac82470ba63fbcca7805165 /2021
download: Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
324 files changed, 2408 insertions, 0 deletions
diff --git a/2021/CH1/EX1.1/Ex1_1.pdf b/2021/CH1/EX1.1/Ex1_1.pdf
new file mode 100755
index 000000000..da0c3864f
--- /dev/null
+++ b/2021/CH1/EX1.1/Ex1_1.pdf
diff --git a/2021/CH1/EX1.1/Ex1_1.sce b/2021/CH1/EX1.1/Ex1_1.sce
new file mode 100755
index 000000000..d6f516147
--- /dev/null
+++ b/2021/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,12 @@
+// Finding Specific weight,Density,Specific Gravity
+//Given
+V=1/1000;                    //volume in m^3
+w=9.6;                  //weight in Newton
+g=9.81;                 //gravitational force in m/s^2
+//To Find
+spwt=w/V;               //Specific weoght in N/m^3
+rho=spwt/g;                //density in kg/m^3
+spgr=rho/1000;          //Specific gravity no units
+disp("specific weight = "+string(spwt)+" N/m^3");
+disp("density = "+string(rho)+" kg/m^3");
+disp("specific gravity = "+string(spgr)+" no unit");
diff --git a/2021/CH1/EX1.3/EX1_3.sce b/2021/CH1/EX1.3/EX1_3.sce
new file mode 100755
index 000000000..94dbac2c2
--- /dev/null
+++ b/2021/CH1/EX1.3/EX1_3.sce
@@ -0,0 +1,10 @@
+//Finding of Viscosity
+//Given
+dy=0.025*10^-3;            //distance in meter
+du=0.5;                    //velocity  in m/s 
+tau=1.471;                 //shear stress in N/m^2
+//To Find
+mu=tau*dy/du;              //viscosity in Ns/m^2 
+mu1=mu*10;                 // Viscosity in poise
+disp("viscosity ="+string(mu)+" in Ns/m^2");
+disp("Viscosity ="+string(mu1)+" in poise")  
diff --git a/2021/CH1/EX1.3/Ex1_3.pdf b/2021/CH1/EX1.3/Ex1_3.pdf
new file mode 100755
index 000000000..07c234a9a
--- /dev/null
+++ b/2021/CH1/EX1.3/Ex1_3.pdf
diff --git a/2021/CH1/EX1.4/EX1_4.sce b/2021/CH1/EX1.4/EX1_4.sce
new file mode 100755
index 000000000..5d871b424
--- /dev/null
+++ b/2021/CH1/EX1.4/EX1_4.sce
@@ -0,0 +1,11 @@
+//Finding of Diameter of water droplet
+//Given
+st=0.716;                   //Surface Tension in N/m
+p=0.147*10^4;               //Pressure in N/m^2
+//To Find 
+d=4*st/p;                   //Dia in meter 
+d1=d*10^2;                  //Dia in centimeter 
+d2=d*10^3;                  // Dia in millimeter
+disp("Diameter ="+string(d)+" meter");
+disp("Diameter ="+string(d1)+" Centi meter");
+disp("Diameter ="+string(d2)+" Milli meter");
diff --git a/2021/CH1/EX1.4/Ex1_4.pdf b/2021/CH1/EX1.4/Ex1_4.pdf
new file mode 100755
index 000000000..2af5131b5
--- /dev/null
+++ b/2021/CH1/EX1.4/Ex1_4.pdf
diff --git a/2021/CH1/EX1.5/EX1_5.sce b/2021/CH1/EX1.5/EX1_5.sce
new file mode 100755
index 000000000..3ac1df2af
--- /dev/null
+++ b/2021/CH1/EX1.5/EX1_5.sce
@@ -0,0 +1,10 @@
+//Finding of Shear Stress
+//Given
+//du/dy = vg
+vg=.25;            //Velocity gradient in m/sec/meter
+nu=6.30*10^-4;         //Kinematic viscosity in m^2/sec
+rho=1268.4;             //Mass density in Kg/m^3
+mu=rho*nu;              //Dynamic Viscosity
+//To Find
+tau=mu*vg;          //Shear stress in N/m^2
+disp("Shear Stress ="+string(tau)+" N/m^2");
diff --git a/2021/CH1/EX1.5/Ex1_5.pdf b/2021/CH1/EX1.5/Ex1_5.pdf
new file mode 100755
index 000000000..af5003c76
--- /dev/null
+++ b/2021/CH1/EX1.5/Ex1_5.pdf
diff --git a/2021/CH1/EX1.6/EX1_6.sce b/2021/CH1/EX1.6/EX1_6.sce
new file mode 100755
index 000000000..e23be7ef2
--- /dev/null
+++ b/2021/CH1/EX1.6/EX1_6.sce
@@ -0,0 +1,7 @@
+//Finding of increase of Pressure
+//Given
+k=2.07*10^6;             // Bulk Modulus in KN/m^2
+dv=0.01;              //Change in Volume  
+//To Find
+p=k*(dv);            // Change in pressure
+disp(" Increase in Pressure ="+string(p)+" KN/m^2");
diff --git a/2021/CH1/EX1.6/Ex1_6.pdf b/2021/CH1/EX1.6/Ex1_6.pdf
new file mode 100755
index 000000000..239a86b4b
--- /dev/null
+++ b/2021/CH1/EX1.6/Ex1_6.pdf
diff --git a/2021/CH1/EX1.7/EX1_7.pdf b/2021/CH1/EX1.7/EX1_7.pdf
new file mode 100755
index 000000000..2c0f79030
--- /dev/null
+++ b/2021/CH1/EX1.7/EX1_7.pdf
diff --git a/2021/CH1/EX1.7/EX1_7.sce b/2021/CH1/EX1.7/EX1_7.sce
new file mode 100755
index 000000000..d29a1e4a9
--- /dev/null
+++ b/2021/CH1/EX1.7/EX1_7.sce
@@ -0,0 +1,12 @@
+//Finding of Cappilary rise
+//Given
+d=0.03*10^-2;                //Diameter in meter
+st=0.0735;                   //Surface Tension in N/m
+x=0;                     //contact angle in degree
+w=1000*9.81;
+//To Find
+h=(4*st)*cos(x)/(w*d);
+h1=h*10^2;
+disp(h);
+disp("Capillary rise ="+string(h)+" meter");
+disp("Capillary rise ="+string(h1)+" Centi meter");
diff --git a/2021/CH1/EX1.8/EX1_8.sce b/2021/CH1/EX1.8/EX1_8.sce
new file mode 100755
index 000000000..d34428c4c
--- /dev/null
+++ b/2021/CH1/EX1.8/EX1_8.sce
@@ -0,0 +1,14 @@
+//Finding of Kinematic Viscosity
+//Given
+tau=0.2158;                 //Shear stress in N/m^2
+vg=0.218;                //Velocity Gradient in sec^-1
+rho=959.5;                  //Density in Kg/m^3;
+//To Find 
+mu=tau*1/vg;
+disp("Dynamic Viscosity ="+string(mu)+" Ns/m^2");
+nu=mu/rho;
+disp("Kinematic Viscosity ="+string(nu)+" m^2/sec");
+nu1=nu*10^4;
+disp(" Knematic Viscosity ="+string(nu1)+" cm^2/sec");
+nu2=nu1*10^-4;
+disp("Kinematic Viscosity ="+string(nu2)+" strokes");
diff --git a/2021/CH1/EX1.8/Ex1_8.pdf b/2021/CH1/EX1.8/Ex1_8.pdf
new file mode 100755
index 000000000..fa4b180c5
--- /dev/null
+++ b/2021/CH1/EX1.8/Ex1_8.pdf
diff --git a/2021/CH10/EX10.1/EX10_1.pdf b/2021/CH10/EX10.1/EX10_1.pdf
new file mode 100755
index 000000000..f05fa26fe
--- /dev/null
+++ b/2021/CH10/EX10.1/EX10_1.pdf
diff --git a/2021/CH10/EX10.1/EX10_1.sce b/2021/CH10/EX10.1/EX10_1.sce
new file mode 100755
index 000000000..3770ea895
--- /dev/null
+++ b/2021/CH10/EX10.1/EX10_1.sce
@@ -0,0 +1,12 @@
+//Finding of velocity of flow and discharge
+//Given
+c=50;
+sb=1/3000;
+R=10/9;
+a=10;
+//To Find
+b=R*sb;
+v=c*sqrt(b);
+q=a*v;
+disp("Velocity of flow ="+string(v)+" m/sec");
+disp("Discharge ="+string(q)+" m^3/sec");
diff --git a/2021/CH10/EX10.11/EX10_11.pdf b/2021/CH10/EX10.11/EX10_11.pdf
new file mode 100755
index 000000000..03dd19c63
--- /dev/null
+++ b/2021/CH10/EX10.11/EX10_11.pdf
diff --git a/2021/CH10/EX10.11/EX10_11.sce b/2021/CH10/EX10.11/EX10_11.sce
new file mode 100755
index 000000000..c6d46e3a5
--- /dev/null
+++ b/2021/CH10/EX10.11/EX10_11.sce
@@ -0,0 +1,13 @@
+//Finding of Slope 
+//Given
+z=1;
+y=0.225;
+c=50;
+q=0.04;
+//To Find
+A=z*y^2;
+P=2*sqrt(2)*y;
+x=sqrt(0.225/(2*sqrt(2)));
+sb=q/(A*c*x);
+sb1=sb^(2);
+disp("Slope ="+string(sb1)+" no units");
diff --git a/2021/CH10/EX10.12/EX10_12.pdf b/2021/CH10/EX10.12/EX10_12.pdf
new file mode 100755
index 000000000..369b66a8f
--- /dev/null
+++ b/2021/CH10/EX10.12/EX10_12.pdf
diff --git a/2021/CH10/EX10.12/EX10_12.sce b/2021/CH10/EX10.12/EX10_12.sce
new file mode 100755
index 000000000..463ea4730
--- /dev/null
+++ b/2021/CH10/EX10.12/EX10_12.sce
@@ -0,0 +1,14 @@
+//Finding of C and f
+//Given
+n=0.012;
+d=0.5;
+w=2;
+g=9.81;
+//To Find
+A=w*d;
+P=2+(w*d);
+R=P/A;
+C=(1/n)*(R)^(1/6);
+f=sqrt((8*g)/(C^2));
+disp(" C="+string(C)+" m/sec");
+disp("f ="+string(f)+" no units");
diff --git a/2021/CH10/EX10.13/EX10_13.pdf b/2021/CH10/EX10.13/EX10_13.pdf
new file mode 100755
index 000000000..0a3bf5c8f
--- /dev/null
+++ b/2021/CH10/EX10.13/EX10_13.pdf
diff --git a/2021/CH10/EX10.13/EX10_13.sce b/2021/CH10/EX10.13/EX10_13.sce
new file mode 100755
index 000000000..3ba3e5e5b
--- /dev/null
+++ b/2021/CH10/EX10.13/EX10_13.sce
@@ -0,0 +1,12 @@
+//Finding of Normal Depth
+//Given
+w=6;
+q=5;
+sb=0.006;
+n=0.014;
+B=6;
+//To Find
+a=(q/(B^(8/3)*sb^(1/2)))^(3/5);
+b=(1+(0.855)*((q/B^(8/3)*sb^(1/2)))^(3/5));
+y=a*b;
+disp("Normal Depth ="+string(y)+" meter");
diff --git a/2021/CH10/EX10.14/EX10_14.pdf b/2021/CH10/EX10.14/EX10_14.pdf
new file mode 100755
index 000000000..7c7b65c86
--- /dev/null
+++ b/2021/CH10/EX10.14/EX10_14.pdf
diff --git a/2021/CH10/EX10.14/EX10_14.sce b/2021/CH10/EX10.14/EX10_14.sce
new file mode 100755
index 000000000..47ca84b56
--- /dev/null
+++ b/2021/CH10/EX10.14/EX10_14.sce
@@ -0,0 +1,15 @@
+//Finding of velocity , Dischage
+//Given
+z=1.5;
+sb=0.0003;
+B=10;
+n=0.012;
+y=3;
+//To Find
+A=(B+(z*y))*y;
+P=B+(2*y)*sqrt(1+z^2);
+R=A/P;
+v=(1/n)*R^(2/3)*sb^(1/2);
+q=A*v;
+disp("Velocity ="+string(v)+" m/sec^2");
+disp("Discharge ="+string(q)+" m^3/sec");
diff --git a/2021/CH10/EX10.16/EX10_16.pdf b/2021/CH10/EX10.16/EX10_16.pdf
new file mode 100755
index 000000000..ea6ab71ec
--- /dev/null
+++ b/2021/CH10/EX10.16/EX10_16.pdf
diff --git a/2021/CH10/EX10.16/EX10_16.sce b/2021/CH10/EX10.16/EX10_16.sce
new file mode 100755
index 000000000..8e8b4d168
--- /dev/null
+++ b/2021/CH10/EX10.16/EX10_16.sce
@@ -0,0 +1,12 @@
+//Finding of specific energy
+//Given
+B=4;
+y=2.5;
+q=8;
+g=9.81;
+//To Find
+A=B*y;
+v=q/A;
+E=y+(v^2/(2*g));
+disp(" Specific Energy ="+string(E)+" meter of water");
+
diff --git a/2021/CH10/EX10.17/EX10_17.pdf b/2021/CH10/EX10.17/EX10_17.pdf
new file mode 100755
index 000000000..22756c532
--- /dev/null
+++ b/2021/CH10/EX10.17/EX10_17.pdf
diff --git a/2021/CH10/EX10.17/EX10_17.sce b/2021/CH10/EX10.17/EX10_17.sce
new file mode 100755
index 000000000..763153c3f
--- /dev/null
+++ b/2021/CH10/EX10.17/EX10_17.sce
@@ -0,0 +1,13 @@
+//Finding of Critical depth,velocity,Minimum Specific energy
+//Given
+Q=18;
+B=6;
+q=Q/B;
+g=9.81;
+//To Find
+y=(q^2/g)^(1/3);
+v=q/y;
+E=(3/2)*y;
+disp(" Critical depth ="+string(y)+" meter");
+disp(" Critical velocity ="+string(v)+" meter");
+disp(" Minimum Specific Energy ="+string(E)+" meter");
diff --git a/2021/CH10/EX10.2/EX10_2.pdf b/2021/CH10/EX10.2/EX10_2.pdf
new file mode 100755
index 000000000..b72ca8960
--- /dev/null
+++ b/2021/CH10/EX10.2/EX10_2.pdf
diff --git a/2021/CH10/EX10.2/EX10_2.sce b/2021/CH10/EX10.2/EX10_2.sce
new file mode 100755
index 000000000..003c6f0f5
--- /dev/null
+++ b/2021/CH10/EX10.2/EX10_2.sce
@@ -0,0 +1,14 @@
+//Finding of bed slope and conveyance of channel
+//Given
+q=0.15;
+B=.70;
+y=.40;
+C=60;
+A=B*y;
+P=B+(2*y);
+R=(A/P);
+//To Find
+sb=((q^2)*(P))/((A^3)*C^2)
+K=A*C*sqrt(R);
+disp("Bed of slope ="+string(sb)+" no units");
+disp("conveyance of channel ="+string(K)+" m^3/sec");
diff --git a/2021/CH10/EX10.22/EX10_22.pdf b/2021/CH10/EX10.22/EX10_22.pdf
new file mode 100755
index 000000000..c56fabebf
--- /dev/null
+++ b/2021/CH10/EX10.22/EX10_22.pdf
diff --git a/2021/CH10/EX10.22/EX10_22.sce b/2021/CH10/EX10.22/EX10_22.sce
new file mode 100755
index 000000000..d559fd459
--- /dev/null
+++ b/2021/CH10/EX10.22/EX10_22.sce
@@ -0,0 +1,14 @@
+//Finding of Water surface Slope
+//Given
+sb=1/4000;
+sf=.00004;
+T=10;
+B=10;
+g=9.81;
+y=1.5;
+v=1;
+//To Find
+A=B*y;
+q=A*v;
+z=(sb-sf)/(1-((q^2*T)/(g*A^3)));
+disp("Water surface slope ="+string(z)+"no units ")
diff --git a/2021/CH10/EX10.23/EX10_23.pdf b/2021/CH10/EX10.23/EX10_23.pdf
new file mode 100755
index 000000000..ff781b8cb
--- /dev/null
+++ b/2021/CH10/EX10.23/EX10_23.pdf
diff --git a/2021/CH10/EX10.23/EX10_23.sce b/2021/CH10/EX10.23/EX10_23.sce
new file mode 100755
index 000000000..8317e1869
--- /dev/null
+++ b/2021/CH10/EX10.23/EX10_23.sce
@@ -0,0 +1,10 @@
+//Finding of discharge at section-1
+//Given
+T=30;
+dy=0.06;
+dt=3600;
+dx=1000;
+q2=35;
+//To Find
+q1=q2+((T*dy)/dt)*dx;
+disp("Discharge at section-1 ="+string(q1)+" m^3/sec");
diff --git a/2021/CH10/EX10.3/EX10_3.pdf b/2021/CH10/EX10.3/EX10_3.pdf
new file mode 100755
index 000000000..012b29b01
--- /dev/null
+++ b/2021/CH10/EX10.3/EX10_3.pdf
diff --git a/2021/CH10/EX10.3/EX10_3.sce b/2021/CH10/EX10.3/EX10_3.sce
new file mode 100755
index 000000000..b5c71c150
--- /dev/null
+++ b/2021/CH10/EX10.3/EX10_3.sce
@@ -0,0 +1,13 @@
+//Finding of discharge through trapezoidal channel
+//Given
+B=6;
+z=1/3;
+C=60;
+y=3;
+sb=1/5000;
+//To Find
+A=(B+z*y)*y;
+P=B+(2*y*sqrt(1+z^2));
+R=A/P;
+q=A*C*sqrt(R*sb);
+disp("Discharge through Trapezoidal channel ="+string(q)+" m^3/sec");
diff --git a/2021/CH10/EX10.4/EX10_4.pdf b/2021/CH10/EX10.4/EX10_4.pdf
new file mode 100755
index 000000000..ab15c8666
--- /dev/null
+++ b/2021/CH10/EX10.4/EX10_4.pdf
diff --git a/2021/CH10/EX10.4/EX10_4.sce b/2021/CH10/EX10.4/EX10_4.sce
new file mode 100755
index 000000000..09506f4df
--- /dev/null
+++ b/2021/CH10/EX10.4/EX10_4.sce
@@ -0,0 +1,14 @@
+//Finding of Bottom slope, Conveyance
+//Given
+q=0.1;
+B=0.6;
+y=0.3;
+A=B*y;
+n=0.013;
+P=1.2;
+R=A/P;
+//To Find
+b=((q^2)*(P))/((A^3)*B^2)
+K=A*B*sqrt(R);
+disp("Bed of slope ="+string(b)+" no units");
+disp("conveyance of channel ="+string(K)+" m^3/sec");
diff --git a/2021/CH10/EX10.5/EX10_5.pdf b/2021/CH10/EX10.5/EX10_5.pdf
new file mode 100755
index 000000000..e417963c8
--- /dev/null
+++ b/2021/CH10/EX10.5/EX10_5.pdf
diff --git a/2021/CH10/EX10.5/EX10_5.sce b/2021/CH10/EX10.5/EX10_5.sce
new file mode 100755
index 000000000..422d0468f
--- /dev/null
+++ b/2021/CH10/EX10.5/EX10_5.sce
@@ -0,0 +1,13 @@
+//Finding of Bed slope of Trapezoidal channel
+//Given
+B=6;
+y=3;
+z=3/4;
+q=30;
+A=(B+(z*y))*y;
+P=B+(2*y)*sqrt(1+z^2);
+R=(A/P);
+n=0.0158;
+//To Find
+sb=((q^2)*n^2)/((A^2)*(R^(4/3)));
+disp("Bed slope of Trapezoidal channel ="+string(sb)+" no units");
diff --git a/2021/CH10/EX10.6/EX10_6.pdf b/2021/CH10/EX10.6/EX10_6.pdf
new file mode 100755
index 000000000..b63d970fd
--- /dev/null
+++ b/2021/CH10/EX10.6/EX10_6.pdf
diff --git a/2021/CH10/EX10.6/EX10_6.sce b/2021/CH10/EX10.6/EX10_6.sce
new file mode 100755
index 000000000..a2497e786
--- /dev/null
+++ b/2021/CH10/EX10.6/EX10_6.sce
@@ -0,0 +1,14 @@
+//Finding of discharge through triangular channel
+//Given
+y=4;
+theta=60;
+b=theta/2;
+n=0.0182;
+sb=1/1000;
+T=2*tan(b)*y;z=tan(b); disp(z);
+A=0.5*T*y;
+P=2*sqrt(y^2+(y*tan(b))^2);
+R=A/P;disp(A);disp(P);disp(R);
+//To Find
+q=A*(1/n)*(R)^2/3*(sb)^1/2;
+disp("discharge through triangular channel ="+string(q)+" m^3/sec");
diff --git a/2021/CH10/EX10.7/EX10_7.pdf b/2021/CH10/EX10.7/EX10_7.pdf
new file mode 100755
index 000000000..0af9147fa
--- /dev/null
+++ b/2021/CH10/EX10.7/EX10_7.pdf
diff --git a/2021/CH10/EX10.7/EX10_7.sce b/2021/CH10/EX10.7/EX10_7.sce
new file mode 100755
index 000000000..b3bedb50c
--- /dev/null
+++ b/2021/CH10/EX10.7/EX10_7.sce
@@ -0,0 +1,10 @@
+//Finding of Diameter of circular channel
+//Given
+q=1;
+n=0.02;
+sb=1/10000;
+//To Find
+q=((%pi/8)*(1/n)*(1/sb)^(1/2)*(1/4)^(2/3));
+D=(8/%pi)*n*(1/sb)^(-0.5)*(4)^(2/3);
+D1=(D)^(3/8);
+disp("Diameter of Circular Pipe ="+string(D)+" meter");
diff --git a/2021/CH10/EX10.8/EX10_8.pdf b/2021/CH10/EX10.8/EX10_8.pdf
new file mode 100755
index 000000000..eb9032a01
--- /dev/null
+++ b/2021/CH10/EX10.8/EX10_8.pdf
diff --git a/2021/CH10/EX10.8/EX10_8.sce b/2021/CH10/EX10.8/EX10_8.sce
new file mode 100755
index 000000000..7b28eab71
--- /dev/null
+++ b/2021/CH10/EX10.8/EX10_8.sce
@@ -0,0 +1,15 @@
+//Finding of Dimemsions
+//Given
+q=0.5;
+sb=1/3000;
+c=60;
+n=0.015;
+//To Find
+y=q/(2*c*(1/2)^(0.5)*(sb)^(1/2));
+y1=y^(2/5);
+b=2*y1;
+y2=q/(2*(1/n)*(1/2)^(2/3)*(sb)^(1/2));
+y3=(y2)^(3/8);
+b1=2*y3;
+disp("Economical Dimensions ="+string(b)+" meter");
+disp("Economical Dimensions ="+string(b1)+" meter");
diff --git a/2021/CH11/EX11.1/EX11_1.pdf b/2021/CH11/EX11.1/EX11_1.pdf
new file mode 100755
index 000000000..5f2918813
--- /dev/null
+++ b/2021/CH11/EX11.1/EX11_1.pdf
diff --git a/2021/CH11/EX11.1/EX11_1.sce b/2021/CH11/EX11.1/EX11_1.sce
new file mode 100755
index 000000000..43962c63a
--- /dev/null
+++ b/2021/CH11/EX11.1/EX11_1.sce
@@ -0,0 +1,14 @@
+//Finding of Pressure Difference
+//Given
+mu=0.09;
+spgr=0.8;
+rho=800;
+D=0.08;
+L=15;
+//To Find
+A=(%pi/4)*D^2;
+q=(50/10)*(1/rho);
+v=q/A;
+p=(128*mu*q*L)/(%pi*D^4);
+p1=p/10000;
+disp(" Pressure Difference ="+string(p1)+" N/cm^2"); 
diff --git a/2021/CH11/EX11.2/EX11_2.pdf b/2021/CH11/EX11.2/EX11_2.pdf
new file mode 100755
index 000000000..d8fc8069a
--- /dev/null
+++ b/2021/CH11/EX11.2/EX11_2.pdf
diff --git a/2021/CH11/EX11.2/EX11_2.sce b/2021/CH11/EX11.2/EX11_2.sce
new file mode 100755
index 000000000..7cea72a88
--- /dev/null
+++ b/2021/CH11/EX11.2/EX11_2.sce
@@ -0,0 +1,15 @@
+//Finding of Pressure Drop
+//Given
+mu=0.15;
+spgr=.9;
+rho=900;
+D=.055;
+L=325;
+R=D/2;
+q=.0037;
+//To Find
+P=(128*mu*q*L)/(%pi*D^4);
+p1=P/100;
+x=(p1/L)*R;
+x1=x*10^4;
+disp("Pressure Drop ="+string(x1)+" N/m^2")
diff --git a/2021/CH11/EX11.3/EX11_3.sce b/2021/CH11/EX11.3/EX11_3.sce
new file mode 100755
index 000000000..2b16868cd
--- /dev/null
+++ b/2021/CH11/EX11.3/EX11_3.sce
@@ -0,0 +1,15 @@
+//Finding of Pressure Gradienr,Avg velocity,Reynolds number
+//Given
+mu=.5;
+spgr=1.2;
+rho=1200;
+D=.1;
+x=147.15;
+//To Find
+dp=-(x*4)/D;
+dp1=-dp;
+v=(1/(32*mu))*(-dp)*D^2;
+R=(rho*v*D)/mu;
+disp("Pressure Gradient ="+string(dp1)+" N/m^3");
+disp("Average Velocity ="+string(v)+" N/m^3");
+disp("Reynolds Number ="+string(R)+" N/m^3");
diff --git a/2021/CH11/EX11.3/Ex11_3.pdf b/2021/CH11/EX11.3/Ex11_3.pdf
new file mode 100755
index 000000000..0cb7a4f96
--- /dev/null
+++ b/2021/CH11/EX11.3/Ex11_3.pdf
diff --git a/2021/CH11/EX11.4/EX11_4.pdf b/2021/CH11/EX11.4/EX11_4.pdf
new file mode 100755
index 000000000..3623165a3
--- /dev/null
+++ b/2021/CH11/EX11.4/EX11_4.pdf
diff --git a/2021/CH11/EX11.4/EX11_4.sce b/2021/CH11/EX11.4/EX11_4.sce
new file mode 100755
index 000000000..c632fa987
--- /dev/null
+++ b/2021/CH11/EX11.4/EX11_4.sce
@@ -0,0 +1,12 @@
+//Finding of Power Required
+//Given
+L=100;
+D=0.1;
+q=0.01;
+mu=0.8;
+//To Find
+A=(%pi/4)*D^2;
+v=q/A;
+p=(32*q*mu*v*L)/D^2;
+P=p/100;
+disp("Power Required ="+string(P)+" KiloWatts");
diff --git a/2021/CH11/EX11.5/EX11_5.sce b/2021/CH11/EX11.5/EX11_5.sce
new file mode 100755
index 000000000..d405bb153
--- /dev/null
+++ b/2021/CH11/EX11.5/EX11_5.sce
@@ -0,0 +1,16 @@
+//Finding of Pressuure gradient.Avg velocity,Discharge,Shear at wall
+//Given
+mu=0.02;
+B=0.01;
+b=1;
+v=2;
+//To Find
+A=B*b;
+dp=-((16*mu)/B^2);
+dp1=-dp;
+V=(B^2/(12*mu))*(-dp);
+q=A*V;
+x=(-dp*(B/2))
+disp(" Pressuure Gradient ="+string(dp1)+" N/m^2 per meter");
+disp(" Avg velocity ="+string(V)+" m/sec");
+disp(" Shear at wall ="+string(x)+" N//m^2");
diff --git a/2021/CH11/EX11.5/Ex11_5.pdf b/2021/CH11/EX11.5/Ex11_5.pdf
new file mode 100755
index 000000000..4dd3ec556
--- /dev/null
+++ b/2021/CH11/EX11.5/Ex11_5.pdf
diff --git a/2021/CH11/EX11.6/EX11_6.sce b/2021/CH11/EX11.6/EX11_6.sce
new file mode 100755
index 000000000..31fd7385a
--- /dev/null
+++ b/2021/CH11/EX11.6/EX11_6.sce
@@ -0,0 +1,14 @@
+//Finding of Pressure Gradient , Shear at wall
+//Given
+D=15;
+f=0.05;
+r=4;
+tau=0.01962;
+//To Find
+R=64/f;
+dp=-(tau*(2/r));
+dp1=-dp;
+r1=D/2;
+tau2=(tau*r1)/r;
+disp("Pressure Gradient ="+string(dp1)+" N/m^3");
+disp(" Shear at wall ="+string(tau2)+" N/cm^2");
diff --git a/2021/CH11/EX11.6/Ex11_6.pdf b/2021/CH11/EX11.6/Ex11_6.pdf
new file mode 100755
index 000000000..cad60bd82
--- /dev/null
+++ b/2021/CH11/EX11.6/Ex11_6.pdf
diff --git a/2021/CH12/EX12.1/EX12_1.pdf b/2021/CH12/EX12.1/EX12_1.pdf
new file mode 100755
index 000000000..9c5aae2a3
--- /dev/null
+++ b/2021/CH12/EX12.1/EX12_1.pdf
diff --git a/2021/CH12/EX12.1/EX12_1.sce b/2021/CH12/EX12.1/EX12_1.sce
new file mode 100755
index 000000000..00a62a9a9
--- /dev/null
+++ b/2021/CH12/EX12.1/EX12_1.sce
@@ -0,0 +1,13 @@
+//Finding the type of boundary
+//Given
+ks=0.20*10^-3;
+tau=7.848;
+nu=0.01*10^-4;
+rho=1000;
+//To Find
+v=sqrt(tau/rho);
+R=(v*ks)/nu;
+disp("R ="+string(R)+" no units");
+if(R>4);
+if(R<60);
+disp("Flow is Transitional");    
diff --git a/2021/CH12/EX12.4/EX12_4.pdf b/2021/CH12/EX12.4/EX12_4.pdf
new file mode 100755
index 000000000..70a65ea04
--- /dev/null
+++ b/2021/CH12/EX12.4/EX12_4.pdf
diff --git a/2021/CH12/EX12.4/EX12_4.sce b/2021/CH12/EX12.4/EX12_4.sce
new file mode 100755
index 000000000..8f58bd217
--- /dev/null
+++ b/2021/CH12/EX12.4/EX12_4.sce
@@ -0,0 +1,19 @@
+//Finding of  Power Lost 
+//Given
+D=0.6;
+L=1000;
+Q=0.6;
+ks=0.003;
+rho=1000;
+g=9.81;
+c=50;
+//To Find
+//For Turbulant Flow
+A=(%pi/4)*D^2;
+a=2*log(3.71*(D/ks));
+b=sqrt(a);
+v=Q/A;
+f=1/c;
+Hf=(f*L*v^2)/(2*g*D);
+P=(rho*g*Q*Hf)/1000;
+disp(" Power Lost ="+string(P)+" Kilowatt");
diff --git a/2021/CH12/EX12.5/EX12_5.pdf b/2021/CH12/EX12.5/EX12_5.pdf
new file mode 100755
index 000000000..ca902bd91
--- /dev/null
+++ b/2021/CH12/EX12.5/EX12_5.pdf
diff --git a/2021/CH12/EX12.5/EX12_5.sce b/2021/CH12/EX12.5/EX12_5.sce
new file mode 100755
index 000000000..d6fe9289d
--- /dev/null
+++ b/2021/CH12/EX12.5/EX12_5.sce
@@ -0,0 +1,24 @@
+//Finding of Friction Factor
+//Given
+D=0.1;
+ks=0.0025;
+v=2;
+v1=10^-6;
+//To Find
+//case-1
+R=(v*D)/v1;
+fa=(1.785*log10(R))-1.424;
+a=(fa)^2;
+f1=1/a;
+//case-2
+fb=2*log10((3.71*D)/ks);
+b=(fb)^2;
+f2=1/b;
+//Case-3
+fc=-(2*log10((ks/3.71*D)+(5.186/R^(0.89))));
+c=(fc)^2;
+f3=1/c;
+disp(" Friction Factor for");
+disp("Smooth Turbulent flow ="+string(f1)+" no units ");
+disp("Rough Turbulent flow ="+string(f2)+" no units ");
+disp("Smooth and Rough Turbulent flow ="+string(f3)+" no units ");
diff --git a/2021/CH13/EX13.2/EX13_2.pdf b/2021/CH13/EX13.2/EX13_2.pdf
new file mode 100755
index 000000000..ed3c26bac
--- /dev/null
+++ b/2021/CH13/EX13.2/EX13_2.pdf
diff --git a/2021/CH13/EX13.2/EX13_2.sce b/2021/CH13/EX13.2/EX13_2.sce
new file mode 100755
index 000000000..2523477a8
--- /dev/null
+++ b/2021/CH13/EX13.2/EX13_2.sce
@@ -0,0 +1,21 @@
+//Finding of Boundary layer thickness , Drag Force
+//Given
+x=1;
+L=1.5;
+b=1.2;
+vs=0.25;
+mu=0.001;
+rho=1000;
+x2=1.2;
+L2=1.2;
+//To Find
+A=L*b;
+R=(rho*vs*x)/mu;
+t=(5.477*x)/sqrt(R);
+tau=(0.365*mu*vs*sqrt(R))/x;
+R1=(rho*vs*L)/mu;
+Cd=1.46/sqrt(R1);
+Fd=(1/2)*Cd*rho*(vs)^2*A;
+disp("Boundary Layer Thickness ="+string(t)+" meter");
+disp("Drag Force ="+string(Fd)+" Newtons");
+
diff --git a/2021/CH13/EX13.3/EX13_3.pdf b/2021/CH13/EX13.3/EX13_3.pdf
new file mode 100755
index 000000000..be61e5d5e
--- /dev/null
+++ b/2021/CH13/EX13.3/EX13_3.pdf
diff --git a/2021/CH13/EX13.3/EX13_3.sce b/2021/CH13/EX13.3/EX13_3.sce
new file mode 100755
index 000000000..8866b2181
--- /dev/null
+++ b/2021/CH13/EX13.3/EX13_3.sce
@@ -0,0 +1,21 @@
+//Finding of //Finding of Boundary layer thickness , Drag Force
+//Given
+x=1.5;
+L=2;
+b=1.4;
+vs=0.2;
+mu=0.001;
+rho=1000;
+//To Find
+A=L*b;
+R=(rho*vs*x)/mu;
+t=(4.64*x)/sqrt(R);
+t1=t*1000;
+tau=(0.323*mu*vs*sqrt(R))/x;
+R1=(rho*vs*L)/mu;
+Cd=1.292/sqrt(R);
+Fd=(1/2)*Cd*rho*(vs)^2*(2*A);
+disp("Co-efficient of Drag ="+string(Cd)+" no units")
+disp("Boundary Layer Thickness ="+string(t1)+" millimeter");
+disp("Drag Force ="+string(Fd)+" Newtons");
+
diff --git a/2021/CH14/EX14.10/EX14_10.pdf b/2021/CH14/EX14.10/EX14_10.pdf
new file mode 100755
index 000000000..b6324dbd2
--- /dev/null
+++ b/2021/CH14/EX14.10/EX14_10.pdf
diff --git a/2021/CH14/EX14.10/EX14_10.sce b/2021/CH14/EX14.10/EX14_10.sce
new file mode 100755
index 000000000..64baba6b6
--- /dev/null
+++ b/2021/CH14/EX14.10/EX14_10.sce
@@ -0,0 +1,9 @@
+//Finding of "n" of the Model
+//Given
+Lm=1;
+Lp=64;
+Np=0.02;
+//To Find
+A=sqrt(Lp/Lm);
+Nm=A*(Lm/Lp)^(2/3)*Np;
+disp("Manning n of the model ="+string(Nm)+" No units");
diff --git a/2021/CH14/EX14.11/EX14_11.pdf b/2021/CH14/EX14.11/EX14_11.pdf
new file mode 100755
index 000000000..852745124
--- /dev/null
+++ b/2021/CH14/EX14.11/EX14_11.pdf
diff --git a/2021/CH14/EX14.11/EX14_11.sce b/2021/CH14/EX14.11/EX14_11.sce
new file mode 100755
index 000000000..3ad8a717a
--- /dev/null
+++ b/2021/CH14/EX14.11/EX14_11.sce
@@ -0,0 +1,12 @@
+//Finding of Qm.Nm
+//Given
+Qp=3000;
+Np=0.025;
+L1=1000;
+L2=100;
+//To Find
+B=sqrt(L2);
+Qm=Qp/(L1*L2*B);
+Nm=(Qp/Qm)*1/(((L1*L2*(L2)^(2/3))/Np)*B*sqrt(1/L1));
+disp(" Qm ="+string(Qm)+" m^3/sec");
+disp(" Nm ="+string(Nm)+" No units");
diff --git a/2021/CH14/EX14.12/EX14_12.pdf b/2021/CH14/EX14.12/EX14_12.pdf
new file mode 100755
index 000000000..72f276baf
--- /dev/null
+++ b/2021/CH14/EX14.12/EX14_12.pdf
diff --git a/2021/CH14/EX14.12/EX14_12.sce b/2021/CH14/EX14.12/EX14_12.sce
new file mode 100755
index 000000000..994e49caa
--- /dev/null
+++ b/2021/CH14/EX14.12/EX14_12.sce
@@ -0,0 +1,10 @@
+//Finding of (Vm/Vp) and (Np/Nm)
+//Given
+L1=1/5000;
+L2=1/256;
+Qr=1/(2*10^7);
+//To Find
+Vr=(1/L1)*(1/L2)*Qr;
+Nr=Vr*((L2)^-(2/3)*(L2)^-(1/2)*(1/L1)^-(1/2));
+disp("Vm/Vp ="+string(Vr)+" m/sec");
+disp("Np/Nm ="+string(Nr)+" No units");
diff --git a/2021/CH14/EX14.13/EX14_13.pdf b/2021/CH14/EX14.13/EX14_13.pdf
new file mode 100755
index 000000000..ce1c754e2
--- /dev/null
+++ b/2021/CH14/EX14.13/EX14_13.pdf
diff --git a/2021/CH14/EX14.13/EX14_13.sce b/2021/CH14/EX14.13/EX14_13.sce
new file mode 100755
index 000000000..75db6e6ba
--- /dev/null
+++ b/2021/CH14/EX14.13/EX14_13.sce
@@ -0,0 +1,15 @@
+//Finding of Qm,Lm,Hm
+//Given
+Lp=16;
+Lm=1;
+Hp=4;
+L1=150;
+H1=7.2;
+H2=16;
+//To Find
+Hm=H1*(Lm/Lp);
+lm=L1*(Lm/Lp);
+Qm=(Lp/Lm)^2*(Hp/H2)^(1/2);
+disp("Lm ="+string(lm)+" meter");
+disp("Hm ="+string(Hm)+" meter");
+disp("Qm ="+string(Qm)+" m^3/sec");
diff --git a/2021/CH14/EX14.6/EX14_6.pdf b/2021/CH14/EX14.6/EX14_6.pdf
new file mode 100755
index 000000000..fb2d5e155
--- /dev/null
+++ b/2021/CH14/EX14.6/EX14_6.pdf
diff --git a/2021/CH14/EX14.6/EX14_6.sce b/2021/CH14/EX14.6/EX14_6.sce
new file mode 100755
index 000000000..3a32b0c3d
--- /dev/null
+++ b/2021/CH14/EX14.6/EX14_6.sce
@@ -0,0 +1,11 @@
+//Finding of velocity , discharge of prototype
+//Given
+qm=2;
+vm=1.5;
+lp=36;
+lm=1;
+//To Find
+vp=sqrt(lp/lm)*vm;
+qp=(lp/lm)^2*(vp/vm)*qm;
+disp("Velocity of Prototype ="+string(vp)+" m/sec");
+disp("Dischage of Prototype ="+string(qp)+" m^3/sec");
diff --git a/2021/CH14/EX14.7/EX14_7.pdf b/2021/CH14/EX14.7/EX14_7.pdf
new file mode 100755
index 000000000..ea45a1864
--- /dev/null
+++ b/2021/CH14/EX14.7/EX14_7.pdf
diff --git a/2021/CH14/EX14.7/EX14_7.sce b/2021/CH14/EX14.7/EX14_7.sce
new file mode 100755
index 000000000..56caf62f1
--- /dev/null
+++ b/2021/CH14/EX14.7/EX14_7.sce
@@ -0,0 +1,15 @@
+//Finding of Velocity of Prototype
+//Given
+vm=30;
+lm=100;
+lp=1;
+Am=0.018*10^-4;
+Ap=0.012*10^-4;
+rho1=1030;
+rho2=1.24;
+Fm=60;
+//To Find
+vp=(Ap/Am)*(lp/lm)*vm;
+Fp=Fm*(lm/lp)^2*(vp/vm)^2*(rho1/rho2);
+disp("Velocity of Prototype ="+string(vp)+" m/sec");
+disp("Resistance of Prototype ="+string(Fp)+" Newton");
diff --git a/2021/CH14/EX14.8/EX14_8.pdf b/2021/CH14/EX14.8/EX14_8.pdf
new file mode 100755
index 000000000..a043551bb
--- /dev/null
+++ b/2021/CH14/EX14.8/EX14_8.pdf
diff --git a/2021/CH14/EX14.8/EX14_8.sce b/2021/CH14/EX14.8/EX14_8.sce
new file mode 100755
index 000000000..582e96543
--- /dev/null
+++ b/2021/CH14/EX14.8/EX14_8.sce
@@ -0,0 +1,14 @@
+//Finding of Velocity of Model
+//Given
+vp=20;
+lm=1;
+lp=15;
+rho1=1024;
+rho2=1000;
+Fp=600;
+Fm=0.12;
+//To Find
+vm=sqrt(lm/lp)*vp;
+Fp=Fm*(lm/lp)^2*(vp/vm)^2*(rho1/rho2);
+disp("Velocity of Prototype ="+string(vm)+" m/sec");
+disp("Resistance of Prototype ="+string(Fp)+" Newton");
diff --git a/2021/CH14/EX14.9/EX14_9.pdf b/2021/CH14/EX14.9/EX14_9.pdf
new file mode 100755
index 000000000..c43c6cd25
--- /dev/null
+++ b/2021/CH14/EX14.9/EX14_9.pdf
diff --git a/2021/CH14/EX14.9/EX14_9.sce b/2021/CH14/EX14.9/EX14_9.sce
new file mode 100755
index 000000000..fb8ead2ab
--- /dev/null
+++ b/2021/CH14/EX14.9/EX14_9.sce
@@ -0,0 +1,11 @@
+//Finding of discharge through Model
+//Given
+A=50;
+B=10;
+C=sqrt(10);
+Qp=1.5;
+//To Find
+D=A*B;disp(D);
+Qm=(D)*(1/C);
+Qm1=Qp/Qm;
+disp(" Discharge Through Model ="+string(Qm1)+" m^3/sec");
diff --git a/2021/CH15/EX15.1/EX15_1.pdf b/2021/CH15/EX15.1/EX15_1.pdf
new file mode 100755
index 000000000..33b273b6b
--- /dev/null
+++ b/2021/CH15/EX15.1/EX15_1.pdf
diff --git a/2021/CH15/EX15.1/EX15_1.sce b/2021/CH15/EX15.1/EX15_1.sce
new file mode 100755
index 000000000..e09698926
--- /dev/null
+++ b/2021/CH15/EX15.1/EX15_1.sce
@@ -0,0 +1,17 @@
+//Finding of Velocities
+//Given
+T1=293;
+T2=293;
+P1=40;
+P2=35;
+R=287;
+A1=30*10^-4;
+A2=15*10^-4;
+Q=0.15;
+//To Find
+rho1=P1/(R*T1);
+V1=Q/(A1*rho1*10000);
+rho2=P2/(R*T2);
+V2=Q/(A2*rho2*10000);
+disp("Velocity at Section-1 ="+string(V1)+" m/sec");
+disp("Velocity at Section-2 ="+string(V2)+" m/sec");
diff --git a/2021/CH15/EX15.2/EX15_2.pdf b/2021/CH15/EX15.2/EX15_2.pdf
new file mode 100755
index 000000000..0276e4777
--- /dev/null
+++ b/2021/CH15/EX15.2/EX15_2.pdf
diff --git a/2021/CH15/EX15.2/EX15_2.sce b/2021/CH15/EX15.2/EX15_2.sce
new file mode 100755
index 000000000..7d1cc6ae6
--- /dev/null
+++ b/2021/CH15/EX15.2/EX15_2.sce
@@ -0,0 +1,9 @@
+//Finding of Speed of Sound waves
+//Given
+k=1.4;
+R=287;
+T=293;
+//To Find
+C=sqrt(k*R*T);
+C1=C*(18/5);
+disp("Speed of Sound waves ="+string(C1)+" Km/hr");
diff --git a/2021/CH15/EX15.3/EX15_3.pdf b/2021/CH15/EX15.3/EX15_3.pdf
new file mode 100755
index 000000000..ebe43f34a
--- /dev/null
+++ b/2021/CH15/EX15.3/EX15_3.pdf
diff --git a/2021/CH15/EX15.3/EX15_3.sce b/2021/CH15/EX15.3/EX15_3.sce
new file mode 100755
index 000000000..50065ac48
--- /dev/null
+++ b/2021/CH15/EX15.3/EX15_3.sce
@@ -0,0 +1,12 @@
+//Finding of Mach Number
+//Given
+k=1.4;
+R=287;
+T=288;
+V=900;
+//To Find
+C=sqrt(k*R*T);
+C1=C*(18/5);
+disp("Speed of Sound waves ="+string(C1)+" Km/hr");
+M=V/C1;
+disp("Mach Number = "+string(M)+" No units");
diff --git a/2021/CH15/EX15.4/EX15_4.pdf b/2021/CH15/EX15.4/EX15_4.pdf
new file mode 100755
index 000000000..a37227162
--- /dev/null
+++ b/2021/CH15/EX15.4/EX15_4.pdf
diff --git a/2021/CH15/EX15.4/EX15_4.sce b/2021/CH15/EX15.4/EX15_4.sce
new file mode 100755
index 000000000..dc30e0942
--- /dev/null
+++ b/2021/CH15/EX15.4/EX15_4.sce
@@ -0,0 +1,11 @@
+//Finding of Speed
+//Given
+k=1.4;
+R=287;
+T=233;
+M=1.8;
+//To Find
+C=sqrt(k*R*T);
+C1=C*(18/5);
+V=C1*M;
+disp("Speed of Aeroplane ="+string(V)+" Km/hr");
diff --git a/2021/CH15/EX15.5/EX15_5.pdf b/2021/CH15/EX15.5/EX15_5.pdf
new file mode 100755
index 000000000..af44a9e6c
--- /dev/null
+++ b/2021/CH15/EX15.5/EX15_5.pdf
diff --git a/2021/CH15/EX15.5/EX15_5.sce b/2021/CH15/EX15.5/EX15_5.sce
new file mode 100755
index 000000000..b3f8c73e5
--- /dev/null
+++ b/2021/CH15/EX15.5/EX15_5.sce
@@ -0,0 +1,12 @@
+//Finding of Velocity of Projectile
+//Given
+theta=30;
+k=1.4;
+R=287;
+T=268;
+//To Find
+Ma=sin(theta);
+C=sqrt(k*R*T);
+V=Ma*C;
+disp("Velocity of Projectile ="+string(V)+" m/sec");
+
diff --git a/2021/CH15/EX15.6/EX15_6.pdf b/2021/CH15/EX15.6/EX15_6.pdf
new file mode 100755
index 000000000..f440f26fd
--- /dev/null
+++ b/2021/CH15/EX15.6/EX15_6.pdf
diff --git a/2021/CH15/EX15.6/EX15_6.sce b/2021/CH15/EX15.6/EX15_6.sce
new file mode 100755
index 000000000..ada5ea2fa
--- /dev/null
+++ b/2021/CH15/EX15.6/EX15_6.sce
@@ -0,0 +1,12 @@
+//Finding of Mach Number and Mach Angle
+//Given
+k=1.4;
+R=287;
+T=263;
+V=1200;
+//To Find
+C=sqrt(k*R*T);
+Ma=V/C;
+alpha=asind(1/Ma);
+disp("Mach Number ="+string(Ma)+" No units");
+disp("Mach Angle ="+string(alpha)+" Degrees");
diff --git a/2021/CH15/EX15.7/EX15_7.pdf b/2021/CH15/EX15.7/EX15_7.pdf
new file mode 100755
index 000000000..e95d9b262
--- /dev/null
+++ b/2021/CH15/EX15.7/EX15_7.pdf
diff --git a/2021/CH15/EX15.7/EX15_7.sce b/2021/CH15/EX15.7/EX15_7.sce
new file mode 100755
index 000000000..ed8142058
--- /dev/null
+++ b/2021/CH15/EX15.7/EX15_7.sce
@@ -0,0 +1,21 @@
+//Finding of Mach's Number
+//Given
+k=1.4;
+R=287;
+T=273;
+T1=273-15;
+v=900;
+p1=8*10^4;
+//To Find
+V=v*(5/18);
+C=sqrt(k*R*T);
+Ma=V/C;
+ps=p1*((1+((k-1)/2)*Ma^2)^(k/(k-1)));
+Ps=ps*10^-4;
+Ts=T1*((1+((k-1)/2)*Ma^2));
+rho=ps/(R*T);
+t=Ts-T;
+disp("Mach Number ="+string(Ma)+" No Units");
+disp("Density ="+string(rho)+" Kg/m^3");
+disp("Pressure ="+string(Ps)+" N/cm^2");
+disp("Temperature ="+string(t)+" celcius")
diff --git a/2021/CH15/EX15.8/EX15_8.pdf b/2021/CH15/EX15.8/EX15_8.pdf
new file mode 100755
index 000000000..baec1bf3f
--- /dev/null
+++ b/2021/CH15/EX15.8/EX15_8.pdf
diff --git a/2021/CH15/EX15.8/EX15_8.sce b/2021/CH15/EX15.8/EX15_8.sce
new file mode 100755
index 000000000..b849a0af0
--- /dev/null
+++ b/2021/CH15/EX15.8/EX15_8.sce
@@ -0,0 +1,12 @@
+//Finding of velocity at the outlet of a nozzle
+//Given
+k=1.4;
+P1=294.3;
+P2=137.34;
+T1=303;
+R=287;
+//To Find
+rho=P1/(R*T1);
+V2=sqrt((2*k/(k-1))*(P1/rho)*(1-(P2/P1)^((k-1)/k)));
+disp("velocity at the outlet of a nozzle ="+string(V2)+" m/sec");
+
diff --git a/2021/CH15/EX15.9/EX15_9.pdf b/2021/CH15/EX15.9/EX15_9.pdf
new file mode 100755
index 000000000..b02a3a371
--- /dev/null
+++ b/2021/CH15/EX15.9/EX15_9.pdf
diff --git a/2021/CH15/EX15.9/EX15_9.sce b/2021/CH15/EX15.9/EX15_9.sce
new file mode 100755
index 000000000..d1ece5c0d
--- /dev/null
+++ b/2021/CH15/EX15.9/EX15_9.sce
@@ -0,0 +1,16 @@
+//Finding of Mass Flow Rate
+//Given
+D1=0.4;
+D2=0.2;
+P1=27.468*10^4;
+P2=25.506*10^4;
+T1=293;
+k=1.4;
+R=287;
+//To Find
+A1=(%pi/4)*D1^2;
+A2=(%pi/4)*D2^2;
+rho1=P1/(R*T1);
+rho2=((rho1^(1.4)*P2)/P1)^(1/1.4);
+m=rho2*A2*sqrt((2*k/(k-1))*(P1/rho1)*(1-(P2/P1)^((k-1/k)))/(1-(P2/P1)^(2/k))*(A2/A1)^2);
+disp("Mass Flow Rate ="+string(m)+" Kg/sec");
diff --git a/2021/CH16/EX16.1/EX16_1.pdf b/2021/CH16/EX16.1/EX16_1.pdf
new file mode 100755
index 000000000..58435291d
--- /dev/null
+++ b/2021/CH16/EX16.1/EX16_1.pdf
diff --git a/2021/CH16/EX16.1/EX16_1.sce b/2021/CH16/EX16.1/EX16_1.sce
new file mode 100755
index 000000000..422713697
--- /dev/null
+++ b/2021/CH16/EX16.1/EX16_1.sce
@@ -0,0 +1,20 @@
+//Finding of Lift,Drag,Power Required
+//Given
+A=4;
+V=40*(5/18);
+Cd1=0.8;
+Cd2=0.2;
+rho=1.25;
+//To Find
+FL=Cd1*A*rho*((V^2)/2);
+Fd=Cd2*A*rho*((V^2)/2);
+F=sqrt(FL^2+Fd^2);
+P=Fd*V;
+P1=P/1000;
+theta=(FL/Fd);
+theta1=(tan(theta))^-1;
+disp("Lift Force ="+string(FL)+" Newton");
+disp("Power Required ="+string(P1)+" Kilo Watts");
+disp("Drag Force ="+string(Fd)+" Newton");
+disp("Resultant Force ="+string(F)+" Newton");
+disp("Angle of Flow Direction ="+string(theta1)+" degrees");
diff --git a/2021/CH16/EX16.3/EX16_3.pdf b/2021/CH16/EX16.3/EX16_3.pdf
new file mode 100755
index 000000000..feeb26e7a
--- /dev/null
+++ b/2021/CH16/EX16.3/EX16_3.pdf
diff --git a/2021/CH16/EX16.3/EX16_3.sce b/2021/CH16/EX16.3/EX16_3.sce
new file mode 100755
index 000000000..b875b980a
--- /dev/null
+++ b/2021/CH16/EX16.3/EX16_3.sce
@@ -0,0 +1,11 @@
+//Finding of Diameter
+//Given
+W=80*9.81;
+Fd=80*9.81;
+V=25;
+Cd=0.5;
+rho=1.25;
+//To Find
+D=(2*Fd)/(Cd*rho*(V^2)*(4/%pi));
+D1=sqrt(D);
+disp("Diameter ="+string(D1)+" meter");
diff --git a/2021/CH16/EX16.4/EX16_4.pdf b/2021/CH16/EX16.4/EX16_4.pdf
new file mode 100755
index 000000000..bc8813abd
--- /dev/null
+++ b/2021/CH16/EX16.4/EX16_4.pdf
diff --git a/2021/CH16/EX16.4/EX16_4.sce b/2021/CH16/EX16.4/EX16_4.sce
new file mode 100755
index 000000000..4353e04db
--- /dev/null
+++ b/2021/CH16/EX16.4/EX16_4.sce
@@ -0,0 +1,14 @@
+//Finding of Coefficient of Lift ,Drag
+//Given
+A=25;
+P=588.6*(7/10);
+FL=19620;
+V=200*(5/18);
+rho=1000;
+FD=7416;
+//To Find
+FD=(P*1000)/(V);
+Cd=(FD*2)/(rho*A*(V^2));
+Cl=(FL*2)/(rho*A*(V^2));
+disp("Coefficient Of Lift ="+string(Cl)+" No Units");
+disp("Coefficient Of Drag ="+string(Cd)+" No Units");
diff --git a/2021/CH16/EX16.5/EX16_5.pdf b/2021/CH16/EX16.5/EX16_5.pdf
new file mode 100755
index 000000000..f768780d2
--- /dev/null
+++ b/2021/CH16/EX16.5/EX16_5.pdf
diff --git a/2021/CH16/EX16.5/EX16_5.sce b/2021/CH16/EX16.5/EX16_5.sce
new file mode 100755
index 000000000..40101c8c4
--- /dev/null
+++ b/2021/CH16/EX16.5/EX16_5.sce
@@ -0,0 +1,11 @@
+//Finding of Weight
+//Given
+D=0.05;
+v=1.5*10^-4;
+V=10;
+rho=1.25;
+Cd=0.5;
+//TO Find
+A=(%pi/4)*D^2;
+Fd=Cd*rho*A*((V^2)/2);
+disp("Weight of the ball ="+string(Fd)+" Newtons");
diff --git a/2021/CH16/EX16.6/EX16_6.pdf b/2021/CH16/EX16.6/EX16_6.pdf
new file mode 100755
index 000000000..d8a33d39b
--- /dev/null
+++ b/2021/CH16/EX16.6/EX16_6.pdf
diff --git a/2021/CH16/EX16.6/EX16_6.sce b/2021/CH16/EX16.6/EX16_6.sce
new file mode 100755
index 000000000..74a692514
--- /dev/null
+++ b/2021/CH16/EX16.6/EX16_6.sce
@@ -0,0 +1,42 @@
+//Finding of Circulation,Theoretical Drag,Actual Drag,Lift,Resultant,Direction
+//Given
+V=20;
+D=2;
+A=2*10;
+R=D/2;
+N=300;
+L=10;
+Cd=0.65;
+Cl=3.4;
+rho=1000;
+//To find
+Vp=(%pi*D*N)/60;
+//case 1
+C=2*(%pi)*R*Vp;
+disp("Circulation ="+string(C)+" m^2/sce");
+//case 2
+Fl=rho*V*L*C;
+disp("Theoretical Lift ="+string(Fl)+" Newtons");
+//case 3
+si=C/(4*(%pi)*V*R);
+theta1=(180+si);
+theta2=(360-si);
+disp("theta ="+string(theta1)+" Degrees");
+disp("theta ="+string(theta2)+" Degrees");
+//case 3
+FL=0.5*rho*A*V^2*Cl;
+disp("Lift Force ="+string(FL)+" Newtons");
+//case 4
+FD=0.5*rho*A*(V^2)*Cd;
+disp("Drag Force ="+string(FD)+" Newtons");
+//case 5
+F=sqrt((FL^2)+(FD^2));
+disp("Resultant Force ="+string(F)+" Newtons");
+//case 6
+theta=1/tan(FL/FD);
+disp("Direction ="+string(theta)+" Degrees");
+//case 7
+C1=4*(%pi)*V*R;
+Vp=C1/(2*(%pi)*R);
+N=(Vp*60)/(2*(%pi));
+disp("Speed ="+string(N)+" rpm");
diff --git a/2021/CH17/EX17.1/EX17_1.pdf b/2021/CH17/EX17.1/EX17_1.pdf
new file mode 100755
index 000000000..97e831d5f
--- /dev/null
+++ b/2021/CH17/EX17.1/EX17_1.pdf
diff --git a/2021/CH17/EX17.1/EX17_1.sce b/2021/CH17/EX17.1/EX17_1.sce
new file mode 100755
index 000000000..7ddd305ae
--- /dev/null
+++ b/2021/CH17/EX17.1/EX17_1.sce
@@ -0,0 +1,9 @@
+//Finding of Force exerted
+//Given
+rho=1000;
+d=0.04;
+V=25;
+//To Find
+A=(%pi/4)*d^2;
+P=rho*A*V^2;
+disp("Force Exerted ="+string(P)+" Newtons");
diff --git a/2021/CH17/EX17.11/EX17_11.pdf b/2021/CH17/EX17.11/EX17_11.pdf
new file mode 100755
index 000000000..3c665f161
--- /dev/null
+++ b/2021/CH17/EX17.11/EX17_11.pdf
diff --git a/2021/CH17/EX17.11/EX17_11.sce b/2021/CH17/EX17.11/EX17_11.sce
new file mode 100755
index 000000000..a57589ca3
--- /dev/null
+++ b/2021/CH17/EX17.11/EX17_11.sce
@@ -0,0 +1,18 @@
+//Finding of Propelling Force,Work Done ,Efficiency
+//Given
+Cv=0.97;
+g=9.81;
+H=6;
+rho=1000;
+u=4;
+d=0.15;
+//To Find
+V=Cv*sqrt(2*g*H);
+A=(%pi/4)*d^2;
+P=rho*A*(V+u)*V;
+W=P*u;
+E=(2*u*V)/(u+V)^2;
+E1=E*100;
+disp("Propelling Force ="+string(P)+" Newtons");
+disp("Work Done ="+string(W)+" N-m");
+disp("Efficiency ="+string(E1)+" Percentage");
diff --git a/2021/CH17/EX17.12/EX17_12.pdf b/2021/CH17/EX17.12/EX17_12.pdf
new file mode 100755
index 000000000..2a1b56b36
--- /dev/null
+++ b/2021/CH17/EX17.12/EX17_12.pdf
diff --git a/2021/CH17/EX17.12/EX17_12.sce b/2021/CH17/EX17.12/EX17_12.sce
new file mode 100755
index 000000000..ebc04a892
--- /dev/null
+++ b/2021/CH17/EX17.12/EX17_12.sce
@@ -0,0 +1,12 @@
+//Finding of Propelling Force, Efficiency
+//Given
+u=35*(5/18);
+V=25;
+a=0.04;
+rho=1000;
+//To Find
+P=rho*a*(V+u)*V;
+E=(2*u)/(V+(2*u));
+E1=E*100;
+disp("Propelling Force ="+string(P)+" Newtons");
+disp("Efficiency ="+string(E1)+" No Units");
diff --git a/2021/CH17/EX17.2/EX17_2.pdf b/2021/CH17/EX17.2/EX17_2.pdf
new file mode 100755
index 000000000..1e46d330e
--- /dev/null
+++ b/2021/CH17/EX17.2/EX17_2.pdf
diff --git a/2021/CH17/EX17.2/EX17_2.sce b/2021/CH17/EX17.2/EX17_2.sce
new file mode 100755
index 000000000..793761e7a
--- /dev/null
+++ b/2021/CH17/EX17.2/EX17_2.sce
@@ -0,0 +1,11 @@
+//Finding of Discharge
+//Given
+rho=1000;
+d=0.05;
+P=1226.25;
+//To Find
+A=(%pi/4)*d^2;
+V=P/(rho*A);
+V1=sqrt(V);
+Q=A*V1;
+disp("Discharge ="+string(Q)+" m^3/sec");
diff --git a/2021/CH17/EX17.3/EX17_3.pdf b/2021/CH17/EX17.3/EX17_3.pdf
new file mode 100755
index 000000000..67a33a42b
--- /dev/null
+++ b/2021/CH17/EX17.3/EX17_3.pdf
diff --git a/2021/CH17/EX17.3/EX17_3.sce b/2021/CH17/EX17.3/EX17_3.sce
new file mode 100755
index 000000000..92e5dd991
--- /dev/null
+++ b/2021/CH17/EX17.3/EX17_3.sce
@@ -0,0 +1,9 @@
+//Finding of Force Exerted
+//Given
+rho=1000;
+d=0.15;
+V=25;
+//To Find
+A=(%pi/4)*d^2;
+P=rho*A*V^2*sin(%pi/6);
+disp("Force Exerted ="+string(P)+" Newtons");
diff --git a/2021/CH17/EX17.4/EX17_4.pdf b/2021/CH17/EX17.4/EX17_4.pdf
new file mode 100755
index 000000000..cbf178a88
--- /dev/null
+++ b/2021/CH17/EX17.4/EX17_4.pdf
diff --git a/2021/CH17/EX17.4/EX17_4.sce b/2021/CH17/EX17.4/EX17_4.sce
new file mode 100755
index 000000000..dbda005f1
--- /dev/null
+++ b/2021/CH17/EX17.4/EX17_4.sce
@@ -0,0 +1,10 @@
+//Finding of Force Exerted
+//Given
+rho=1000;
+d=0.04;
+V=35;
+theta=180-125;
+//To FInd
+A=(%pi/4)*d^2;
+Fx=2*rho*A*V^2;
+disp("Force Exerted ="+string(Fx)+" Newtons");
diff --git a/2021/CH17/EX17.5/EX17_5.pdf b/2021/CH17/EX17.5/EX17_5.pdf
new file mode 100755
index 000000000..a59473d6b
--- /dev/null
+++ b/2021/CH17/EX17.5/EX17_5.pdf
diff --git a/2021/CH17/EX17.5/EX17_5.sce b/2021/CH17/EX17.5/EX17_5.sce
new file mode 100755
index 000000000..cccda93c6
--- /dev/null
+++ b/2021/CH17/EX17.5/EX17_5.sce
@@ -0,0 +1,13 @@
+//Finding of Force Exerted
+//Given
+rho=1000;
+d=0.07;
+V=25;
+theta=20;
+theta2=15;
+//To Find
+A=(%pi/4)*d^2
+Fx=rho*A*V^2*(sin(%pi/9)+cos(%pi/12));
+Fy=rho*A*V^2*(sin(%pi/9)-sin(%pi/12));
+disp("Fx ="+string(Fx)+" Newtons");
+disp("Fy ="+string(Fy)+" Newtons");
diff --git a/2021/CH17/EX17.6/EX17_6.pdf b/2021/CH17/EX17.6/EX17_6.pdf
new file mode 100755
index 000000000..01af74b1e
--- /dev/null
+++ b/2021/CH17/EX17.6/EX17_6.pdf
diff --git a/2021/CH17/EX17.6/EX17_6.sce b/2021/CH17/EX17.6/EX17_6.sce
new file mode 100755
index 000000000..faf70cbb5
--- /dev/null
+++ b/2021/CH17/EX17.6/EX17_6.sce
@@ -0,0 +1,11 @@
+//Finding of inclination
+rho=1000;
+d=0.03;
+V=16;
+w=125;
+//To Find
+A=(%pi/4)*d^2;
+P=rho*A*V^2;
+Q=P*(16/32);
+theta=asin((rho*A*V^2)/w);
+disp("Inclination ="+string(theta)+" degrees");
diff --git a/2021/CH17/EX17.7/EX17_7.sce b/2021/CH17/EX17.7/EX17_7.sce
new file mode 100755
index 000000000..4328ed513
--- /dev/null
+++ b/2021/CH17/EX17.7/EX17_7.sce
@@ -0,0 +1,13 @@
+//Finding of Vane Angle
+//Given
+V=40;
+u=20;
+alpha=30;
+b=90;
+u1=20;
+//TO Find
+theta=atand((V*sin(%pi/6))/((V*cos(%pi/6))-u));
+Vr=((V*sin(%pi/6))/(sin(theta)));
+pi=acosd(u1/Vr);
+disp("Vane angle at Inlet ="+string(theta)+" Degrees");
+disp("Vane angle at Outlet ="+string(pi)+" Degrees");
diff --git a/2021/CH17/EX17.7/Ex17_7.pdf b/2021/CH17/EX17.7/Ex17_7.pdf
new file mode 100755
index 000000000..d6b3750f2
--- /dev/null
+++ b/2021/CH17/EX17.7/Ex17_7.pdf
diff --git a/2021/CH18/EX18.1/EX18_1.pdf b/2021/CH18/EX18.1/EX18_1.pdf
new file mode 100755
index 000000000..de332417e
--- /dev/null
+++ b/2021/CH18/EX18.1/EX18_1.pdf
diff --git a/2021/CH18/EX18.1/EX18_1.sce b/2021/CH18/EX18.1/EX18_1.sce
new file mode 100755
index 000000000..88b465538
--- /dev/null
+++ b/2021/CH18/EX18.1/EX18_1.sce
@@ -0,0 +1,19 @@
+//Finding of Power delivered , Efficiency
+//Given
+u=35;
+Q=1;
+theta=10; 
+H=270;
+Cv=0.98;
+g=9.81;
+rho=1000;
+//To Find
+V=Cv*sqrt(2*g*H);
+Vr=V-u;
+Vw1=Vr*cos(%pi/18)-u;
+W=rho*(Q*(V+Vw1)*u);
+P=W/1000;
+E=(2*(V+Vw1)*u)/V^2;
+E1=E*100;
+disp("Power delivered ="+string(P)+" Kilo watts");
+disp("Hydraulic Efficiency ="+string(E1)+" percentage");
diff --git a/2021/CH18/EX18.10/EX18_10.sce b/2021/CH18/EX18.10/EX18_10.sce
new file mode 100755
index 000000000..9cf783063
--- /dev/null
+++ b/2021/CH18/EX18.10/EX18_10.sce
@@ -0,0 +1,20 @@
+//Finding of Vane Angle ,Head ,Velocity ,Efficiency
+//Given
+u=12;
+D=0.8;
+D1=1;
+Vw1=0;
+Hout=1;
+Vw=12;
+Vf=3;
+g=9.81;
+//To Find
+a=atand(Vf/Vw);
+V=sqrt(Vw^2+Vf^2);
+u1=(D1/D)*u;
+V1=u1*tan(%pi/9);
+H=((V1^2/(2*g))+1)+((Vw*u)/g);
+E=((Vw*u)/(g*H))*100;
+disp("Absolute Velocity ="+string(V)+" m/sec");
+disp("Vane Angle ="+string(a)+" degrees");
+disp("Efficiency ="+string(E)+" Percentage");
diff --git a/2021/CH18/EX18.10/Ex18_10.pdf b/2021/CH18/EX18.10/Ex18_10.pdf
new file mode 100755
index 000000000..544fbecd5
--- /dev/null
+++ b/2021/CH18/EX18.10/Ex18_10.pdf
diff --git a/2021/CH18/EX18.11/EX18_11.sce b/2021/CH18/EX18.11/EX18_11.sce
new file mode 100755
index 000000000..c00582d0e
--- /dev/null
+++ b/2021/CH18/EX18.11/EX18_11.sce
@@ -0,0 +1,24 @@
+//Finding of Angle ,Diameter ,Width
+//Given
+E=0.75;
+P=147.15;
+H=8;
+N=200;
+Vw1=0;
+Cv=0.3;
+g=9.81;
+rho=1000;
+//To Find
+u=Cv*sqrt(2*g*H);
+Vf=0.96*sqrt(2*g*H);
+E1=0.8;
+Vw=(E1*g*H)/u;
+a=atand(Vf/Vw);
+b=atand(Vf/(Vw-u));
+D=(60*u)/(%pi*N);
+Q=(E1*P*1000)/(E*rho*Vw*u);
+B=Q/(%pi*D*Vf);
+disp("Diameter ="+string(D)+" meter");
+disp(" Inlet angle ="+string(a)+" degrees");
+disp(" Outlet angle ="+string(b)+" degrees");
+disp("Width ="+string(B)+" meter");
diff --git a/2021/CH18/EX18.11/Ex18_11.pdf b/2021/CH18/EX18.11/Ex18_11.pdf
new file mode 100755
index 000000000..6ce0d39c2
--- /dev/null
+++ b/2021/CH18/EX18.11/Ex18_11.pdf
diff --git a/2021/CH18/EX18.14/EX18_14.pdf b/2021/CH18/EX18.14/EX18_14.pdf
new file mode 100755
index 000000000..752b2e60a
--- /dev/null
+++ b/2021/CH18/EX18.14/EX18_14.pdf
diff --git a/2021/CH18/EX18.14/EX18_14.sce b/2021/CH18/EX18.14/EX18_14.sce
new file mode 100755
index 000000000..0e257a19f
--- /dev/null
+++ b/2021/CH18/EX18.14/EX18_14.sce
@@ -0,0 +1,19 @@
+//Finding of Pressure Head ,Efficiency
+//Given
+Di=0.8;
+Do=1.2;
+V2=3;
+L=8;
+y=2;
+Hs=6;
+g=9.81;
+//To Find
+Q=(%pi/4)*Do^2*V2;
+V1=Q/((%pi/4)*Di^2);
+a=(V1^2/(2*g))-(V2^2/(2*g));
+b=0.25*(V2^2/(2*g));
+P=10.3-Hs-a-b;
+E=(a-b)/(V1^2/(2*g));
+E1=E*100;disp(V1);
+disp("Pressure Head ="+string(P)+" meter of water");
+disp("Efficiency ="+string(E1)+" Percentage");
diff --git a/2021/CH18/EX18.15/EX18_15.pdf b/2021/CH18/EX18.15/EX18_15.pdf
new file mode 100755
index 000000000..c24747afa
--- /dev/null
+++ b/2021/CH18/EX18.15/EX18_15.pdf
diff --git a/2021/CH18/EX18.15/EX18_15.sce b/2021/CH18/EX18.15/EX18_15.sce
new file mode 100755
index 000000000..a7502d7e6
--- /dev/null
+++ b/2021/CH18/EX18.15/EX18_15.sce
@@ -0,0 +1,11 @@
+//Finding of Speed ,Power developed
+//Given
+P1=8000;
+N1=90;
+H1=25;
+H2=15;
+//To Find
+N2=N1*(sqrt(H2)/sqrt(H1));
+P2=(P1*(H2)^(3/2))/(H1)^(3/2);
+disp("Speed ="+string(N2)+" rpm");
+disp("Power Developed ="+string(P2)+" Kilo watts");
diff --git a/2021/CH18/EX18.16/EX18_16.pdf b/2021/CH18/EX18.16/EX18_16.pdf
new file mode 100755
index 000000000..57a987ed0
--- /dev/null
+++ b/2021/CH18/EX18.16/EX18_16.pdf
diff --git a/2021/CH18/EX18.16/EX18_16.sce b/2021/CH18/EX18.16/EX18_16.sce
new file mode 100755
index 000000000..05fc0ba5e
--- /dev/null
+++ b/2021/CH18/EX18.16/EX18_16.sce
@@ -0,0 +1,13 @@
+//Finding of Specific speed,Power generated
+//Given
+H=30;
+N=300;
+Q=10;
+E=0.9;
+g=9.81;
+rho=1000;
+//To Find
+P=(E*rho*g*Q*H)/1000;
+Ns=(N*sqrt(P))/(H)^(5/4);
+disp("Power Developed ="+string(P)+" Kilowatts");
+disp("Specific Speed ="+string(Ns)+" rpm");
diff --git a/2021/CH18/EX18.2/EX18_2.pdf b/2021/CH18/EX18.2/EX18_2.pdf
new file mode 100755
index 000000000..8f69e695c
--- /dev/null
+++ b/2021/CH18/EX18.2/EX18_2.pdf
diff --git a/2021/CH18/EX18.2/EX18_2.sce b/2021/CH18/EX18.2/EX18_2.sce
new file mode 100755
index 000000000..ac0f0e6b8
--- /dev/null
+++ b/2021/CH18/EX18.2/EX18_2.sce
@@ -0,0 +1,22 @@
+//Finding of D,d number of jets
+//Given
+E=0.86;
+Dr=10;
+Cv=0.98;
+a=0.45;
+Sp=735.75*1000;
+H=200;
+g=9.81;
+N=800;
+rho=1000;
+//To Find
+V=Cv*sqrt(2*g*H);
+u=V*a;
+D=(60*u)/(%pi*N);
+d=(D/10);
+Q1=(%pi/4)*(d^2)*V;
+Q2=1/((E*rho*g*H)/Sp);
+j=Q2/Q1;disp(Q2);
+disp("D= "+string(D)+" meter");
+disp("d= "+string(d)+" meter");
+disp("Number of Jets ="+string(j)+" nos");
diff --git a/2021/CH18/EX18.3/EX18_3.pdf b/2021/CH18/EX18.3/EX18_3.pdf
new file mode 100755
index 000000000..05e896a20
--- /dev/null
+++ b/2021/CH18/EX18.3/EX18_3.pdf
diff --git a/2021/CH18/EX18.3/EX18_3.sce b/2021/CH18/EX18.3/EX18_3.sce
new file mode 100755
index 000000000..261ce0d57
--- /dev/null
+++ b/2021/CH18/EX18.3/EX18_3.sce
@@ -0,0 +1,18 @@
+//Finding of Power ,Efficiency
+//Given
+D=0.8;
+N=1000;
+a=15;
+Q=0.15;
+Cv=0.98;
+rho=1000'
+g=9.81;
+H=400;
+//To Find
+u=(%pi*D*N)/60;
+V=Cv*sqrt(2*g*H);
+P=(rho*g*Q*H)/1000;
+E=2*(V-u)*(1+cos(%pi/12))*u
+E1=(E/V^2)*100;
+disp("Power available ="+string(P)+" Kilo watts");
+disp("Hydraulic efficiency ="+string(E1)+" percentage");
diff --git a/2021/CH18/EX18.4/EX18_4.pdf b/2021/CH18/EX18.4/EX18_4.pdf
new file mode 100755
index 000000000..0c8e3b9d2
--- /dev/null
+++ b/2021/CH18/EX18.4/EX18_4.pdf
diff --git a/2021/CH18/EX18.4/EX18_4.sce b/2021/CH18/EX18.4/EX18_4.sce
new file mode 100755
index 000000000..cb4481be9
--- /dev/null
+++ b/2021/CH18/EX18.4/EX18_4.sce
@@ -0,0 +1,22 @@
+//Finding of Power delivered , Efficiency
+//Given
+Q=1.8;
+theta=12; 
+Hg=450;
+H=300;
+hf=Hg/3;
+Cv=0.98;
+g=9.81;
+a=0.46;
+rho=1000;
+//To Find
+V=Cv*sqrt(2*g*H);
+u=a*V;
+Vr=V-u;disp(V);
+Vw1=Vr*cos(%pi/15)-u;
+W=rho*(Q*(V+Vw1)*u);
+P=W/1000;disp(V);
+E=(2*(V+Vw1)*u)/V^2;
+E1=E*100;
+disp("Power delivered ="+string(P)+" Kilo watts");
+disp("Hydraulic Efficiency ="+string(E1)+" percentage");
diff --git a/2021/CH18/EX18.5/EX18_5.pdf b/2021/CH18/EX18.5/EX18_5.pdf
new file mode 100755
index 000000000..addf6a335
--- /dev/null
+++ b/2021/CH18/EX18.5/EX18_5.pdf
diff --git a/2021/CH18/EX18.5/EX18_5.sce b/2021/CH18/EX18.5/EX18_5.sce
new file mode 100755
index 000000000..fffb8f5f2
--- /dev/null
+++ b/2021/CH18/EX18.5/EX18_5.sce
@@ -0,0 +1,19 @@
+//Finding of Power Developed ,Force Exerted
+//Given
+d=0.13;
+a=15;
+H=400;
+Cv=0.97;
+b=0.45;
+g=9.81;
+rho=1000;
+//To Find
+A=(%pi/4)*d^2;
+u=b*sqrt(2*g*H);
+V=0.97*sqrt(2*g*H);
+Vr1=0.8*(V-u);
+Vw1=u-(Vr1*cos(%pi/15));
+Fx=rho*A*V*(V-Vw1);
+P=(Fx*u)/1000;
+disp("Force Exerted ="+string(Fx)+" Newton");
+disp("Power developed ="+string(P)+" Kilo Watts");
diff --git a/2021/CH18/EX18.6/EX18_6.pdf b/2021/CH18/EX18.6/EX18_6.pdf
new file mode 100755
index 000000000..a03533821
--- /dev/null
+++ b/2021/CH18/EX18.6/EX18_6.pdf
diff --git a/2021/CH18/EX18.6/EX18_6.sce b/2021/CH18/EX18.6/EX18_6.sce
new file mode 100755
index 000000000..f05b19f78
--- /dev/null
+++ b/2021/CH18/EX18.6/EX18_6.sce
@@ -0,0 +1,11 @@
+//Finding of Discharge, Width
+//Given
+D=1.2;
+D1=0.6;
+Vf=1.8;
+B=.20;
+//To Find
+Q=(%pi*D*B*Vf);
+B1=((D*B)/D1)*100;
+disp("Discharge ="+string(Q)+" m^3/sec");
+disp("Width ="+string(B1)+" centimeter");
diff --git a/2021/CH18/EX18.7/EX18_7.pdf b/2021/CH18/EX18.7/EX18_7.pdf
new file mode 100755
index 000000000..41324d871
--- /dev/null
+++ b/2021/CH18/EX18.7/EX18_7.pdf
diff --git a/2021/CH18/EX18.7/EX18_7.sce b/2021/CH18/EX18.7/EX18_7.sce
new file mode 100755
index 000000000..f70f060ed
--- /dev/null
+++ b/2021/CH18/EX18.7/EX18_7.sce
@@ -0,0 +1,21 @@
+//Finding of Discharge,Power developed,Efficiency
+//Given
+N=500;
+H=100;
+D=1;
+A=35;
+a=15;
+b=60;
+Vw1=0;
+g=9.81;
+rho=1000;
+//To Find
+u=(%pi*D*N)/60;
+Vw=(tan(%pi/3)*u)/1.464;
+Vf=Vw*tan(%pi/12);
+Q=A*Vf;
+P=(rho*g*Vw*u)/1000;
+E=((Vw*u)/(g*H))*100;
+disp("Discharge ="+string(Q)+" m^3/sec");
+disp("Power Developed ="+string(P)+" Kilo Watts");
+disp("Efficiency ="+string(E)+" Percentage");
diff --git a/2021/CH18/EX18.8/EX18_8.sce b/2021/CH18/EX18.8/EX18_8.sce
new file mode 100755
index 000000000..d3205ca9e
--- /dev/null
+++ b/2021/CH18/EX18.8/EX18_8.sce
@@ -0,0 +1,21 @@
+//Finding of Power developed,Outlet Vane Angle ,Speed
+//Given
+H=100;
+D=.675;
+D1=0.5
+B=0.15;
+B1=.225;
+g=9.81;
+rho=1000;
+Vf=3;
+Vw=3;
+//To Find
+u=Vf/tan(%pi/15);
+N=1/((%pi*D)/(u*60));
+u1=u*(D1/D);disp(Vf);
+Vf1=(D*B*Vf)/(0.15*B1);
+z=atand(Vf1/u1);
+P=(rho*%pi*B*D1*Vf*Vw*u)/10000;
+disp("Speed ="+string(N)+" rpm");
+disp("Power Developed ="+string(P)+" Kilo Watts");
+disp("Outlet Vane Angle ="+string(z)+" degrees");
diff --git a/2021/CH18/EX18.8/Ex18_8.pdf b/2021/CH18/EX18.8/Ex18_8.pdf
new file mode 100755
index 000000000..970e7c7c7
--- /dev/null
+++ b/2021/CH18/EX18.8/Ex18_8.pdf
diff --git a/2021/CH18/EX18.9/EX18_9.pdf b/2021/CH18/EX18.9/EX18_9.pdf
new file mode 100755
index 000000000..ee1facb9d
--- /dev/null
+++ b/2021/CH18/EX18.9/EX18_9.pdf
diff --git a/2021/CH18/EX18.9/EX18_9.sce b/2021/CH18/EX18.9/EX18_9.sce
new file mode 100755
index 000000000..36966cdd1
--- /dev/null
+++ b/2021/CH18/EX18.9/EX18_9.sce
@@ -0,0 +1,22 @@
+//Finding of Work Done ,Vane Angles
+//Given
+D=0.5;
+D1=1;
+Vw1=0;
+Vf=3;
+Vf1=3;
+g=9.81;
+rho=1000;
+N=250;
+//To Find
+u=(%pi*D*N)/60;
+u1=(%pi*D1*N)/60;
+Vw=Vf/tan(%pi/12);
+a=atand(Vf/(Vw-u));
+b=atand(Vf1/u1);
+W=(Vw*u)/g;
+E=(W/10);
+disp("Inlet Vane Angle ="+string(a)+" degrees");
+disp("Outlet Vane Angle ="+string(b)+" degrees");
+disp("Work Done ="+string(W)+" N-m/N");
+disp("Efficiency ="+string(E)+" Percentage");
diff --git a/2021/CH19/EX19.1/EX19_1.pdf b/2021/CH19/EX19.1/EX19_1.pdf
new file mode 100755
index 000000000..5f4b04f42
--- /dev/null
+++ b/2021/CH19/EX19.1/EX19_1.pdf
diff --git a/2021/CH19/EX19.1/EX19_1.sce b/2021/CH19/EX19.1/EX19_1.sce
new file mode 100755
index 000000000..1e5282a61
--- /dev/null
+++ b/2021/CH19/EX19.1/EX19_1.sce
@@ -0,0 +1,17 @@
+//Finding of workdone
+//Given
+D1=.6;
+D=0.3;
+a=20;
+b=30;
+N=1000;
+g=9.81;
+Vw=0;
+rho=1000;
+//To Find
+u=(%pi*D*N)/60;
+u1=(%pi*D1*N)/60;
+Vf=u*tan(%pi/9);
+Vw1=(u1*tan(%pi/6)-Vf)/tan(%pi/6);
+W=(Vw1*u1)/g;
+disp("Work Done ="+string(W)+" N-m/N");
diff --git a/2021/CH19/EX19.10/EX19_10.pdf b/2021/CH19/EX19.10/EX19_10.pdf
new file mode 100755
index 000000000..50c1653a8
--- /dev/null
+++ b/2021/CH19/EX19.10/EX19_10.pdf
diff --git a/2021/CH19/EX19.10/EX19_10.sce b/2021/CH19/EX19.10/EX19_10.sce
new file mode 100755
index 000000000..9f35c5264
--- /dev/null
+++ b/2021/CH19/EX19.10/EX19_10.sce
@@ -0,0 +1,23 @@
+//Finding of Head , Shaft Power
+//Given
+n=3;
+D1=0.4;
+B1=0.025;
+a=30;
+A=0.15;
+A1=0.0267;
+Emano=0.85;
+E=0.75;
+Q=0.06;
+N=1200;
+g=9.81;
+rho=1000;
+//To Find
+V=Q/A1;disp(V);
+u=(%pi*D1*N)/60;
+Vw1=(u*tan(%pi/6)-V)/tan(%pi/6);
+Hm=(Emano*u*Vw1)/g;
+Hm1=3*Hm;
+SP=(rho*g*Q*Hm1)/(1000*E);
+disp("Head ="+string(Hm1)+" meter");
+disp("Shaft Power ="+string(SP)+" Kilo watts");
diff --git a/2021/CH19/EX19.11/EX19_11.pdf b/2021/CH19/EX19.11/EX19_11.pdf
new file mode 100755
index 000000000..9b865735c
--- /dev/null
+++ b/2021/CH19/EX19.11/EX19_11.pdf
diff --git a/2021/CH19/EX19.11/EX19_11.sce b/2021/CH19/EX19.11/EX19_11.sce
new file mode 100755
index 000000000..4f4f3ea39
--- /dev/null
+++ b/2021/CH19/EX19.11/EX19_11.sce
@@ -0,0 +1,11 @@
+//Finding of Number of pumps
+//Given
+H=156;
+N=1000;
+Ns=20;
+Q=0.15;
+//To Find
+Hm=(N*sqrt(Q))/Ns;
+Hm1=(Hm)^(4/3);
+pumps=(H/Hm1);
+disp("Number of Pumps ="+string(pumps)+" Nos");
diff --git a/2021/CH19/EX19.12/EX19_12.pdf b/2021/CH19/EX19.12/EX19_12.pdf
new file mode 100755
index 000000000..41c4241be
--- /dev/null
+++ b/2021/CH19/EX19.12/EX19_12.pdf
diff --git a/2021/CH19/EX19.12/EX19_12.sce b/2021/CH19/EX19.12/EX19_12.sce
new file mode 100755
index 000000000..9828feb2b
--- /dev/null
+++ b/2021/CH19/EX19.12/EX19_12.sce
@@ -0,0 +1,16 @@
+//Finding of Head Discharge , Ratio of Power
+//Given
+Q1=0.035;
+H1=25;
+D1=0.5;
+N1=1200;
+D2=0.3;
+N2=2000;
+//To Find
+H=(D2*N2*sqrt(H1))/(D1*N1);
+H2=H^2;
+Q=(Q1*D2^3*N2)/(D1^3*N1);
+Pr=(D1/D2)^5*(N1/N2)^3;
+disp("Head ="+string(H2)+" meter");
+disp("Discharge ="+string(Q)+" m^3/sec");
+disp("Power Ratio ="+string(Pr)+" No Units");
diff --git a/2021/CH19/EX19.2/EX19_2.pdf b/2021/CH19/EX19.2/EX19_2.pdf
new file mode 100755
index 000000000..9d37bd7c5
--- /dev/null
+++ b/2021/CH19/EX19.2/EX19_2.pdf
diff --git a/2021/CH19/EX19.2/EX19_2.sce b/2021/CH19/EX19.2/EX19_2.sce
new file mode 100755
index 000000000..0b746ef96
--- /dev/null
+++ b/2021/CH19/EX19.2/EX19_2.sce
@@ -0,0 +1,25 @@
+//Finding of vane angle , Work done , Efficiency
+//Given
+D1=0.6;
+D=0.3;
+a=30;
+b=0.05;
+N=1200;
+g=9.81;
+Hm=75;
+Vf=3;
+rho=1000;
+B1=1;
+//To Find
+u=(%pi*D*N)/60;
+u1=(%pi*D1*N)/60;
+Q=%pi*D1*B1*Vf;
+a=atand(Vf/u);disp(u1);
+Vw1=((u1*tan(%pi/6))-Vf)/tan(%pi/6);
+W=(rho*g*Q*u1*Vw1)/g;
+W1=W/1000;
+E=((g*Hm)/(u1*Vw1))*100;
+disp("Vane Angle ="+string(a)+" degrees");
+disp("Work Done ="+string(W1)+" KW/sec");
+disp("Manometric Efficiency ="+string(E)+" Percentage");
+
diff --git a/2021/CH19/EX19.3/EX19_3.pdf b/2021/CH19/EX19.3/EX19_3.pdf
new file mode 100755
index 000000000..c8b7a7df0
--- /dev/null
+++ b/2021/CH19/EX19.3/EX19_3.pdf
diff --git a/2021/CH19/EX19.3/EX19_3.sce b/2021/CH19/EX19.3/EX19_3.sce
new file mode 100755
index 000000000..436959f16
--- /dev/null
+++ b/2021/CH19/EX19.3/EX19_3.sce
@@ -0,0 +1,15 @@
+//Finding of Workdone
+//Given
+D1=0.3;
+D=0.15;
+a=30;
+b=25;
+N=1450;
+g=9.81;
+//To Find
+u=(%pi*D1*N)/60;
+u1=(%pi*D*N)/60;
+Vf=u*tan(%pi/6);
+Vw1=(-u1*tan(%pi/7)+Vf)/tan(%pi/7);
+W=(Vw1*u1)/g;
+disp("Work Done ="+string(W)+" Nm/N");
diff --git a/2021/CH19/EX19.4/EX19_4.pdf b/2021/CH19/EX19.4/EX19_4.pdf
new file mode 100755
index 000000000..f45715508
--- /dev/null
+++ b/2021/CH19/EX19.4/EX19_4.pdf
diff --git a/2021/CH19/EX19.4/EX19_4.sce b/2021/CH19/EX19.4/EX19_4.sce
new file mode 100755
index 000000000..e31458cf1
--- /dev/null
+++ b/2021/CH19/EX19.4/EX19_4.sce
@@ -0,0 +1,18 @@
+//Finding of Vane Angle
+//Given
+N=1450;
+Hm=23;
+D1=0.25;
+B1=0.05;
+Emano=0.75;
+g=9.81;
+Q=1.25;
+//To Find
+u=(%pi*D1*N)/60;
+Vw1=(Emano*u)/(g*Hm);
+z=u-Vw1;
+Vf1=z*tan(%pi/6);
+Vf1=Q/(%pi*D1*B1);
+a=Vf1/(u-Vw1);
+b=atand(a);
+disp("Vane Angle ="+string(b)+" degrees");
diff --git a/2021/CH19/EX19.5/EX19_5.pdf b/2021/CH19/EX19.5/EX19_5.pdf
new file mode 100755
index 000000000..8fa2c97b5
--- /dev/null
+++ b/2021/CH19/EX19.5/EX19_5.pdf
diff --git a/2021/CH19/EX19.5/EX19_5.sce b/2021/CH19/EX19.5/EX19_5.sce
new file mode 100755
index 000000000..b94541e72
--- /dev/null
+++ b/2021/CH19/EX19.5/EX19_5.sce
@@ -0,0 +1,16 @@
+//Finding of Discharge
+//Given
+N=1000;
+Hm=15;
+D1=0.3;
+B1=0.05;
+a=30;
+Emano=0.92;
+g=9.81;
+//To Find
+u=(%pi*D1*N)/60;
+Vw1=(Emano*u)/(g*Hm);
+z=u-Vw1;
+Vf1=z*tan(%pi/6);
+Q=%pi*D1*B1*Vf1;
+disp("Discharge ="+string(Q)+" m^3/sec");
diff --git a/2021/CH19/EX19.6/EX19_6.pdf b/2021/CH19/EX19.6/EX19_6.pdf
new file mode 100755
index 000000000..5ac87498b
--- /dev/null
+++ b/2021/CH19/EX19.6/EX19_6.pdf
diff --git a/2021/CH19/EX19.6/EX19_6.sce b/2021/CH19/EX19.6/EX19_6.sce
new file mode 100755
index 000000000..c67791124
--- /dev/null
+++ b/2021/CH19/EX19.6/EX19_6.sce
@@ -0,0 +1,18 @@
+//Finding of Power Required
+//Given
+Q=0.03;
+H=18.25;
+L=90;
+dp=0.1;
+E=0.75;
+f=0.04;
+g=9.81;
+rho=1000;
+//Given
+V=Q/((%pi/4)*dp^2);
+loss=(f*L*V^2)/(2*g*dp);
+a=V^2/(2*g);
+Hm=H+loss+a;
+SP=(rho*g*Q*Hm)/(E*1000);
+disp("Power required ="+string(SP)+" Kilowatts");
+
diff --git a/2021/CH19/EX19.7/EX19_7.pdf b/2021/CH19/EX19.7/EX19_7.pdf
new file mode 100755
index 000000000..dcd8a972f
--- /dev/null
+++ b/2021/CH19/EX19.7/EX19_7.pdf
diff --git a/2021/CH19/EX19.7/EX19_7.sce b/2021/CH19/EX19.7/EX19_7.sce
new file mode 100755
index 000000000..ae5926a8a
--- /dev/null
+++ b/2021/CH19/EX19.7/EX19_7.sce
@@ -0,0 +1,12 @@
+//Finding of Minimum Speed
+//Given
+Hm=7.5;
+D1=1;
+D=0.5;
+g=9.81;
+//To Find
+u=(4/3)*(Hm*2*g);
+u1=sqrt(u);
+N=(60*u1)/(4*%pi);
+disp(" u="+string(u1)+" m/sec");
+disp("Minimum Speed ="+string(N)+" rpm");
diff --git a/2021/CH19/EX19.8/EX19_8.pdf b/2021/CH19/EX19.8/EX19_8.pdf
new file mode 100755
index 000000000..ffdde0368
--- /dev/null
+++ b/2021/CH19/EX19.8/EX19_8.pdf
diff --git a/2021/CH19/EX19.8/EX19_8.sce b/2021/CH19/EX19.8/EX19_8.sce
new file mode 100755
index 000000000..aa053aea6
--- /dev/null
+++ b/2021/CH19/EX19.8/EX19_8.sce
@@ -0,0 +1,12 @@
+//Finding of Mininmum Speed
+//Given
+D=0.3;
+D1=0.6;
+Vf1=2.5;
+a=45;
+Emano=0.75;
+//To Find
+u=(%pi*D)/60;
+N=(-2.5)*(120*Emano**D1)/(%pi*(D1^2-D^2));
+N1=-N+(1/u);
+disp("Speed ="+string(N1)+" rpm");
diff --git a/2021/CH19/EX19.9/EX19_9.pdf b/2021/CH19/EX19.9/EX19_9.pdf
new file mode 100755
index 000000000..912fea755
--- /dev/null
+++ b/2021/CH19/EX19.9/EX19_9.pdf
diff --git a/2021/CH19/EX19.9/EX19_9.sce b/2021/CH19/EX19.9/EX19_9.sce
new file mode 100755
index 000000000..8ba1719c3
--- /dev/null
+++ b/2021/CH19/EX19.9/EX19_9.sce
@@ -0,0 +1,17 @@
+//Finding of Manometric head
+//Given
+D1=0.4;
+B1=0.025;
+Q=0.06;
+N=1000;
+a=30;
+g=9.81;
+Emano=0.8;
+//To Find
+u=(%pi*D1*N)/60;
+Vf=Q/(%pi*D1*B1);
+Vw1=(-Vf*tan(%pi/6)+u);
+H=(Vw1*u)/g;
+Hm=(Emano*u*Vw1)/g;
+Hm1=2*Hm;
+disp("Head Developed ="+string(Hm1)+" meter");
diff --git a/2021/CH2/EX2.1/EX2_1.sce b/2021/CH2/EX2.1/EX2_1.sce
new file mode 100755
index 000000000..b7af3c3d2
--- /dev/null
+++ b/2021/CH2/EX2.1/EX2_1.sce
@@ -0,0 +1,25 @@
+//Finding of Pressure Intensity
+//Given
+p=73.575*10^4;       //Pressure in N/mm^2
+Patm=76;
+rho=1000;             //Density in kg/m^3
+spgr=13.6;
+g=9.81;                //Gravitational force in m/s^2
+//To Find 
+//Gauge units
+P1=p/(rho*g);
+P2=p/(spgr*rho*g);
+//Absolute units
+P3=(Patm*spgr)/100;
+P4=(Patm*spgr*rho*g)/100;
+P5=P2/10000;
+P6=p+P5;
+P7=75+P5;
+P8=5.5147+(Patm/100);
+disp("Gauge Units ");
+disp("Pressure ="+string(P1)+" meter of water");
+disp("Pressure ="+string(P2)+" meter of mercury");
+disp("Pressure ="+string(P3)+" meter of water");
+disp("Pressure ="+string(P5)+" N/cm^2");
+disp("Pressure ="+string(P7)+" meter of water");
+disp("Pressure ="+string(P8)+" meter of Mercury");
diff --git a/2021/CH2/EX2.1/Ex2_1.pdf b/2021/CH2/EX2.1/Ex2_1.pdf
new file mode 100755
index 000000000..702076751
--- /dev/null
+++ b/2021/CH2/EX2.1/Ex2_1.pdf
diff --git a/2021/CH2/EX2.2/EX2_2.sce b/2021/CH2/EX2.2/EX2_2.sce
new file mode 100755
index 000000000..70806f197
--- /dev/null
+++ b/2021/CH2/EX2.2/EX2_2.sce
@@ -0,0 +1,10 @@
+//Finding of Depth of Water
+//Given
+p=100.5525*10^4;                //pressure intensity in N/m^2
+spgr=1.025;                //Specific gravity
+rho=1000;                  //Density of water in kg/m^3
+g=9.81;                    //Gravitational force due to acceleration in m/sec^2
+w=rho*g;
+//To Find
+h=p/w;
+disp("Depth of Water = "+string(h)+" m");
diff --git a/2021/CH2/EX2.2/Ex2_2.pdf b/2021/CH2/EX2.2/Ex2_2.pdf
new file mode 100755
index 000000000..b2bce50f7
--- /dev/null
+++ b/2021/CH2/EX2.2/Ex2_2.pdf
diff --git a/2021/CH2/EX2.3/EX2_3.pdf b/2021/CH2/EX2.3/EX2_3.pdf
new file mode 100755
index 000000000..51e5f0a0b
--- /dev/null
+++ b/2021/CH2/EX2.3/EX2_3.pdf
diff --git a/2021/CH2/EX2.3/EX2_3.sce b/2021/CH2/EX2.3/EX2_3.sce
new file mode 100755
index 000000000..02d14e386
--- /dev/null
+++ b/2021/CH2/EX2.3/EX2_3.sce
@@ -0,0 +1,13 @@
+//Finding of Height
+//Given
+p=4.9*10^4;           //Pressure intensity in N/mm^2
+rho=1000;             //Density of water in kg/m^3
+g=9.81;               //gravitational force in m/sec^2
+spgr=0.8;             //Specific gravity of oil
+w=rho*g;
+w1=rho*g*spgr;
+//To Find
+h=p/w;
+h1=p/w1;
+disp("Height of water ="+string(h)+" m");
+disp("Height of oil ="+string(h1)+" m");
diff --git a/2021/CH2/EX2.4/EX2_4.sce b/2021/CH2/EX2.4/EX2_4.sce
new file mode 100755
index 000000000..2f5ead1f3
--- /dev/null
+++ b/2021/CH2/EX2.4/EX2_4.sce
@@ -0,0 +1,14 @@
+//Finding of Pressure intensity
+//Given
+spgr=0.8;               //Specific gravity of oil
+spgr1=1.5;                
+rho=1000;               //Density of water in Kg/m^3
+g=9.81;                 //Gravitational force m/s^2
+h1=2;                    //Depth in meter
+//To Find
+w=rho*g*spgr;
+p=w*h1;                  //Pressure at the interface
+disp("Pressure Intensity ="+string(p)+" N/m^2");
+w1=rho*g*spgr1;
+p1=1.5969+w1;
+disp("Pressure Intensity at Bottom ="+string(p1)+" N/m^2");
diff --git a/2021/CH2/EX2.4/Ex2_4.pdf b/2021/CH2/EX2.4/Ex2_4.pdf
new file mode 100755
index 000000000..b7ec90a2d
--- /dev/null
+++ b/2021/CH2/EX2.4/Ex2_4.pdf
diff --git a/2021/CH2/EX2.5/EX2_5.sce b/2021/CH2/EX2.5/EX2_5.sce
new file mode 100755
index 000000000..f3531c5ef
--- /dev/null
+++ b/2021/CH2/EX2.5/EX2_5.sce
@@ -0,0 +1,15 @@
+//Finding of Pressure
+//Given
+spgr1=0.8;              //specific gravity of liquid
+spgr2=13.6;             //specific gravity of mercury 
+h1=0.6;                 //height in left limb in meter
+h2=0.15;                //height in right limb in meter
+g=9.81;                 //gravitaional force in m/s^2
+rho1=spgr1*1000;        //density of liquid in Kg/m^3
+rho2=spgr2*1000;        //density of mecury  in Kg/m^3
+p=10.13;
+//To Find
+p1=(rho2*g*h2)-(rho1*g*h1);
+disp("Pressure at A ="+string(p1)+" N/m^2");
+p2=p1/100+p;
+disp("Absolute Pressure "+string(p2)+" N/cm^2")
diff --git a/2021/CH2/EX2.5/Ex2_5.pdf b/2021/CH2/EX2.5/Ex2_5.pdf
new file mode 100755
index 000000000..e5ecb87dc
--- /dev/null
+++ b/2021/CH2/EX2.5/Ex2_5.pdf
diff --git a/2021/CH2/EX2.6/EX2_6.pdf b/2021/CH2/EX2.6/EX2_6.pdf
new file mode 100755
index 000000000..a4b2717ed
--- /dev/null
+++ b/2021/CH2/EX2.6/EX2_6.pdf
diff --git a/2021/CH2/EX2.6/EX2_6.sce b/2021/CH2/EX2.6/EX2_6.sce
new file mode 100755
index 000000000..c0f4c540c
--- /dev/null
+++ b/2021/CH2/EX2.6/EX2_6.sce
@@ -0,0 +1,17 @@
+//Finding of Pressure difference
+//Given
+P1=10.8*10^4;       //Pressure in N/mm^2
+P2=17.16*10^4;
+rho=1000;            //Density in kg/m^3
+g=9.81;            //Gravitational force in m/s^2
+spgr1=1.594;
+spgr2=13.6;
+spgr3=0.8;
+z1=4;                 //height in meter
+//To Find
+left=P1+((spgr1*rho*g)*z1)+(spgr2*rho*g);
+right=((spgr3*rho*g)*1.5)+P2;
+h=left/(12*right);
+h1=h*100;
+disp("Height ="+string(h)+" meter of mercury");
+disp("Height ="+string(h1)+" centimeter of mercury");
diff --git a/2021/CH2/EX2.7/EX2_7.sce b/2021/CH2/EX2.7/EX2_7.sce
new file mode 100755
index 000000000..2b43bd80a
--- /dev/null
+++ b/2021/CH2/EX2.7/EX2_7.sce
@@ -0,0 +1,16 @@
+//Finding of Pressure 
+//Given
+//Left Limb
+h=0.6;
+rho=1000;
+g=9.81;
+//Right Limb
+h1=0.45;
+spgr1=13.6;
+h2=0.30;
+spgr2=0.88;
+//To Find
+P=(h1*spgr1)+(h2*spgr2)-(h);
+disp("Pressure is ="+string(P)+" cm of water");
+p1=P*rho*g;
+disp("Pressure is = "+string(p1)+" N/m^2");
diff --git a/2021/CH2/EX2.7/Ex2_7.pdf b/2021/CH2/EX2.7/Ex2_7.pdf
new file mode 100755
index 000000000..f7b77164c
--- /dev/null
+++ b/2021/CH2/EX2.7/Ex2_7.pdf
diff --git a/2021/CH2/EX2.8/EX2_8.sce b/2021/CH2/EX2.8/EX2_8.sce
new file mode 100755
index 000000000..f03764dcc
--- /dev/null
+++ b/2021/CH2/EX2.8/EX2_8.sce
@@ -0,0 +1,15 @@
+//Finding of Elevation
+//Given
+//At Sea level
+p=760;    //pressure in mm of mercury
+rho=1000; //Density in kg/m^3
+spgr=13.6; 
+g=9.81; //gravitational force in m/sec^2
+p1=(p/1000)*rho*g*spgr;
+//At Mountain
+p2=735;
+p3=(p2/1000)*rho*g*spgr;
+rho1=1.2;
+//To Find
+h=(p1-p3)/(rho1*g);
+disp("Elevation is ="+string(h)+" meter");
diff --git a/2021/CH2/EX2.8/Ex2_8.pdf b/2021/CH2/EX2.8/Ex2_8.pdf
new file mode 100755
index 000000000..144adce46
--- /dev/null
+++ b/2021/CH2/EX2.8/Ex2_8.pdf
diff --git a/2021/CH2/EX2.9/EX2_9.pdf b/2021/CH2/EX2.9/EX2_9.pdf
new file mode 100755
index 000000000..04241c12a
--- /dev/null
+++ b/2021/CH2/EX2.9/EX2_9.pdf
diff --git a/2021/CH2/EX2.9/EX2_9.sce b/2021/CH2/EX2.9/EX2_9.sce
new file mode 100755
index 000000000..540306f30
--- /dev/null
+++ b/2021/CH2/EX2.9/EX2_9.sce
@@ -0,0 +1,12 @@
+//Finding of Pressure and Temperature
+//Given
+h=18.288;    //Height in  kilometer
+t0=288.15;
+l=6.5;       // Lengh in meter
+p0=101.18;   //Pressure in N/mm^2
+g=9.81;         //gravitational force in m/s^2
+//To Find
+T=t0-(l*h);
+disp("Temperature is ="+string(T)+" kelvin");
+Ps=p0*(1-(l*h/t0))^(g/(287.1*l));
+disp("Pressure is = "+string(Ps)+" KN/m^2");
diff --git a/2021/CH20/EX20.1/EX20_1.pdf b/2021/CH20/EX20.1/EX20_1.pdf
new file mode 100755
index 000000000..0bdd10674
--- /dev/null
+++ b/2021/CH20/EX20.1/EX20_1.pdf
diff --git a/2021/CH20/EX20.1/EX20_1.sce b/2021/CH20/EX20.1/EX20_1.sce
new file mode 100755
index 000000000..70af83e46
--- /dev/null
+++ b/2021/CH20/EX20.1/EX20_1.sce
@@ -0,0 +1,16 @@
+//Finding of theoretical discharge ,Coefficient of Discharge ,Slip
+//Given
+N=30;
+Qac=0.012;
+d=0.25;
+L=0.5;
+//To Find
+A=(%pi/4)*d^2;
+Qth=(A*L*N)/60;
+S=Qth-Qac;
+Cd=Qac/Qth;
+S1=((Qth-Qac)/Qth)*100;
+disp("Theoretical Discharge ="+string(Qth)+" m^3/sec");
+disp("Co efficient of Discharge ="+string(Cd)+" No Units");
+disp("Slip ="+string(S)+" m^3/sec");
+disp("Percentage Slip ="+string(S1)+" No Units");
diff --git a/2021/CH20/EX20.2/EX20_2.pdf b/2021/CH20/EX20.2/EX20_2.pdf
new file mode 100755
index 000000000..56a7f86b6
--- /dev/null
+++ b/2021/CH20/EX20.2/EX20_2.pdf
diff --git a/2021/CH20/EX20.2/EX20_2.sce b/2021/CH20/EX20.2/EX20_2.sce
new file mode 100755
index 000000000..76084a852
--- /dev/null
+++ b/2021/CH20/EX20.2/EX20_2.sce
@@ -0,0 +1,17 @@
+//Finding of Slip ,Power required
+//Given
+N=50;
+Qac=0.015;
+L=0.4;
+D=0.25;
+hd=25;
+hs=4;
+rho=1000;
+g=9.81;
+//To Find
+A=(%pi/4)*D^2;
+Qth=(2*A*L*N)/60;
+S=Qth-Qac;
+P=((2*rho*g*A*L*N)*(hs+hd))/60000;
+disp("Slip ="+string(S)+" m^3.sec");
+disp("Power required ="+string(P)+" Kilo Watts");
diff --git a/2021/CH20/EX20.3/EX20_3.pdf b/2021/CH20/EX20.3/EX20_3.pdf
new file mode 100755
index 000000000..700fc8024
--- /dev/null
+++ b/2021/CH20/EX20.3/EX20_3.pdf
diff --git a/2021/CH20/EX20.3/EX20_3.sce b/2021/CH20/EX20.3/EX20_3.sce
new file mode 100755
index 000000000..31840764a
--- /dev/null
+++ b/2021/CH20/EX20.3/EX20_3.sce
@@ -0,0 +1,17 @@
+//Finding of Pressure Head
+//Given
+D=0.15;
+L=0.3;
+hs=4;
+N=40;
+l=5;
+ds=0.1;
+p=10.3;
+g=9.81;
+//To Find
+A=(%pi/4)*D^2;
+a=(%pi/4)*(ds)^2;
+r=L/2;
+Z=(2*%pi*N)/60;
+ha=(l/g)*(A/a)*r*Z^2;
+disp("Pressure Head ="+string(ha)+" meter");
diff --git a/2021/CH20/EX20.4/EX20_4.pdf b/2021/CH20/EX20.4/EX20_4.pdf
new file mode 100755
index 000000000..4a80d281f
--- /dev/null
+++ b/2021/CH20/EX20.4/EX20_4.pdf
diff --git a/2021/CH20/EX20.4/EX20_4.sce b/2021/CH20/EX20.4/EX20_4.sce
new file mode 100755
index 000000000..375819144
--- /dev/null
+++ b/2021/CH20/EX20.4/EX20_4.sce
@@ -0,0 +1,24 @@
+//Finding of Qth ,Pth ,ha
+//Given
+D=0.15;
+L=0.3;
+N=50;
+H=25;
+ld=22;
+dd=0.1;
+Qac=0.0042;
+rho=1000;
+g=9.81;
+//To Find
+A=(%pi/4)*D^2;
+a=(%pi/4)*(dd)^2;
+Z=(2*%pi*N)/60;
+r=L/2;
+Qth=(A*L*N)/60;
+Pth=(rho*g*Qth*H)/1000;
+S1=((Qth-Qac)/Qth)*100;
+ha=(ld/g)*(A/a)*r*Z^2;
+disp("Qth ="+string(Qth)+" m^3/sec");
+disp("Pth ="+string(Pth)+" Kilo Watts");
+disp("ha ="+string(ha)+" meter");
+
diff --git a/2021/CH20/EX20.5/EX20_5.pdf b/2021/CH20/EX20.5/EX20_5.pdf
new file mode 100755
index 000000000..f64c29471
--- /dev/null
+++ b/2021/CH20/EX20.5/EX20_5.pdf
diff --git a/2021/CH20/EX20.5/EX20_5.sce b/2021/CH20/EX20.5/EX20_5.sce
new file mode 100755
index 000000000..099538353
--- /dev/null
+++ b/2021/CH20/EX20.5/EX20_5.sce
@@ -0,0 +1,21 @@
+//Finding of Pmax ,Pressure at begining and end of Stroke
+//Given
+D=0.2;
+L=0.4;
+l=6;
+ds=0.1;
+hs=3.5;
+H=10.3;
+N=35;
+g=9.81;
+//To Find
+A=(%pi/4)*D^2;
+a=(%pi/4)*(ds)^2;
+r=L/2;
+Z=(2*%pi*N)/60;
+Pmax=(l/g)*(A/a)*r*Z^2;
+P=hs+Pmax;
+P1=H-P;
+disp("Pmax ="+string(Pmax)+" meter");
+disp("Pressure at Begining ="+string(P)+" meter");
+disp("Pressure at End ="+string(P1)+" meter");
diff --git a/2021/CH20/EX20.6/EX20_6.pdf b/2021/CH20/EX20.6/EX20_6.pdf
new file mode 100755
index 000000000..066f92882
--- /dev/null
+++ b/2021/CH20/EX20.6/EX20_6.pdf
diff --git a/2021/CH20/EX20.6/EX20_6.sce b/2021/CH20/EX20.6/EX20_6.sce
new file mode 100755
index 000000000..2c9db2af9
--- /dev/null
+++ b/2021/CH20/EX20.6/EX20_6.sce
@@ -0,0 +1,19 @@
+//Finding of Maximum Speed
+//Given
+D=0.125;
+L=0.3;
+hs=4.5;
+ds=0.075;
+l=6.8;
+h=2.6;
+g=9.81;
+H=10.3;
+//TO Find
+A=(%pi/4)*D^2;
+a=(%pi/4)*(ds)^2;
+r=L/2;
+ha=H-h-hs;
+Z=(ha*g*a)/(l*A*r);
+Z1=sqrt(Z);
+N=(Z1*60)/(2*%pi);
+disp("Maximum Speed ="+string(N)+" rpm");
diff --git a/2021/CH21/EX21.10/EX21_10.pdf b/2021/CH21/EX21.10/EX21_10.pdf
new file mode 100755
index 000000000..53d9a9f32
--- /dev/null
+++ b/2021/CH21/EX21.10/EX21_10.pdf
diff --git a/2021/CH21/EX21.10/EX21_10.sce b/2021/CH21/EX21.10/EX21_10.sce
new file mode 100755
index 000000000..fc3acf61e
--- /dev/null
+++ b/2021/CH21/EX21.10/EX21_10.sce
@@ -0,0 +1,10 @@
+//Finding of Diameter
+//Given
+D=25;
+P1=25;
+P2=120;
+//To Find
+A=(%pi/4)*D^2;
+d=(A*P1)/P2;
+d1=sqrt(d);
+disp("Diameter ="+string(d1)+" centimeter");
diff --git a/2021/CH21/EX21.11/EX21_11.pdf b/2021/CH21/EX21.11/EX21_11.pdf
new file mode 100755
index 000000000..101b0aac3
--- /dev/null
+++ b/2021/CH21/EX21.11/EX21_11.pdf
diff --git a/2021/CH21/EX21.11/EX21_11.sce b/2021/CH21/EX21.11/EX21_11.sce
new file mode 100755
index 000000000..fb204145d
--- /dev/null
+++ b/2021/CH21/EX21.11/EX21_11.sce
@@ -0,0 +1,11 @@
+//Finding of load on Ram
+//Given
+D=0.2;
+L=6;
+p=588.6*10^4;
+//To Find
+A=(%pi/4)*D^2;
+W=p*A;
+capacity=W*L;
+disp("Load ="+string(W)+" Newtons");
+disp("Capacity of the accumulator ="+string(capacity)+" N-m");
diff --git a/2021/CH21/EX21.12/EX21_12.sce b/2021/CH21/EX21.12/EX21_12.sce
new file mode 100755
index 000000000..c829e9b10
--- /dev/null
+++ b/2021/CH21/EX21.12/EX21_12.sce
@@ -0,0 +1,13 @@
+//Finding of Pressure of water
+//Given
+W=490500;
+Fr=39240;
+d=40;
+//To Find
+A=(%pi/4)*d^2;
+Wu=W+Fr;
+P1=Wu/A;
+Wd=W-Fr;
+P2=Wd/A;
+disp("Pressure while moving up ="+string(P1)+" N/cm^2");
+disp("Pressure while moving down ="+string(P2)+" N/cm^2");
diff --git a/2021/CH21/EX21.12/Ex21_12.pdf b/2021/CH21/EX21.12/Ex21_12.pdf
new file mode 100755
index 000000000..2080c5a88
--- /dev/null
+++ b/2021/CH21/EX21.12/Ex21_12.pdf
diff --git a/2021/CH21/EX21.13/EX21_13.pdf b/2021/CH21/EX21.13/EX21_13.pdf
new file mode 100755
index 000000000..ec9fd134f
--- /dev/null
+++ b/2021/CH21/EX21.13/EX21_13.pdf
diff --git a/2021/CH21/EX21.13/EX21_13.sce b/2021/CH21/EX21.13/EX21_13.sce
new file mode 100755
index 000000000..5228cdc26
--- /dev/null
+++ b/2021/CH21/EX21.13/EX21_13.sce
@@ -0,0 +1,19 @@
+//Finding of Power delivered
+//Given
+L=(588600-(588600*4)/100);
+d=35;
+l=5;
+q=0.008;
+t1=2.5*60;
+rho=1000;
+g=9.81;
+//To Find
+A=(%pi/4)*d^2;
+P=L/A;
+P1=P*10^4;
+H=P1/(rho*g);
+W1=q*1000*g*H;
+W2=L*l/t1;
+W3=W1+W2;
+W4=W3/1000;
+disp("Power Delivered ="+string(W4)+" Kilo Watts");
diff --git a/2021/CH21/EX21.14/EX21_14.pdf b/2021/CH21/EX21.14/EX21_14.pdf
new file mode 100755
index 000000000..19c9c0d06
--- /dev/null
+++ b/2021/CH21/EX21.14/EX21_14.pdf
diff --git a/2021/CH21/EX21.14/EX21_14.sce b/2021/CH21/EX21.14/EX21_14.sce
new file mode 100755
index 000000000..3c3b73cf6
--- /dev/null
+++ b/2021/CH21/EX21.14/EX21_14.sce
@@ -0,0 +1,10 @@
+//Finding of Efficiency ,Slip
+//Given
+Nb=780;
+Na=800;
+//To Find
+E=Nb/Na;
+E1=E*100;
+S=100-E1;
+disp("Efficiency ="+string(E1)+" percentage");
+disp("Slip ="+string(S)+" Percentage");
diff --git a/2021/CH21/EX21.2/EX21_2.pdf b/2021/CH21/EX21.2/EX21_2.pdf
new file mode 100755
index 000000000..6805f65c6
--- /dev/null
+++ b/2021/CH21/EX21.2/EX21_2.pdf
diff --git a/2021/CH21/EX21.2/EX21_2.sce b/2021/CH21/EX21.2/EX21_2.sce
new file mode 100755
index 000000000..48e64dd48
--- /dev/null
+++ b/2021/CH21/EX21.2/EX21_2.sce
@@ -0,0 +1,17 @@
+//Finding of Force,Power,strokes
+//Given
+d1=0.3;
+d2=0.15;
+W=600;
+d=1.2;
+s=0.25;
+//To Find
+A1=(%pi/4)*d1^2;
+A2=(%pi/4)*d2^2;
+F=(A1/A2)*W;
+W1=W*(d/1200);
+P=W1/1000;
+S=(A1/A2)*(d/s);
+disp("Force ="+string(F)+" Newtons");
+disp("Power required ="+string(P)+" Kilo Watts");
+disp("Number of strokes ="+string(S)+" No units");
diff --git a/2021/CH21/EX21.3/EX21_3.pdf b/2021/CH21/EX21.3/EX21_3.pdf
new file mode 100755
index 000000000..d46e8da75
--- /dev/null
+++ b/2021/CH21/EX21.3/EX21_3.pdf
diff --git a/2021/CH21/EX21.3/EX21_3.sce b/2021/CH21/EX21.3/EX21_3.sce
new file mode 100755
index 000000000..fe96c4eec
--- /dev/null
+++ b/2021/CH21/EX21.3/EX21_3.sce
@@ -0,0 +1,15 @@
+//Finding of Efficiency
+//Given
+W=0.03;
+rho=1000;
+g=9.81;
+w=0.003;
+H1=4;
+H2=18;
+//To Find
+W1=rho*g*W;
+w1=rho*g*w;
+E1=(w1*H2)/(W1*H1)*100;
+E2=(w1*(H2-H1))/((W1-w1)*H1)*100;
+disp("D Aubuissons Efficiency ="+string(E1)+" percentage");
+disp("Rankine Efficiency ="+string(E2)+" percentage");
diff --git a/2021/CH21/EX21.4/EX21_4.pdf b/2021/CH21/EX21.4/EX21_4.pdf
new file mode 100755
index 000000000..71f42721c
--- /dev/null
+++ b/2021/CH21/EX21.4/EX21_4.pdf
diff --git a/2021/CH21/EX21.4/Ex21_4.sce b/2021/CH21/EX21.4/Ex21_4.sce
new file mode 100755
index 000000000..76e659259
--- /dev/null
+++ b/2021/CH21/EX21.4/Ex21_4.sce
@@ -0,0 +1,14 @@
+//Finding of Power, Working ,Idle Period
+//Given
+H=12;
+t=100;
+W=98100;
+v=0.6;
+w=981*12;
+//To Find
+P=w/1000;
+T1=H/v;
+T2=100-T1;
+disp("Power Required ="+string(P)+" Kilo watt");
+disp("Time for working ="+string(T1)+" seconds");
+disp("Idle Time ="+string(T2)+" seconds");
diff --git a/2021/CH21/EX21.7/EX21_7.pdf b/2021/CH21/EX21.7/EX21_7.pdf
new file mode 100755
index 000000000..c601f5673
--- /dev/null
+++ b/2021/CH21/EX21.7/EX21_7.pdf
diff --git a/2021/CH21/EX21.7/EX21_7.sce b/2021/CH21/EX21.7/EX21_7.sce
new file mode 100755
index 000000000..ed92a20d3
--- /dev/null
+++ b/2021/CH21/EX21.7/EX21_7.sce
@@ -0,0 +1,15 @@
+//Finding of Volume ,Diameter
+//Given
+L=44145;
+H=10;
+E=0.55;
+P=490.5*10^4;
+//To Find
+W=L*H;
+F=P*(%pi/4);
+Energy=F*5;
+d=W/(Energy*E);
+d1=d^(.33);
+V=((%pi/4)*d1^2)*(5*d1);
+disp("Diameter ="+string(d1)+" meter");
+disp("Volume ="+string(V)+" m^3");
diff --git a/2021/CH21/EX21.9/EX21_9.pdf b/2021/CH21/EX21.9/EX21_9.pdf
new file mode 100755
index 000000000..2979a5ed8
--- /dev/null
+++ b/2021/CH21/EX21.9/EX21_9.pdf
diff --git a/2021/CH21/EX21.9/EX21_9.sce b/2021/CH21/EX21.9/EX21_9.sce
new file mode 100755
index 000000000..1675c8548
--- /dev/null
+++ b/2021/CH21/EX21.9/EX21_9.sce
@@ -0,0 +1,10 @@
+//Finding of Pressure intensity
+//Given
+P=17;
+D=12;
+d=5;
+//To Find
+A=(%pi/4)*D^2;
+a=(%pi/4)*d^2;
+p=(A/a)*P;
+disp("Pressure Intensity ="+string(p)+" N/cm^2");
diff --git a/2021/CH3/EX3.2/EX3_2.sce b/2021/CH3/EX3.2/EX3_2.sce
new file mode 100755
index 000000000..2311392a3
--- /dev/null
+++ b/2021/CH3/EX3.2/EX3_2.sce
@@ -0,0 +1,13 @@
+//Finding of Total Pressure
+//Given
+d=1.5;
+y1=2;
+rho=1000;
+g=9.81;
+//To Find
+Ig=(%pi*d^4)/64;
+Ay=(%pi/4)*d^2;
+P=Ay*rho*g*y1;
+Ycp=(Ig/Ay)+y1;
+disp("P= "+string(P)+" Newtons");
+disp("Ycp ="+string(Ycp)+" meter");
diff --git a/2021/CH3/EX3.2/Ex3_2.pdf b/2021/CH3/EX3.2/Ex3_2.pdf
new file mode 100755
index 000000000..d2d35243b
--- /dev/null
+++ b/2021/CH3/EX3.2/Ex3_2.pdf
diff --git a/2021/CH3/EX3.3/EX3_3.sce b/2021/CH3/EX3.3/EX3_3.sce
new file mode 100755
index 000000000..bdb6cd25b
--- /dev/null
+++ b/2021/CH3/EX3.3/EX3_3.sce
@@ -0,0 +1,14 @@
+//Finding of Totoal Pressure , Depth of centre
+//given
+d=2.5;
+rho=1000;
+g=9.81;
+y1=2;
+//To Find
+Ig=(%pi*d^4)/64;
+Ay=(%pi/4)*d^2;
+P=Ay*rho*g*y1;
+a=4/6.25;
+Ycp=((Ig*a)/(Ay*y1))+y1;
+disp("P= "+string(P)+" Newtons");
+disp("Ycp ="+string(Ycp)+" meter");
diff --git a/2021/CH3/EX3.3/Ex3_3.pdf b/2021/CH3/EX3.3/Ex3_3.pdf
new file mode 100755
index 000000000..327f6102c
--- /dev/null
+++ b/2021/CH3/EX3.3/Ex3_3.pdf
diff --git a/2021/CH3/EX3.4/EX3_4.sce b/2021/CH3/EX3.4/EX3_4.sce
new file mode 100755
index 000000000..e36e5802e
--- /dev/null
+++ b/2021/CH3/EX3.4/EX3_4.sce
@@ -0,0 +1,15 @@
+//Finding of Total Pressure , Depth of pressure
+//Given
+b=5;
+h=5;
+spgr=0.8;
+rho=800;
+g=9.81;
+y1=(1+(5/3));
+//To Find
+Ig=(b*h^3)/36;
+Ay=(1/2)*b*h;
+P=Ay*rho*g*y1;
+Ycp=(Ig/(Ay*y1))+y1;
+disp("P= "+string(P)+" Newtons");
+disp("Ycp ="+string(Ycp)+" meter");
diff --git a/2021/CH3/EX3.4/Ex3_4.pdf b/2021/CH3/EX3.4/Ex3_4.pdf
new file mode 100755
index 000000000..0a04e0fd5
--- /dev/null
+++ b/2021/CH3/EX3.4/Ex3_4.pdf
diff --git a/2021/CH3/EX3.5/EX3_5.sce b/2021/CH3/EX3.5/EX3_5.sce
new file mode 100755
index 000000000..cdbd61484
--- /dev/null
+++ b/2021/CH3/EX3.5/EX3_5.sce
@@ -0,0 +1,19 @@
+//Finding of Total Pressure , Depth of pressure
+//Given
+T=4;
+rho=1000;
+g=9.81;
+l=2;
+b=1/2;
+y1=2;
+y2=1/3;
+//To Find
+A=(6/2)*1;
+A1=(l*b);
+A2=l*5;
+y3=((A1*y1)+(2*A2*y2))/(A1+2*A2);
+P=rho*g*A*y3;disp(y3);
+Ig=(l^2+(4*l*T)+T^2)/(36*(l+T));
+Ycp=(Ig/(A*y3))+y3;
+disp("P= "+string(P)+" Newtons");
+disp("Ycp ="+string(Ycp)+" meter");
diff --git a/2021/CH3/EX3.5/Ex3_5.pdf b/2021/CH3/EX3.5/Ex3_5.pdf
new file mode 100755
index 000000000..5d5614012
--- /dev/null
+++ b/2021/CH3/EX3.5/Ex3_5.pdf
diff --git a/2021/CH3/EX3.6/EX3_6.sce b/2021/CH3/EX3.6/EX3_6.sce
new file mode 100755
index 000000000..09ff5e97e
--- /dev/null
+++ b/2021/CH3/EX3.6/EX3_6.sce
@@ -0,0 +1,18 @@
+//Finding of Total Pressure , Depth of pressure
+//Given
+spgr=0.9;
+rho=900;
+rho1=1000;
+spgr1=0.6;
+g=9.81;
+y1=spgr*(2/3);
+y2=spgr+(spgr1/2);
+y3=spgr+((spgr1/3)*2);
+//To Find
+P1=rho*g*spgr;
+P2=P1+(rho1*spgr1*g);
+P=(0.5*P1*spgr*1.5)+(((P1+P2)/2)*spgr1*1.5);
+disp("P ="+string(P)+" Newton");
+P3=P2-P1;
+Ycp=((P1*y1)+(P2*y2)+(P3*y3))/P;
+disp("Ycp ="+string(Ycp)+" meter");
diff --git a/2021/CH3/EX3.6/Ex3_6.pdf b/2021/CH3/EX3.6/Ex3_6.pdf
new file mode 100755
index 000000000..49f669ec3
--- /dev/null
+++ b/2021/CH3/EX3.6/Ex3_6.pdf
diff --git a/2021/CH3/EX3.7/EX3_7.sce b/2021/CH3/EX3.7/EX3_7.sce
new file mode 100755
index 000000000..f6ddaa16c
--- /dev/null
+++ b/2021/CH3/EX3.7/EX3_7.sce
@@ -0,0 +1,16 @@
+//Finding of Total Pressure
+//Given
+BC=2;
+d=2;
+y1=2.5;
+rho=1000;
+g=9.81;
+//To Find
+Ig=(1*BC^3)/2;
+Ay=((1*BC^3)/2)*y1;
+Px=Ay*rho*g/2;
+Ycp=(Ig/Ay)+y1;
+Py=((2*1.5)*(%pi/4)*d^2)*rho*g;
+disp("Px= "+string(Px)+" Newtons");
+disp("Ycp ="+string(Ycp)+" meter");
+disp("Py= "+string(Py)+" Newtons");
diff --git a/2021/CH3/EX3.7/Ex3_7.pdf b/2021/CH3/EX3.7/Ex3_7.pdf
new file mode 100755
index 000000000..139ab97ab
--- /dev/null
+++ b/2021/CH3/EX3.7/Ex3_7.pdf
diff --git a/2021/CH4/EX4.1/EX4_1.sce b/2021/CH4/EX4.1/EX4_1.sce
new file mode 100755
index 000000000..c7cdd5cb8
--- /dev/null
+++ b/2021/CH4/EX4.1/EX4_1.sce
@@ -0,0 +1,17 @@
+//Finding of water displaced and position of centre buoyancy
+//Given
+l=2;
+h=1.5;
+b=4;
+v=l*b*h;
+spgr=0.7;
+rho=700;
+rho1=1000;
+g=9.81;
+w=rho*g*v;
+//To Find 
+wd=w/(rho1*g);
+disp("Water Displaced = "+string(wd)+" m^3");
+h1=wd/(l*b);
+h2=h1/2;
+disp("Positin of Centre Buoyancy = "+string(h2)+" m");
diff --git a/2021/CH4/EX4.1/Ex4_1.pdf b/2021/CH4/EX4.1/Ex4_1.pdf
new file mode 100755
index 000000000..516509cf3
--- /dev/null
+++ b/2021/CH4/EX4.1/Ex4_1.pdf
diff --git a/2021/CH4/EX4.3/EX4_3.sce b/2021/CH4/EX4.3/EX4_3.sce
new file mode 100755
index 000000000..f03eceb1b
--- /dev/null
+++ b/2021/CH4/EX4.3/EX4_3.sce
@@ -0,0 +1,14 @@
+//Finding of volume and specific gravity
+//given
+w=490.5;             //In Air in Newton
+w1=196.2;            //In Water in Newton
+rho=1000;
+g=9.81;
+rho1=5000;
+//To Find 
+wd=w-w1;
+vd=wd/(rho*g);
+rho2=(w/g)/vd;
+spgr=rho1/(rho*3);
+disp("Volume = "+string(rho2)+" Kg/m^3");
+disp("Specific Gravity  is= "+string(spgr)+" No units")
diff --git a/2021/CH4/EX4.3/Ex4_3.pdf b/2021/CH4/EX4.3/Ex4_3.pdf
new file mode 100755
index 000000000..0ea14f9a4
--- /dev/null
+++ b/2021/CH4/EX4.3/Ex4_3.pdf
diff --git a/2021/CH4/EX4.4/EX4_4.sce b/2021/CH4/EX4.4/EX4_4.sce
new file mode 100755
index 000000000..14e9ae145
--- /dev/null
+++ b/2021/CH4/EX4.4/EX4_4.sce
@@ -0,0 +1,14 @@
+//Finding of Mass,Density,Specific Gravity
+//Given
+v=2*1*3;
+w=3924;
+rho=1000;
+g=9.81;
+wd=rho*g*v;
+w1=w+wd;
+m=w1/g;
+rho1=m/v;
+spgr=rho1/rho;
+disp("Mass is = "+string(m)+" Kg");
+disp("Density is = "+string(rho1)+" Kg/m^3");
+disp("Specific Gravity = "+string(spgr)+" No units")
diff --git a/2021/CH4/EX4.4/Ex4_4.pdf b/2021/CH4/EX4.4/Ex4_4.pdf
new file mode 100755
index 000000000..c75e9c06b
--- /dev/null
+++ b/2021/CH4/EX4.4/Ex4_4.pdf
diff --git a/2021/CH4/EX4.5/EX4_5.sce b/2021/CH4/EX4.5/EX4_5.sce
new file mode 100755
index 000000000..5660bbd67
--- /dev/null
+++ b/2021/CH4/EX4.5/EX4_5.sce
@@ -0,0 +1,11 @@
+//Finding of Density
+//Given
+h1=0.4;
+h2=0.6;
+rho=1000;
+rho1=13600;
+g=9.81;
+wd=rho*0.6;
+md=rho1*0.4;
+rho2=wd+md;
+disp("Density is = "+string(rho2)+" Kg/m^3");
diff --git a/2021/CH4/EX4.5/Ex4_5.pdf b/2021/CH4/EX4.5/Ex4_5.pdf
new file mode 100755
index 000000000..4a9c3b4a5
--- /dev/null
+++ b/2021/CH4/EX4.5/Ex4_5.pdf
diff --git a/2021/CH4/EX4.6/EX4_6.pdf b/2021/CH4/EX4.6/EX4_6.pdf
new file mode 100755
index 000000000..aaa1ed49a
--- /dev/null
+++ b/2021/CH4/EX4.6/EX4_6.pdf
diff --git a/2021/CH4/EX4.6/EX4_6.sce b/2021/CH4/EX4.6/EX4_6.sce
new file mode 100755
index 000000000..41d34e557
--- /dev/null
+++ b/2021/CH4/EX4.6/EX4_6.sce
@@ -0,0 +1,18 @@
+//Finding of  Weight and Metacentric height
+//Given
+l=4;
+b=2;
+h=1;
+d=0.6;
+v=l*b*d;
+rho=1000;
+g=9.81;
+//To Find
+wd=rho*g*v;
+disp("Weight of the body ="+string(wd)+" Newtons");
+I=(l*b^3)/12;
+h1=h/2;
+d1=d/2;
+h2=h1-d1;
+mh=(I/v)-h2;
+disp("Metacentric Height ="+string(mh)+" meter");
diff --git a/2021/CH4/EX4.7/EX4_7.pdf b/2021/CH4/EX4.7/EX4_7.pdf
new file mode 100755
index 000000000..6b2f55ebf
--- /dev/null
+++ b/2021/CH4/EX4.7/EX4_7.pdf
diff --git a/2021/CH4/EX4.7/EX4_7.sce b/2021/CH4/EX4.7/EX4_7.sce
new file mode 100755
index 000000000..c8f469ba8
--- /dev/null
+++ b/2021/CH4/EX4.7/EX4_7.sce
@@ -0,0 +1,15 @@
+//Finding pf Metacentric Height
+//Given
+d=3;
+h=2;
+spgr=0.7;
+h1=h*spgr;
+pi=3.14;
+//To Find
+h2=h/2;
+h3=h1/2;
+h4=h2-h3;
+mh=(pi*d^4)/64;
+vwd=(pi*d^2*h1)/4;
+mg=(mh/vwd)-h4;
+disp("Metacentric Height is ="+string(mg)+" meter");
diff --git a/2021/CH5/EX5.1/EX5_1.pdf b/2021/CH5/EX5.1/EX5_1.pdf
new file mode 100755
index 000000000..58d39cc7a
--- /dev/null
+++ b/2021/CH5/EX5.1/EX5_1.pdf
diff --git a/2021/CH5/EX5.1/EX5_1.sce b/2021/CH5/EX5.1/EX5_1.sce
new file mode 100755
index 000000000..df65f5113
--- /dev/null
+++ b/2021/CH5/EX5.1/EX5_1.sce
@@ -0,0 +1,15 @@
+//Finding of velocity and discharge
+//Given
+d1=0.4;
+r1=d1/2;
+d2=0.2;
+r2=d2/2;
+v1=5;
+pi=3.14;
+//To Find
+a1=(pi*r1^2);
+a2=(pi*r2^2);
+v2=(a1*v1)/a2;
+q2=a2*v2;
+disp("Velocity at section -2 ="+string(v2)+" m/second");
+disp("Discharge ="+string(q2)+"m^3/seconds");
diff --git a/2021/CH5/EX5.10/EX5_10.pdf b/2021/CH5/EX5.10/EX5_10.pdf
new file mode 100755
index 000000000..63da51fbe
--- /dev/null
+++ b/2021/CH5/EX5.10/EX5_10.pdf
diff --git a/2021/CH5/EX5.10/EX5_10.sce b/2021/CH5/EX5.10/EX5_10.sce
new file mode 100755
index 000000000..9342c533b
--- /dev/null
+++ b/2021/CH5/EX5.10/EX5_10.sce
@@ -0,0 +1,11 @@
+//Finding of convective acceleration
+//Given
+v1=2.5;
+v2=16;
+s=3.75;
+//To Find
+a=(v2-v1)/s;
+a1=v1*a;
+a2=v2*a;
+disp("Acceleration at inlet="+string(a1)+" m/s^2");
+disp("Acceleration at outlet="+string(a2)+" m/s^2");
diff --git a/2021/CH5/EX5.2/EX5_2.pdf b/2021/CH5/EX5.2/EX5_2.pdf
new file mode 100755
index 000000000..74b64fe70
--- /dev/null
+++ b/2021/CH5/EX5.2/EX5_2.pdf
diff --git a/2021/CH5/EX5.2/EX5_2.sce b/2021/CH5/EX5.2/EX5_2.sce
new file mode 100755
index 000000000..ffbaca565
--- /dev/null
+++ b/2021/CH5/EX5.2/EX5_2.sce
@@ -0,0 +1,15 @@
+//Finding of Discharge and velocity
+//Given
+d1=0.4;
+d2=0.3;
+d3=0.2;
+pi=3.14;
+//To Find
+q1=(pi/4)*d1^2*3;
+q2=(pi/4)*d2^2*2;
+q3=q1-q2;
+v3=q3/((pi/4)*d3^2);
+disp("Discharge at section- 1="+string(q1)+" m^3/sec");
+disp("Discharge at section- 2="+string(q2)+" m^3/sec");
+disp("Discharge at section- 3="+string(q3)+" m^3/sec");
+disp("velocity at section- 3 ="+string(v3)+" m/sec");
diff --git a/2021/CH6/EX6.3/EX6_3.pdf b/2021/CH6/EX6.3/EX6_3.pdf
new file mode 100755
index 000000000..35e7bbee1
--- /dev/null
+++ b/2021/CH6/EX6.3/EX6_3.pdf
diff --git a/2021/CH6/EX6.3/EX6_3.sce b/2021/CH6/EX6.3/EX6_3.sce
new file mode 100755
index 000000000..6d5b240e7
--- /dev/null
+++ b/2021/CH6/EX6.3/EX6_3.sce
@@ -0,0 +1,15 @@
+//Finding of Discharge
+//Given
+d1=30;
+d2=15;
+hom=10;
+cod=0.98;
+pi=3.14;
+g=9.81;
+//To Find
+a1=(pi/4)*d1^2;
+a2=(pi/4)*d2^2;
+h=hom*(12.6);
+q=(cod*a1*a2*(2*100*g*h)^(1/2))/((a1^2-a2^2)^(1/2));
+q1=q/1000;
+disp("Discharge ="+string(q1)+" m^3/sec");
diff --git a/2021/CH6/EX6.4/EX6_4.pdf b/2021/CH6/EX6.4/EX6_4.pdf
new file mode 100755
index 000000000..749f83173
--- /dev/null
+++ b/2021/CH6/EX6.4/EX6_4.pdf
diff --git a/2021/CH6/EX6.4/EX6_4.sce b/2021/CH6/EX6.4/EX6_4.sce
new file mode 100755
index 000000000..ca8bd8cf2
--- /dev/null
+++ b/2021/CH6/EX6.4/EX6_4.sce
@@ -0,0 +1,15 @@
+//Finding of Discharge
+//Given
+d1=30;
+d2=15;
+hom=30;
+cod=0.98;
+g=9.81;
+pi=3.14;
+//To Find 
+a1=(pi/4)*d1^2;
+a2=(pi/4)*d2^2;
+h=hom*(12.6);
+q=(cod*a1*a2*(2*100*g*h)^(1/2))/((a1^2-a2^2)^(1/2));
+q1=q/1000;
+disp("Discharge ="+string(q1)+" m^3/sec");
diff --git a/2021/CH6/EX6.5/EX6_5.pdf b/2021/CH6/EX6.5/EX6_5.pdf
new file mode 100755
index 000000000..e4a2a5c9f
--- /dev/null
+++ b/2021/CH6/EX6.5/EX6_5.pdf
diff --git a/2021/CH6/EX6.5/EX6_5.sce b/2021/CH6/EX6.5/EX6_5.sce
new file mode 100755
index 000000000..31f5b34af
--- /dev/null
+++ b/2021/CH6/EX6.5/EX6_5.sce
@@ -0,0 +1,15 @@
+//Finding of Rate of flow
+//Given
+d1=30;
+d2=15;
+hom=30;
+cod=0.98;
+g=9.81;
+pi=3.14;
+//To Find 
+a1=(pi/4)*d1^2;
+a2=(pi/4)*d2^2;
+h=hom*((13.6/0.8)-1);
+q=(cod*a1*a2*(2*100*g*h)^(1/2))/((a1^2-a2^2)^(1/2));
+q1=q/1000;
+disp("Discharge ="+string(q1)+" m^3/sec");
diff --git a/2021/CH6/EX6.6/EX6_6.pdf b/2021/CH6/EX6.6/EX6_6.pdf
new file mode 100755
index 000000000..631c533e7
--- /dev/null
+++ b/2021/CH6/EX6.6/EX6_6.pdf
diff --git a/2021/CH6/EX6.6/EX6_6.sce b/2021/CH6/EX6.6/EX6_6.sce
new file mode 100755
index 000000000..62928e77d
--- /dev/null
+++ b/2021/CH6/EX6.6/EX6_6.sce
@@ -0,0 +1,16 @@
+//Finding of Pressure Difference
+//Given
+g=9.81;
+spgr=0.9;
+spgr1=13.6;
+rho=1000;
+rho1=spgr*1000*g;
+zd=0.3;
+gd=25;
+x=(spgr1/spgr)-1;
+x1=((gd*x)/100)+zd;
+//To find
+pd=x1*rho1;disp(x1);
+disp("Pressure Difference ="+string(pd)+" N/m^2");
+pd1=pd/10000;
+disp("Pressure Difference ="+string(pd1)+" N/cm^2");
diff --git a/2021/CH6/EX6.7/EX6_7.pdf b/2021/CH6/EX6.7/EX6_7.pdf
new file mode 100755
index 000000000..a9a678c83
--- /dev/null
+++ b/2021/CH6/EX6.7/EX6_7.pdf
diff --git a/2021/CH6/EX6.7/EX6_7.sce b/2021/CH6/EX6.7/EX6_7.sce
new file mode 100755
index 000000000..1d3da7b72
--- /dev/null
+++ b/2021/CH6/EX6.7/EX6_7.sce
@@ -0,0 +1,19 @@
+//Finding of Rate Of Flow
+//Given
+d1=15;
+d2=30;
+hm=50;
+spgr=0.9;
+spgr1=13.6;
+cod=0.64;
+g=9.81;
+//To Find
+A0=(%pi/4)*d1^2;
+A1=(%pi/4)*d2^2;
+h=((spgr1/spgr)-1)*hm;
+x=(A0*A1)/sqrt(A1^2-A0^2);
+y=sqrt(2*g*h);
+q=cod*x*y;
+disp(" Discharge ="+string(q)+" cm^3/sec");
+q1=q/100;
+disp(" Discharge ="+string(q1)+" ltr/sec");
diff --git a/2021/CH6/EX6.8/EX6_8.pdf b/2021/CH6/EX6.8/EX6_8.pdf
new file mode 100755
index 000000000..596096575
--- /dev/null
+++ b/2021/CH6/EX6.8/EX6_8.pdf
diff --git a/2021/CH6/EX6.8/EX6_8.sce b/2021/CH6/EX6.8/EX6_8.sce
new file mode 100755
index 000000000..d25d934d5
--- /dev/null
+++ b/2021/CH6/EX6.8/EX6_8.sce
@@ -0,0 +1,12 @@
+//Fiding of Discharge
+//Given
+d1=0.3;
+pd=0.06;
+g=9.81;
+cv=0.98;
+//To Find
+vc=sqrt(2*g*pd)*cv;
+V=0.8*vc;
+A=(%pi/4)*d1^2;
+q=V*A;
+disp(" Discharge ="+string(q)+" m^3/sec");
diff --git a/2021/CH7/EX7.1/EX7_1.pdf b/2021/CH7/EX7.1/EX7_1.pdf
new file mode 100755
index 000000000..5da6d1a72
--- /dev/null
+++ b/2021/CH7/EX7.1/EX7_1.pdf
diff --git a/2021/CH7/EX7.1/EX7_1.sce b/2021/CH7/EX7.1/EX7_1.sce
new file mode 100755
index 000000000..23f790555
--- /dev/null
+++ b/2021/CH7/EX7.1/EX7_1.sce
@@ -0,0 +1,14 @@
+//Finding of Loss of Head
+//Given
+q1=200;
+d1=150;
+d2=300;
+g=9.81;
+//To Find
+v1=200*(4/%pi)*(100/150)^2;
+disp(v1);
+v2=200*(4/%pi)*(100/300)^2;
+disp(v2);
+h=((v1-v2)^2)/20*g;
+h1=h/1000;
+disp(" Loss of Head ="+string(h1)+" meter of water"); 
diff --git a/2021/CH7/EX7.12/EX7_12.pdf b/2021/CH7/EX7.12/EX7_12.pdf
new file mode 100755
index 000000000..8a5627ad3
--- /dev/null
+++ b/2021/CH7/EX7.12/EX7_12.pdf
diff --git a/2021/CH7/EX7.12/EX7_12.sce b/2021/CH7/EX7.12/EX7_12.sce
new file mode 100755
index 000000000..939addd05
--- /dev/null
+++ b/2021/CH7/EX7.12/EX7_12.sce
@@ -0,0 +1,18 @@
+//Finding of Maximum Power Outlet
+//Given
+d=0.4;
+l=400;
+H=420;
+rho=1000;
+f=0.025;
+g=9.81;
+pi=3.14;
+//To Find
+h=H/3;
+h1=(f*l*100)/(2*g*d);
+v=sqrt(h/h1);disp(h); disp(h1);
+a=(pi/4)*d^2;
+q=a*v;
+h3=H-h;disp(h3);
+p=(rho*g*q*h3)/1000;
+disp(" Maximum Power Outlet ="+string(p)+" KW");
diff --git a/2021/CH7/EX7.13/EX7_13.pdf b/2021/CH7/EX7.13/EX7_13.pdf
new file mode 100755
index 000000000..c30764e1c
--- /dev/null
+++ b/2021/CH7/EX7.13/EX7_13.pdf
diff --git a/2021/CH7/EX7.13/Ex7_13.sce b/2021/CH7/EX7.13/Ex7_13.sce
new file mode 100755
index 000000000..75abb0290
--- /dev/null
+++ b/2021/CH7/EX7.13/Ex7_13.sce
@@ -0,0 +1,12 @@
+//Finding of Rise of Pressure
+//Given
+l=2000;
+d=0.6;
+v=2;
+c=1420;
+t=20;
+rho=1000;
+//To Find
+p=(rho*l*v)/t
+p1=p/10000;
+disp("Rise of Pressure ="+string(p1)+"N/cm^2");
diff --git a/2021/CH7/EX7.14/EX7_14.pdf b/2021/CH7/EX7.14/EX7_14.pdf
new file mode 100755
index 000000000..f20c844af
--- /dev/null
+++ b/2021/CH7/EX7.14/EX7_14.pdf
diff --git a/2021/CH7/EX7.14/Ex7_14.sce b/2021/CH7/EX7.14/Ex7_14.sce
new file mode 100755
index 000000000..6a62c0219
--- /dev/null
+++ b/2021/CH7/EX7.14/Ex7_14.sce
@@ -0,0 +1,13 @@
+//Finding of Rise of Pressure
+//Given
+l=2000;
+d=0.6;
+v=2;
+t=20;
+k=19.62*10^8;
+rho=1000;
+//To Find
+c=sqrt(k/rho);
+p=rho*v*c
+p1=p/10000;
+disp("Rise of Pressure ="+string(p1)+"N/cm^2");
diff --git a/2021/CH7/EX7.15/EX7_15.pdf b/2021/CH7/EX7.15/EX7_15.pdf
new file mode 100755
index 000000000..ea4abac31
--- /dev/null
+++ b/2021/CH7/EX7.15/EX7_15.pdf
diff --git a/2021/CH7/EX7.15/Ex7_15.sce b/2021/CH7/EX7.15/Ex7_15.sce
new file mode 100755
index 000000000..4e46fe062
--- /dev/null
+++ b/2021/CH7/EX7.15/Ex7_15.sce
@@ -0,0 +1,20 @@
+//Finding of Rise of Pressure and circumferential, Longitudinal stress
+//Given
+l=2000;
+d=0.6;
+v=2;
+t=0.01;
+k=19.62*10^8;
+rho=1000;
+E=19.62*10^10;
+//To Find
+a=(1/k)+(d/(t*E));
+b=(1/rho)*a;
+c=sqrt(b);
+p=2/c;
+p1=p/10000;
+fc=((p*d)/(2*t))/10000;
+fl=((p*d)/(4*t))/10000;
+disp("Rise of Pressure ="+string(p1)+"N/cm^2");
+disp("Circumferential stress ="+string(fc)+" N/m^2");
+disp("Longitudinal stress ="+string(fl)+" N/m^2");
diff --git a/2021/CH7/EX7.16/EX7_16.pdf b/2021/CH7/EX7.16/EX7_16.pdf
new file mode 100755
index 000000000..88f837f59
--- /dev/null
+++ b/2021/CH7/EX7.16/EX7_16.pdf
diff --git a/2021/CH7/EX7.16/Ex7_16.sce b/2021/CH7/EX7.16/Ex7_16.sce
new file mode 100755
index 000000000..b3e1e4c0c
--- /dev/null
+++ b/2021/CH7/EX7.16/Ex7_16.sce
@@ -0,0 +1,21 @@
+//Finding of Maximum rise of water level,velocity,Time of occurance
+//Given
+d=4;
+d1=1;
+l=150;
+q=2;
+g=9.81;
+//To Find
+a1=(%pi/4)*d^2;
+a2=(%pi/4)*d1^2;
+v=q/a2;
+a=(l*a2)/(g*a1);
+b=sqrt(a);
+h=v*b;
+c=(l*a1)/(g*a2);
+d=sqrt(c);
+t=(%pi/2)*d;
+v1=v*(a2/a1);
+disp("Maximum rise of water ="+string(h)+" meter");
+disp("Time taken ="+string(t)+" seconds");
+disp("Maximum Velocity ="+string(v1)+" m/sec");
diff --git a/2021/CH7/EX7.3/EX7_3.pdf b/2021/CH7/EX7.3/EX7_3.pdf
new file mode 100755
index 000000000..4b67949ab
--- /dev/null
+++ b/2021/CH7/EX7.3/EX7_3.pdf
diff --git a/2021/CH7/EX7.3/EX7_3.sce b/2021/CH7/EX7.3/EX7_3.sce
new file mode 100755
index 000000000..c4d3cd358
--- /dev/null
+++ b/2021/CH7/EX7.3/EX7_3.sce
@@ -0,0 +1,14 @@
+//Finding of Discharge
+//Given
+d=0.3;
+l=400;
+f=0.00225;
+h=5;
+g=9.81;
+//To Find
+x=(h*2*g*d)/(f*l);
+v=sqrt(x);
+disp(v);
+A=(%pi/4)*d^2;
+q=A*v;disp(A);
+disp(" Discharge ="+string(q)+" m^3/sec");
diff --git a/2021/CH7/EX7.6/EX7_6.pdf b/2021/CH7/EX7.6/EX7_6.pdf
new file mode 100755
index 000000000..692b3c4d1
--- /dev/null
+++ b/2021/CH7/EX7.6/EX7_6.pdf
diff --git a/2021/CH7/EX7.6/EX7_6.sce b/2021/CH7/EX7.6/EX7_6.sce
new file mode 100755
index 000000000..7b091d6b9
--- /dev/null
+++ b/2021/CH7/EX7.6/EX7_6.sce
@@ -0,0 +1,15 @@
+//Finding of Head
+//Giiven
+f=0.032;
+l=400;
+d=0.3;
+q=0.3;
+g=9.81;
+//TO find
+A=(%pi/4)*d^2
+V=q/A;
+v1=(V^2);
+x=1.5+(f*l/d);
+y=v1/(2*g);
+H=x*y;
+disp("Difference in water level ="+string(H)+" meter")
diff --git a/2021/CH7/EX7.9/EX7_9.pdf b/2021/CH7/EX7.9/EX7_9.pdf
new file mode 100755
index 000000000..815cc7d47
--- /dev/null
+++ b/2021/CH7/EX7.9/EX7_9.pdf
diff --git a/2021/CH7/EX7.9/EX7_9.sce b/2021/CH7/EX7.9/EX7_9.sce
new file mode 100755
index 000000000..cea9fbe23
--- /dev/null
+++ b/2021/CH7/EX7.9/EX7_9.sce
@@ -0,0 +1,17 @@
+//Finding of Equivalent Diameter
+//Given
+L=1400;
+L1=800;
+L2=400;
+L3=200;
+D1=0.6;
+D2=0.4;
+D3=0.2;
+//To Find
+a=L1/(D1)^5;disp(a);
+b=L2/(D2)^5;disp(b);
+c=L3/(D3)^5;disp(c);
+d=(a+b+c);disp(d);
+d1=d^1/5;
+D=L/d1;
+disp( "Diameter ="+string(D)+" meter");
diff --git a/2021/CH8/EX8.1/EX8_1.pdf b/2021/CH8/EX8.1/EX8_1.pdf
new file mode 100755
index 000000000..7a03e5d93
--- /dev/null
+++ b/2021/CH8/EX8.1/EX8_1.pdf
diff --git a/2021/CH8/EX8.1/Ex8_1.sce b/2021/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..56980043f
--- /dev/null
+++ b/2021/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,14 @@
+//Finding of Actual Discharge,velocity
+//Given
+d=0.05;
+H=12;
+Cd=0.6;
+Cv=0.98;
+g=9.81;
+//To Find
+a=(%pi/4)*d^2;
+v=sqrt(2*g*H);
+q=Cd*a*v;
+V=Cv*v;
+disp("Actual Discharge ="+string(q)+" m^3/sec");
+disp("Actual Velocity ="+string(V)+" m/sec");
diff --git a/2021/CH8/EX8.2/EX8_2.pdf b/2021/CH8/EX8.2/EX8_2.pdf
new file mode 100755
index 000000000..ce7c90457
--- /dev/null
+++ b/2021/CH8/EX8.2/EX8_2.pdf
diff --git a/2021/CH8/EX8.2/Ex8_2.sce b/2021/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..dc4f2912a
--- /dev/null
+++ b/2021/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,12 @@
+//Finding of Coefficient of Discharge
+//Given
+d=0.03;
+H=1.5;
+Ad=2.35*10^-3;
+g=9.81;
+//To Find
+a=(%pi/4)*d^2;
+b=sqrt(2*g*H);
+Td=b*a;
+Cd=Ad/Td;
+disp(" Coefficient of Discharge ="+string(Cd)+" No units");
diff --git a/2021/CH8/EX8.3/EX8_3.pdf b/2021/CH8/EX8.3/EX8_3.pdf
new file mode 100755
index 000000000..009f1cf76
--- /dev/null
+++ b/2021/CH8/EX8.3/EX8_3.pdf
diff --git a/2021/CH8/EX8.3/Ex8_3.sce b/2021/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..ebe1fa600
--- /dev/null
+++ b/2021/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,12 @@
+//Finding of Cv,Cc
+//Given
+H=0.6;
+x=0.1;
+y=0.0045;
+Cd=0.6;
+//To Find
+a=sqrt(4*y*H);
+Cv=x/a;
+Cc=Cd/Cv;
+disp("Cv ="+string(Cv)+" No units");
+disp("Cc ="+string(Cc)+" No units");
diff --git a/2021/CH8/EX8.4/EX8_4.pdf b/2021/CH8/EX8.4/EX8_4.pdf
new file mode 100755
index 000000000..1e081a7a0
--- /dev/null
+++ b/2021/CH8/EX8.4/EX8_4.pdf
diff --git a/2021/CH8/EX8.4/Ex8_4.sce b/2021/CH8/EX8.4/Ex8_4.sce
new file mode 100755
index 000000000..7d2eaeecc
--- /dev/null
+++ b/2021/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,24 @@
+//Finding of Cv,Cd,Cc
+//Given
+H=5;
+d1=0.1;
+d2=2;
+t=30;
+h=0.45;
+x=1;
+g=9.81;
+y=0.052;
+H=5;
+//To Find
+A1=(%pi/4)*d2^2;
+Aq=(A1*h)/t;
+A2=(%pi/4)*d1^2;
+b=sqrt(2*g*H);
+Tq=A2*b;
+Cd=Aq/Tq;
+c=sqrt(4*y*H);
+Cv=x/c;
+Cc=Cd/Cv;
+disp(" Cd ="+string(Cd)+" No units");
+disp(" Cv ="+string(Cv)+" No units");
+disp(" Cc ="+string(Cc)+" No units");
diff --git a/2021/CH8/EX8.6/EX8_6.pdf b/2021/CH8/EX8.6/EX8_6.pdf
new file mode 100755
index 000000000..2d0fa20f9
--- /dev/null
+++ b/2021/CH8/EX8.6/EX8_6.pdf
diff --git a/2021/CH8/EX8.6/Ex8_6.sce b/2021/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..7211abd8f
--- /dev/null
+++ b/2021/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,11 @@
+//Finding of Dischage through Rectangular orifice
+//Given
+H1=4;
+H2=6;
+Cd=0.62;
+g=9.81;
+//To Find
+a=H2^(3/2)-H1^(3/2);
+b=sqrt(2*g);
+q=Cd*2*b*a;
+disp("Dischage through Rectangular orifice ="+string(q)+" m^3/sec");
diff --git a/2021/CH8/EX8.7/EX8_7.pdf b/2021/CH8/EX8.7/EX8_7.pdf
new file mode 100755
index 000000000..3dc8d3dc0
--- /dev/null
+++ b/2021/CH8/EX8.7/EX8_7.pdf
diff --git a/2021/CH8/EX8.7/Ex8_7.sce b/2021/CH8/EX8.7/Ex8_7.sce
new file mode 100755
index 000000000..0efead319
--- /dev/null
+++ b/2021/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,13 @@
+//Finding of Dischage through a fully submersed Orifice
+//Given
+b=2;
+H=0.8;
+H1=2.5;
+H2=3;
+Cd=0.6;
+g=9.81;
+//To Find
+a=sqrt(2*g*H);
+q=Cd*a*b*(H2-H1);
+disp("Dischage through a fully submersed Orifice ="+string(q)+" m^3/sec");
+
diff --git a/2021/CH8/EX8.8/EX8_8.pdf b/2021/CH8/EX8.8/EX8_8.pdf
new file mode 100755
index 000000000..59a81523f
--- /dev/null
+++ b/2021/CH8/EX8.8/EX8_8.pdf
diff --git a/2021/CH8/EX8.8/Ex8_8.sce b/2021/CH8/EX8.8/Ex8_8.sce
new file mode 100755
index 000000000..49ba60776
--- /dev/null
+++ b/2021/CH8/EX8.8/Ex8_8.sce
@@ -0,0 +1,14 @@
+//Dischage through Orifice
+//Given
+b=1.5;
+H1=3.2;
+H2=2;
+H3=2.4;
+Cd=0.62;
+g=9.81;
+//To Find
+a=H3^(3/2)-H2^(3/2);
+q1=(2/3)*Cd*b*sqrt(2*g)*a;
+q2=Cd*b*(H1-H3)*sqrt(2*g*H3);
+q3=q1+q2;
+disp("Dischage through Orifice ="+string(q3)+" m^3/sec");
diff --git a/2021/CH8/EX8.9/EX8_9.pdf b/2021/CH8/EX8.9/EX8_9.pdf
new file mode 100755
index 000000000..bd15bc9a6
--- /dev/null
+++ b/2021/CH8/EX8.9/EX8_9.pdf
diff --git a/2021/CH8/EX8.9/Ex8_9.sce b/2021/CH8/EX8.9/Ex8_9.sce
new file mode 100755
index 000000000..bababab8e
--- /dev/null
+++ b/2021/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,18 @@
+//Finding of Time Taken
+//Given
+d1=3;
+d2=0.4;
+H1=4;
+H2=2;
+g=9.81;
+Cd=0.6;
+//To Find
+A=(%pi/4)*d1^2;
+a=(%pi/4)*d2^2;
+//To empty from 4-2 meter
+c=sqrt(H1)-sqrt(H2);
+T=(2*A*c)/(Cd*a*sqrt(2*g));
+disp("Time Taken To empty from 4 to 2 meter ="+string(T)+" seconds");
+//To empty the tank fully
+T1=(2*A*sqrt(H1))/(Cd*a*sqrt(2*g));
+disp("Time Taken To empty the tank ="+string(T1)+" seconds");
diff --git a/2021/CH9/EX9.1/EX9_1.pdf b/2021/CH9/EX9.1/EX9_1.pdf
new file mode 100755
index 000000000..2dad8565d
--- /dev/null
+++ b/2021/CH9/EX9.1/EX9_1.pdf
diff --git a/2021/CH9/EX9.1/Ex9_1.sce b/2021/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..079694c33
--- /dev/null
+++ b/2021/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,9 @@
+//Finding of Discharge through rectangular Notch
+//Given
+H=0.4;
+L=3;
+Cd=0.6;
+g=9.81;
+//To Find
+q=(2/3)*Cd*L*sqrt(2*g)*H^(3/2);
+disp("Discharge through rectangular Notch ="+string(q)+" m^3/sec");
diff --git a/2021/CH9/EX9.10/EX9_10.pdf b/2021/CH9/EX9.10/EX9_10.pdf
new file mode 100755
index 000000000..35d5c27fc
--- /dev/null
+++ b/2021/CH9/EX9.10/EX9_10.pdf
diff --git a/2021/CH9/EX9.10/Ex9_10.sce b/2021/CH9/EX9.10/Ex9_10.sce
new file mode 100755
index 000000000..fc6d9ff5c
--- /dev/null
+++ b/2021/CH9/EX9.10/Ex9_10.sce
@@ -0,0 +1,8 @@
+//Finding of Discharge
+//Given
+L=5.4;
+n=6;
+H=0.45;
+//To Find
+q=1.84*(L-(0.1*n*H))*H^(3/2);
+disp("Discharge ="+string(q)+" m^3/sec");
diff --git a/2021/CH9/EX9.2/EX9_2.pdf b/2021/CH9/EX9.2/EX9_2.pdf
new file mode 100755
index 000000000..f5146796c
--- /dev/null
+++ b/2021/CH9/EX9.2/EX9_2.pdf
diff --git a/2021/CH9/EX9.2/Ex9_2.sce b/2021/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..fc440b5e2
--- /dev/null
+++ b/2021/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,11 @@
+//Finding of Height
+//Given
+q=1.5;
+Cd=0.6;
+L=5;
+g=9.81;
+//To Find
+H=q/((2/3)*Cd*L*sqrt(2*g));
+H1=H^(2/3);
+Z=q-H1;disp(H1);
+disp("Height ="+string(Z)+" meter");
diff --git a/2021/CH9/EX9.3/EX9_3.pdf b/2021/CH9/EX9.3/EX9_3.pdf
new file mode 100755
index 000000000..de399b15e
--- /dev/null
+++ b/2021/CH9/EX9.3/EX9_3.pdf
diff --git a/2021/CH9/EX9.3/Ex9_3.sce b/2021/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..57d55cd59
--- /dev/null
+++ b/2021/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,13 @@
+//Finding of Position of Apex of Notch
+//Given
+q=0.20;
+d=1;
+theta=90;
+Cd=0.62;
+g=9.81;
+//To Find
+b=(theta/2);
+H=q/((8/15)*Cd*sqrt(2*g)*tan(b));
+H1=H^(2/5);
+p=d-H1;
+disp("Position of Apex of Notch ="+string(p)+" meter");
diff --git a/2021/CH9/EX9.4/EX9_4.pdf b/2021/CH9/EX9.4/EX9_4.pdf
new file mode 100755
index 000000000..7d0d63ad8
--- /dev/null
+++ b/2021/CH9/EX9.4/EX9_4.pdf
diff --git a/2021/CH9/EX9.4/Ex9_4.sce b/2021/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..014783dce
--- /dev/null
+++ b/2021/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,16 @@
+//Finding of discharge through Trapezoidal Notch
+//Given
+H=0.3;
+Cd1=0.62;
+Cd2=0.6;
+d=0.4;
+w1=1.2;
+w2=0.5;
+h=0.4;
+g=9.81;
+//To Find
+theta=((w1-w2)/2)/h;disp(theta);
+q1=((2/3)*Cd1*sqrt(2*g)*H^(3/2));
+q2=((8/15)*Cd2*sqrt(2*g)*theta*H^(5/2));
+q=q1+q2;disp(q1);disp(q2);
+disp("discharge through Trapezoidal Notch ="+string(q)+" m^3/sec");
diff --git a/2021/CH9/EX9.5/EX9_5.pdf b/2021/CH9/EX9.5/EX9_5.pdf
new file mode 100755
index 000000000..07d50302f
--- /dev/null
+++ b/2021/CH9/EX9.5/EX9_5.pdf
diff --git a/2021/CH9/EX9.5/Ex9_5.sce b/2021/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..3e3596865
--- /dev/null
+++ b/2021/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,13 @@
+//Finding of Percentage Error in Discharge
+//Given
+Cd=0.6;
+q=40000;
+L=0.5;
+H=0.2;
+g=9.81;
+//To Find
+H1=q/((2/3)*Cd*L*sqrt(2*g));
+H2=H1^(2/3);
+H3=H2/100;
+dq=(3/2)*(H/H3)*100;
+disp("Percentage Error in Discharge ="+string(dq)+" Percentage");
diff --git a/2021/CH9/EX9.8/EX9_8.pdf b/2021/CH9/EX9.8/EX9_8.pdf
new file mode 100755
index 000000000..030dc695d
--- /dev/null
+++ b/2021/CH9/EX9.8/EX9_8.pdf
diff --git a/2021/CH9/EX9.8/Ex9_8.sce b/2021/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..e3b64dba2
--- /dev/null
+++ b/2021/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,8 @@
+//Finding of Discharge over a Cipolletti weir
+//Given
+L=1.8;
+H=1.2;
+Cd=0.632;
+//To Find
+q=1.866*L*H^(3/2);
+disp("Discharge over a Cipolletti weir ="+string(q)+" m^3/sec");
diff --git a/2021/CH9/EX9.9/EX9_9.pdf b/2021/CH9/EX9.9/EX9_9.pdf
new file mode 100755
index 000000000..af11c2ad9
--- /dev/null
+++ b/2021/CH9/EX9.9/EX9_9.pdf
diff --git a/2021/CH9/EX9.9/Ex9_9.sce b/2021/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..a6a876e60
--- /dev/null
+++ b/2021/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,13 @@
+//Finding of Dischage
+//Given
+Cd1=0.6;
+Cd2=0.8;
+L=3.5;
+g=9.81;
+H1=0.3;
+H2=0.15;
+//To Find
+q1=((2/3)*Cd1*L*sqrt(2*g)*(H1-H2)^(3/2));
+q2=Cd2*L*H2*sqrt(2*g*(H1-H2));
+q3=q1+q2;
+disp("Dischage ="+string(q3)+" m^3/sec");
author	priyanka	2015-06-24 15:03:17 +0530
committer	priyanka	2015-06-24 15:03:17 +0530
commit	b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
tree	ab291cffc65280e58ac82470ba63fbcca7805165 /2021
download	Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip