initial commit / add all books

23 files changed, 507 insertions, 0 deletions

diff --git a/2132/CH5/EX5.1/Example5_1.sce b/2132/CH5/EX5.1/Example5_1.sce
new file mode 100755
index 000000000..6f54cda80
--- /dev/null
+++ b/2132/CH5/EX5.1/Example5_1.sce
@@ -0,0 +1,15 @@
+//Example 5.1

+clc;

+clear;

+close;

+format('v',9);

+//Given data :

+m=2000;//litre or kg(1litre water =1kg)

+M=m/60;//kg/s

+p=4.5;//bar

+p=p*10^5;//N/m^2

+g=9.81;//constant

+w=g*1000;//N/m^3

+H=p/w;//m

+Power=M*g*H/1000;//kW

+disp(Power,"Power required in kW : ");

diff --git a/2132/CH5/EX5.10/Example5_10.sce b/2132/CH5/EX5.10/Example5_10.sce
new file mode 100755
index 000000000..4c6794a9e
--- /dev/null
+++ b/2132/CH5/EX5.10/Example5_10.sce
@@ -0,0 +1,24 @@
+//Example 5.10

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+Q=200*10^-3;//m^3/s

+D1=250/1000;//meter

+D2=200/1000;//meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1=Q/A1;//m/s

+v2=Q/A2;//m/s

+Z1=2;//meter

+Z2=8;//meter

+g=9.81;//constant

+w=1000;//kg/m^3

+p1=w*(Z1-v1^2/2/g);//kg/m^2

+p2=v1^2/2/g*w+p1+Z2*w-v2^2/2/g*w-4*w;//kg/m^2(by Bernolli's theorem)

+p1=p1*g;//N/m^2

+p2=p2*g;//N/m^2

+disp(p1,"Pressure intensity at point P in N/m^2 : ");

+disp(p2,"Pressure intensity at point Q in N/m^2 : ");        

+//Answer in the book is not accurate.

diff --git a/2132/CH5/EX5.11/Example5_11.sce b/2132/CH5/EX5.11/Example5_11.sce
new file mode 100755
index 000000000..8622c839d
--- /dev/null
+++ b/2132/CH5/EX5.11/Example5_11.sce
@@ -0,0 +1,25 @@
+//Example 5.11

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+slope=1/10;

+Z1=0;//meter

+Z2=40*slope;//meter

+p1=1.5/10^-4;//kg/cm^2

+v2=4.1;//m/s

+D1=600/1000;//meter

+D2=300/1000;//meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1=A2*v2/A1;//m/s

+g=9.81;//constant

+w=1000;//kg/m^3

+p2=(p1/w+v1^2/2/g+Z1-v2^2/2/g-Z2)*w;//kg/m^2(by Bernolli's theorem)

+p2=p2*10^-4;//kg/cm^2

+Q1=A1*v1;//m^3/sec

+Q1=Q1*1000;//litre/sec

+disp(p2,"Pressure intensity at point Q in kg/cm^2 : ");        

+disp(Q1,"Discharge of pipe in litres/sec : ");

+//Answer in the book is not accurate. calculation for A1 & A2 is wrong.

diff --git a/2132/CH5/EX5.12/Example5_12.sce b/2132/CH5/EX5.12/Example5_12.sce
new file mode 100755
index 000000000..9f568efa7
--- /dev/null
+++ b/2132/CH5/EX5.12/Example5_12.sce
@@ -0,0 +1,20 @@
+//Example 5.12

+clc;

+clear;

+close;

+format('v',6);

+//Given data :

+D1=180/1000;//meter

+D2=90/1000;//meter

+g=9.81;//gravity constant

+S=0.8;//sp. gravity of oil

+Sm=13.6;//sp. gravity of mercury

+x=300/1000;//meter

+K=0.97;//coeff. of meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=A1*A2*sqrt(2*g)/sqrt(A1^2-A2^2)

+h=x*(Sm/S-1);//meter of oil

+Q=K*C*sqrt(h);//m^3/sec

+Q=Q*1000;//litre/sec

+disp(Q,"Discharge of oil in litres/sec : ");

diff --git a/2132/CH5/EX5.13/Example5_13.sce b/2132/CH5/EX5.13/Example5_13.sce
new file mode 100755
index 000000000..8ac1f8e46
--- /dev/null
+++ b/2132/CH5/EX5.13/Example5_13.sce
@@ -0,0 +1,18 @@
+//Example 5.13

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+D1byD2=1/0.7;

+D1=320/1000;//meter

+D2=320*0.7/1000;//meter

+g=9.81;//gravity constant

+Q=30.6/60;//m^3/sec

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=A1*sqrt(2*g)/sqrt((D1byD2)^4-1);

+h=1.2;//meter of water

+K=Q/C/sqrt(h);//Coeff. of meter

+disp(K,"Coefficient of meter : ");

+//Answer in the book is wrong.

diff --git a/2132/CH5/EX5.14/Example5_14.sce b/2132/CH5/EX5.14/Example5_14.sce
new file mode 100755
index 000000000..ea38e24a0
--- /dev/null
+++ b/2132/CH5/EX5.14/Example5_14.sce
@@ -0,0 +1,19 @@
+//Example 5.14

+clc;

+clear;

+close;

+format('v',5);

+//Given data :

+D1=320/1000;//meter

+D2=224/1000;//meter

+g=9.81;//gravity constant

+Q=25000/1000/60;//m^3/sec

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=0.4984;//venturi constant

+K=0.92;//Coeff. of meter

+h=(Q/K/C)^2

+S=1;//sp. gravity

+Sm=13.6;//sp. gravity

+x=h/(Sm/S-1);//meter of water

+disp(x*1000,"Deflection in  manometer(mm) : ");

diff --git a/2132/CH5/EX5.15/Example5_15.sce b/2132/CH5/EX5.15/Example5_15.sce
new file mode 100755
index 000000000..53068d13e
--- /dev/null
+++ b/2132/CH5/EX5.15/Example5_15.sce
@@ -0,0 +1,23 @@
+//Example 5.15

+clc;

+clear;

+close;

+format('v',6);

+//Given data :

+D1=120/1000;//meter

+D2=120*0.55/1000;//meter

+g=9.81;//gravity constant

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+Q=30/1000;//m^3/sec

+C=A1*sqrt(2*g)/sqrt((D1/D2)^4-1);//venturi constant

+K=0.94;//Coeff. of meter

+h=(Q/K/C)^2;//meter

+Z1=0;//meter

+Z2=0.3;//meter

+S=0.79;//sp. gravity

+w=1000*S;//kg/m^3

+delta_p=(h+Z1-Z2)*w;//kg/m^2

+delta_p=delta_p*g;//N/m^2

+disp(delta_p,"Pressure difference in N/m^2 : ");

+//answer is wrong in the book.

diff --git a/2132/CH5/EX5.16/Example5_16.sce b/2132/CH5/EX5.16/Example5_16.sce
new file mode 100755
index 000000000..fdca4d8c8
--- /dev/null
+++ b/2132/CH5/EX5.16/Example5_16.sce
@@ -0,0 +1,19 @@
+//Example 5.16

+clc;

+clear;

+close;

+format('v',6);

+//Given data :

+D1=160/1000;//meter

+D2=60/1000;//meter

+g=9.81;//gravity constant

+S=0.8;//sp. gravity

+Sm=13.6;//sp. gravity of mercury

+Q=0.05;//m^3/sec

+K=0.98;//Coeff. of meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=A1*sqrt(2*g)/sqrt((A1/A2)^2-1);//venturi constant

+h=(Q/K/C)^2;//meter

+x=h/(Sm/S-1);//meter

+disp(x,"Deflection in meter : ");

diff --git a/2132/CH5/EX5.17/Example5_17.sce b/2132/CH5/EX5.17/Example5_17.sce
new file mode 100755
index 000000000..295f7f743
--- /dev/null
+++ b/2132/CH5/EX5.17/Example5_17.sce
@@ -0,0 +1,20 @@
+//Example 5.17

+clc;

+clear;

+close;

+format('v',4);

+//Given data :

+D1=200/1000;//meter

+D2=100/1000;//meter

+x=220/1000;//meter

+g=9.81;//gravity constant

+K=0.98;//Coeff. of meter

+S=1;//sp. gravity

+Sm=13.6;//sp. gravity of mercury

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=A1*sqrt(2*g)/sqrt((A1/A2)^2-1);//venturi constant

+h=x*(Sm/S-1);//meter

+Q=K*C*sqrt(h);//m^3/sec

+Q=Q*1000;//litres/sec

+disp(Q,"Rate of flow in litres/sec : ");

diff --git a/2132/CH5/EX5.18/Example5_18.sce b/2132/CH5/EX5.18/Example5_18.sce
new file mode 100755
index 000000000..cc64aecb0
--- /dev/null
+++ b/2132/CH5/EX5.18/Example5_18.sce
@@ -0,0 +1,21 @@
+//Example 5.18

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+D1=40/100;//meter

+D2=15/100;//meter

+x=25/100;//meter

+g=9.81;//gravity constant

+K=0.98;//Coeff. of meter

+S=1;//sp. gravity

+Sm=13.6;//sp. gravity of mercury

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+C=A1*A2*sqrt(2*g)/sqrt(A1^2-A2^2);//venturi constant

+h=x*(Sm/S-1);//meter

+Q=K*C*sqrt(h);//m^3/sec

+Q=Q*1000*3600;//litres/hour

+disp(Q,"Flow of water in litres/hour : ");

+//Answer in the book is wrong.

diff --git a/2132/CH5/EX5.19/Example5_19.sce b/2132/CH5/EX5.19/Example5_19.sce
new file mode 100755
index 000000000..f03adb9b0
--- /dev/null
+++ b/2132/CH5/EX5.19/Example5_19.sce
@@ -0,0 +1,22 @@
+//Example 5.19

+clc;

+clear;

+close;

+format('v',4);

+//Given data :

+D1=15/100;//meter

+D2=7.5/100;//meter

+g=9.81;//gravity constant

+p1=4*g*10^4;//N/m^2

+p2=1.5*g*10^4;//kg/cm^2

+w=9.81;//kg/m^2

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1BYv2=A2/A1;

+//v1^2/2/g+p1/w=v2^2/2/g+p2/w

+//v1^2=v2^2-50*g

+v2=sqrt(50*g/(1-v1BYv2^2));//m/s

+Q=A2*v2;//m^3/sec

+Q=Q*1000;//litres/sec

+disp(Q,"Flow of water in litres/sec : ");

+//Answer is wrong in the book.

diff --git a/2132/CH5/EX5.2/Example5_2.sce b/2132/CH5/EX5.2/Example5_2.sce
new file mode 100755
index 000000000..9f2a095d5
--- /dev/null
+++ b/2132/CH5/EX5.2/Example5_2.sce
@@ -0,0 +1,16 @@
+//Example 5.2

+clc;

+clear;

+close;

+format('v',9);

+//Given data :

+v1=400*10^-3;//m/s

+d1=300/1000;//meter

+d2=450/1000;//meter

+A1=%pi*d1^2/4;//m^2

+A2=%pi*d2^2/4;//m^2

+Q1=A1*v1*1000;//litres/sec(1m^3=1000litres)

+disp(Q1,"Discharge of pipe in litres/sec : ");

+v2=(Q1/1000)/A2;//m/s(Q1=Q2)

+disp(v2,"Mean velocity of flow in m/s : ");

+//Answer of discharge is wrong in the book.

diff --git a/2132/CH5/EX5.20/Example5_20.sce b/2132/CH5/EX5.20/Example5_20.sce
new file mode 100755
index 000000000..08905eba9
--- /dev/null
+++ b/2132/CH5/EX5.20/Example5_20.sce
@@ -0,0 +1,17 @@
+//Example 5.20

+clc;

+clear;

+close;

+format('v',6);

+//Given data :

+D1=20/100;//meter

+D2=15/100;//meter

+A1=%pi/4*D1^2;//m^2

+A2=%pi/4*D2^2;//m^2

+v1=2;//m/s

+v2=A1*v1/A2;//m/s

+disp(v2,"Velocity at another section in m/s : ");

+FlowRate=A1*v1;//m^3/s

+FlowRate=FlowRate*1000;//litres/s

+disp(FlowRate,"Flow Rate in litres/sec : ");

+//Answer of velocity in the book is not accurate.

diff --git a/2132/CH5/EX5.21/Example5_21.sce b/2132/CH5/EX5.21/Example5_21.sce
new file mode 100755
index 000000000..c072f4065
--- /dev/null
+++ b/2132/CH5/EX5.21/Example5_21.sce
@@ -0,0 +1,21 @@
+//Example 5.21

+clc;

+clear;

+close;

+format('v',4);

+//Given data :

+rd=0.75;//relative density

+D=12.5/100;//meter

+p=1;//bar

+p=p*1.02;//kg/cm^2

+p=p*9.81*10^4/1000;//kPa

+g=9.81;//gravity constant

+w=g*rd;//N/m^3

+pH=p/w;//meter

+Z=2.5;//meter

+Et=20;//Nm

+v=sqrt((Et-p/w-Z)*2*g);//m/s

+A=%pi/4*D^2;//m^2

+Q=A*v;//m^3/sec

+Q=A*v*1000;//litres/sec

+disp(Q,"Flow Rate of oil in litres/sec : ");

diff --git a/2132/CH5/EX5.22/Example5_22.sce b/2132/CH5/EX5.22/Example5_22.sce
new file mode 100755
index 000000000..3b6206f3e
--- /dev/null
+++ b/2132/CH5/EX5.22/Example5_22.sce
@@ -0,0 +1,21 @@
+//Example 5.22

+clc;

+clear;

+close;

+format('v',5);

+//Given data :

+rd=0.75;//relative density

+d1=0.3;//meter

+d2=0.1;//meter

+Q=50/1000;//m^3/sec

+A1=%pi/4*d1^2;//m^2

+A2=%pi/4*d2^2;//m^2

+v1=Q/A1;//m/s

+v2=A1*v1/A2;//m/s

+p1=200;//kN/m^2

+p2=100;//kN/m^2

+w=9.81;//kN/m^3

+g=9.81;//gravity constant

+Z1=0;//meter

+Z2=Z1+p1/w+v1^2/2/g-p2/w-v2^2/2/g;//meter

+disp(Z2,"Z in meter : ");

diff --git a/2132/CH5/EX5.23/Example5_23.sce b/2132/CH5/EX5.23/Example5_23.sce
new file mode 100755
index 000000000..7231682df
--- /dev/null
+++ b/2132/CH5/EX5.23/Example5_23.sce
@@ -0,0 +1,22 @@
+//Example 5.23

+clc;

+clear;

+close;

+format('v',6);

+//Given data :

+D1=300/1000;//meter

+D2=150/1000;//meter

+Q=50/1000;//m^3/sec

+A1=%pi/4*D1^2;//m^2

+A2=%pi/4*D2^2;//m^2

+delpBYw=3;//p1/w-p2/w=3;//m

+v1BYv2=A2/A1;

+Z1=0;//meter

+Z2=0;//meter

+g=9.81;//gravity constant

+//HeadLoss=1/8*v^2/2/g

+//Z1+p1/w+v1^2/2/g=Z2+p2/w+v2^2/2/g+HeadLoss

+v2=sqrt((Z1-Z2+delpBYw)/(1/2/g-v1BYv2^2/2/g+1/8/2/g));//m/s

+Q=A2*v2;//m^3/s

+Q=Q*1000;//litres/sec

+disp(Q,"Discharge in pipe in litres/sec : ");

diff --git a/2132/CH5/EX5.3/Example5_3.sce b/2132/CH5/EX5.3/Example5_3.sce
new file mode 100755
index 000000000..78755320f
--- /dev/null
+++ b/2132/CH5/EX5.3/Example5_3.sce
@@ -0,0 +1,19 @@
+//Example 5.3

+clc;

+clear;

+close;

+format('v',7);

+//Given data :

+PotentialHead=2;//meter of fluid

+disp("Potential Head is "+string(PotentialHead)+" meter of fluid.");

+v=5;//m/s

+g=9.81;//constant

+VelocityHead=v^2/2/g;//m

+disp("Velocity Head is "+string(VelocityHead)+" meter of fluid.");

+w=g*1000;//N/m^3

+S=0.8;//sp. gravity of fluid

+p=200;//kPa

+PressureHead=p*10^3/w/S;//meter of fluid

+disp("Pressure Head is "+string(PressureHead)+" meter of fluid.");

+TotalHead=PotentialHead+VelocityHead+PressureHead;//meter of fluid

+disp("Total Head is "+string(TotalHead)+" meter of fluid.");

diff --git a/2132/CH5/EX5.4/Example5_4.sce b/2132/CH5/EX5.4/Example5_4.sce
new file mode 100755
index 000000000..d51f20fe6
--- /dev/null
+++ b/2132/CH5/EX5.4/Example5_4.sce
@@ -0,0 +1,15 @@
+//Example 5.4

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+p=0.8/10^-4;//kg/m^2

+datumH=4;//meter

+v=0.8;//m/s

+g=9.81;//constant

+VelocityH=v^2/2/g;//m

+w=1000;//kg/m^3

+PressureH=p/w;//meter of fluid

+TotalH=datumH+VelocityH+PressureH;//meter of fluid

+disp("Total Energy is "+string(TotalH)+" meter.");

diff --git a/2132/CH5/EX5.5/Example5_5.sce b/2132/CH5/EX5.5/Example5_5.sce
new file mode 100755
index 000000000..d11fbc687
--- /dev/null
+++ b/2132/CH5/EX5.5/Example5_5.sce
@@ -0,0 +1,29 @@
+//Example 5.5

+clc;

+clear;

+close;

+format('v',7);

+//Given data :

+D1=800/1000;//m^2

+D2=600/1000;//m^2

+p1=100;//kPa

+p2=40;//kPa

+v1=4000*10^-3;//m/s

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+Z1=4;//meter

+Z2=7;//meter

+rho=1;//sp. gravity

+g=9.81;//constant

+PHeadA=p1/rho/g;//meter of fluid

+PHeadB=p2/rho/g;//meter of fluid

+v2=A1*v1/A2;//m/s

+VHeadA=v1^2/2/g;//meter

+VHeadB=v2^2/2/g;//meter

+E1=Z1+PHeadA+VHeadA;//meter

+E2=Z2+PHeadB+VHeadB;//meter

+if E1>E2 then

+    disp("Total Energy at A("+string(E1)+" meter) is greater than total energy at B("+string(E2)+" meter). Flow of water is from A to B.");

+else

+        disp("Total Energy at B("+string(E2)+" meter) is greater than total energy at A("+string(E1)+" meter). Flow of water is from B to A.");

+end

diff --git a/2132/CH5/EX5.6/Example5_6.sce b/2132/CH5/EX5.6/Example5_6.sce
new file mode 100755
index 000000000..45144cb25
--- /dev/null
+++ b/2132/CH5/EX5.6/Example5_6.sce
@@ -0,0 +1,29 @@
+//Example 5.6

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+D1=1.25;//meter

+D2=0.625;//meter

+slope=100;

+L=300;///meter

+g=9.81;//constant

+Z12=L/slope;//meter

+Q=100;//litres/sec

+Q=Q*10^-3;//m^3/sec

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1=Q/A1;//m/s

+v2=Q/A2;//m/s

+p1=100;//kN/m^2

+//Higher End : 

+w=9.81;//kN/m^3

+Phead=p1/w;//meter

+Vhead=v1^2/2/g;//meter

+//Lower End : 

+w=9.81;//kN/m^3

+//Phead=p1/w;//meter

+Vhead=v2^2/2/g;//meter

+p2=(Z12+v1^2/2/g+p1/w-v2^2/2/g)*w;//kN/m^2(By Bernoulli's theorem)

+disp(p2,"Pressure at the lower end in kN per m^2 : ");

diff --git a/2132/CH5/EX5.7/Example5_7.sce b/2132/CH5/EX5.7/Example5_7.sce
new file mode 100755
index 000000000..5682da970
--- /dev/null
+++ b/2132/CH5/EX5.7/Example5_7.sce
@@ -0,0 +1,34 @@
+//Example 5.7

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+Z1=0;//meter

+Z2=5;//meter

+Q=300*10^-3;//m/s

+D1=0.3;//meter

+D2=0.6;//meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1=Q/A1;//m/s

+v2=Q/A2;//m/s

+p1=100;//kN/m^2

+p2=600;//kN/m^2

+g=9.81;//constant

+Vhead11=v1^2/2/g;//meter

+Vhead22=v2^2/2/g;//meter

+Phead11=p1/g;//meter

+Phead22=p2/g;//meter

+E1_11=Z1+Vhead11+Phead11;//meter

+E2_22=Z2+Vhead22+Phead22;//meter

+if E1_11>E2_22 then

+    disp("Energy at section 1-1("+string(E1_11)+" meter) is greater than energy at section 2-2("+string(E2_22)+" meter). Flow of water is from section 1-1 to 2-2.");

+    HeadLoss=E1_11-E2_22;//meter

+    disp(HeadLoss,"Head Loss in meter : ");

+else

+        disp("Energy at section 2-2("+string(E2_22)+" meter) is greater than energy at section 1-1("+string(E1_11)+" meter). Flow of water is from section 2-2 to 1-1.");  

+            HeadLoss=E2_22-E1_11;//meter

+        disp(HeadLoss,"Head Loss in meter : ");

+end

+

diff --git a/2132/CH5/EX5.8/Example5_8.sce b/2132/CH5/EX5.8/Example5_8.sce
new file mode 100755
index 000000000..fa7cdc082
--- /dev/null
+++ b/2132/CH5/EX5.8/Example5_8.sce
@@ -0,0 +1,23 @@
+//Example 5.8

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+D=400/1000;//meter

+v1=20;//m/s

+Z1=28;//meter

+Z2=31;//meter

+p1=4/10^-4;//kg/m^2

+p2=3/10^-4;//kg/m^2

+g=9.81;//constant

+w=1000;//kg/m^3

+Vhead1=v1^2/2/g;//meter

+Phead1=p1/w;//meter

+Vhead2=Vhead1;//meter

+Phead2=p2/w;//meter

+E1=Z1+Vhead1+Phead1;//meter

+E2=Z2+Vhead2+Phead2;//meter

+HL=E1-E2;//meter

+disp(HL,"Loss of head between P & Q in meter : ");

+        

diff --git a/2132/CH5/EX5.9/Example5_9.sce b/2132/CH5/EX5.9/Example5_9.sce
new file mode 100755
index 000000000..82595d8b8
--- /dev/null
+++ b/2132/CH5/EX5.9/Example5_9.sce
@@ -0,0 +1,35 @@
+//Example 5.9

+clc;

+clear;

+close;

+format('v',8);

+//Given data :

+Z1=0;//meter

+Z2=4;//meter

+rho=0.8;//sp. gravity

+

+Q=250*10^-3;//m/s or cumec

+D1=250/1000;//meter

+D2=500/1000;//meter

+A1=%pi*D1^2/4;//m^2

+A2=%pi*D2^2/4;//m^2

+v1=Q/A1;//m/s

+v2=Q/A2;//m/s

+p1=0.1*10^3;//N/m^2

+p2=0.06*10^3;//N/m^2

+g=9.81;//constant

+Vhead1=v1^2/2/g;//meter

+Phead1=p1/rho/g;//meter

+Vhead2=v2^2/2/g;//meter

+Phead2=p2/rho/g;//meter

+H1=Z1+Vhead1+Phead1;//meter

+H2=Z2+Vhead2+Phead2;//meter

+if H1>H2 then

+    disp("Total head at A("+string(H1)+" meter) is greater than total head at B("+string(H2)+" meter). Flow will take place from A-B.");

+    HeadLoss=H1-H2;//meter

+    disp(HeadLoss,"Head Loss in meter : ");

+else

+    disp("Total head at B("+string(H2)+" meter) is greater than total head at A("+string(H1)+" meter). Flow will take place from B-A.");

+    HeadLoss=H2-H1;//meter

+    disp(HeadLoss,"Head Loss in meter : ");end

+