7 files changed, 204 insertions, 0 deletions
diff --git a/416/CH11/EX11.1/exp11_1.png b/416/CH11/EX11.1/exp11_1.png
new file mode 100755
index 000000000..9aa8ddfd6
--- /dev/null
+++ b/416/CH11/EX11.1/exp11_1.png
diff --git a/416/CH11/EX11.1/exp11_1pp.sce b/416/CH11/EX11.1/exp11_1pp.sce
new file mode 100755
index 000000000..ca1f515fa
--- /dev/null
+++ b/416/CH11/EX11.1/exp11_1pp.sce
@@ -0,0 +1,95 @@
+clc
+clear
+disp("example 11 1")
+wd=[0 5 8 12 13 17 21 24] //given week days
+wlld=[100 150 250 100 250 350 150] //given load in week days
+wld=[wlld 0]
+we=[0 5 17 21 24]//given week ends
+wed=[100 150 200 150]//given load in week ends
+wed=[wed 0]
+h=90//head
+f=50 //flow
+et=0.97//is available for 97 persent
+eff=0.9 //efficiency
+tl=0.05 //transmission loss
+pa=735.5*f*h*eff/75 //power available
+nap=pa*(1-tl) //net available power
+he=nap*24/1000 //hydro energy for 24 in MW
+he1=round(he/100)*100
+[m,n]=size(wd)
+[x,y]=min(wlld)
+[q,r]=max(wlld)
+for i=1:n-1
+    fl(i)=wd(i+1)-wd(i)
+end
+[o,p]=size(we)
+for i=1:p-1
+    fll(i)=we(i+1)-we(i)
+end
+for j=x:10:q
+    pp=wlld-j
+    for l=1:n-1
+    if pp(l)<0 then
+        pp(l)=0
+    end
+    end
+    heq=pp*fl
+    heq=round(heq/100)*100
+    if heq==he1 then
+        break
+    end
+end//rearrangeing for plot
+subplot(211) 
+plot2d2(wd,wld)
+xtitle("chronological load curve for week day for example 11.1","hour of day","load MW")
+subplot(212)
+plot2d2(we,wed)
+xtitle("chronological load curve for weekend day for example 11.1","hour of day","load MW")
+
+printf("power available from the hydro plant for %dMW of the time is %.2fMW",et*100,pa/1000)
+printf("\nnet available hydra power after taking transmission loss into account %.2fMW",nap/100)
+printf("\nhydro energy available during 24 hours %.2fMW",he)
+printf("\nthe magnitude of hydro power is %dMW \ntotal capacity of hydro plant required on week days %dMW ",q-j,(q-j)/(1-tl))
+printf("capacity of thermal plant on week days %dMW",q)
+printf("\nthe schedule for hydro plant is on week days")
+for i=1:n
+    if wd(i)>12 then
+        wd(i)=wd(i)-12
+    end
+end
+disp(wd)
+disp(round(pp/(1-tl)))
+disp("the schedule for thermal plant is on week days")
+disp(wd)
+disp(wlld-pp)
+[m,n]=size(we)
+[x,y]=min(wed)
+[q,r]=max(wed)
+for j=x:10:q
+    pp=wed-j
+    for l=1:n-1
+    if pp(l)<0 then
+        pp(l)=0
+    end
+end
+pp(n)=[]
+    heq=pp*fll
+    heq=floor(heq/100)*100
+    if heq==he1 then
+        break
+    end
+end
+printf("\nthe magnitude of hydro power is %dMW \ntotal capacity of hydro plant required on week ends %dMW ",q-j,(q-j)/(1-tl))
+printf("capacity of thermal plant on week ends %dMW",q)
+printf("\nthe schedule for hydro plant is on week ends")
+for i=1:n
+    if we(i)>12 then
+        we(i)=we(i)-12
+    end
+end
+disp(we)
+disp(round(pp/(1-tl)))
+disp("the schedule for thermal plant is on week days")
+disp(we)
+pp(n)=0
+disp(wed-pp)
diff --git a/416/CH11/EX11.1/exp11_1pp.txt b/416/CH11/EX11.1/exp11_1pp.txt
new file mode 100755
index 000000000..e9344f81f
--- /dev/null
+++ b/416/CH11/EX11.1/exp11_1pp.txt
@@ -0,0 +1,31 @@
+  
+ example 11 1   
+power available from the hydro plant for 97MW of the time is 39.72MW
+net available hydra power after taking transmission loss into account 377.31MW
+hydro energy available during 24 hours 905.55MW
+the magnitude of hydro power is 140MW 
+total capacity of hydro plant required on week days 147MW capacity of thermal plant on week days 350MW
+the schedule for hydro plant is on week days 
+    0.    5.    8.    12.    1.    5.    9.    12.  
+ 
+    0.    0.    42.    0.    42.    147.    0.  
+ 
+ the schedule for thermal plant is on week days   
+ 
+    0.    5.    8.    12.    1.    5.    9.    12.  
+ 
+    100.    150.    210.    100.    210.    210.    150.  
+
+the magnitude of hydro power is 90MW 
+total capacity of hydro plant required on week ends 94MW capacity of thermal plant on week ends 200MW
+the schedule for hydro plant is on week ends 
+    0.    5.    5.    9.    12.  
+ 
+    0.    42.    95.    42.  
+ 
+ the schedule for thermal plant is on week days   
+ 
+    0.    5.    5.    9.    12.  
+ 
+    100.    110.    110.    110.    0.  
+ 
+\ No newline at end of file
diff --git a/416/CH11/EX11.2/exp11_2pp.sce b/416/CH11/EX11.2/exp11_2pp.sce
new file mode 100755
index 000000000..1bb5de4f7
--- /dev/null
+++ b/416/CH11/EX11.2/exp11_2pp.sce
@@ -0,0 +1,38 @@
+clc
+clear
+disp("example 11.2")
+//given
+l1=700;t1=14;l2=500;t2=10
+ac=24;bc=0.02//variables of cost equation
+aw=6;bw=0.0025 //variables of watere quantity equation
+b22=0.0005 //loss coefficient
+r2=2.5
+lam=1:0.001:40
+gg=1;q=1
+for lam=25:0.001:40
+    a=[2*bc 0;0 r2*bw*2+2*b22*lam]
+    b=[lam-ac;lam-aw*r2]
+    p=inv(a)*b
+    g=round(p(1)+p(2))
+    l=round(l1+b22*p(2)^2)
+    lq=round(l2+b22*p(2)^2)
+    if g>=l then
+        printf("\nfor load condition %dMW \n then, \n \t lamda %f \t p1=%fMW \n \t p2=%fMW \t pl=%fMW",l1,lam,p(1),p(2),2*b22*p(2))
+        break
+    end
+end
+for lam=25:0.001:40
+    a=[2*bc 0;0 r2*bw*2+2*b22*lam]
+    b=[lam-ac;lam-aw*r2]
+    pq=inv(a)*b
+    g=round(pq(1)+pq(2))
+    lq=round(l2+b22*pq(2)^2)
+
+    if g>=lq then
+        printf("\nfor load condition %dMW \n then, \n \t lamda %f \t p1=%fMW \n \t p2=%fMW \t pl=%fMW",l2,lam,pq(1),pq(2),2*b22*pq(2))
+        break
+    end
+end
+dwu=[(aw+bw*p(2))*p(2)*t1+t2*(aw+bw*pq(2))*pq(2)]*3600
+doc=[(ac+bc*p(1))*p(1)*t1+(ac+bc*pq(1))*pq(1)*t2]
+printf("\ndaily water used %fm^3 \ndaily operating cost of thermal plant Rs%f",dwu,doc)
diff --git a/416/CH11/EX11.2/exp11_2pp.txt b/416/CH11/EX11.2/exp11_2pp.txt
new file mode 100755
index 000000000..342f0be48
--- /dev/null
+++ b/416/CH11/EX11.2/exp11_2pp.txt
@@ -0,0 +1,13 @@
+ 
+ example 11.2   
+
+for load condition 700MW 
+ then, 
+ 	 lamda 37.918000 	 p1=347.950000MW 
+ 	 p2=454.559879MW 	 pl=0.454560MW
+for load condition 500MW 
+ then, 
+ 	 lamda 31.705000 	 p1=192.625000MW 
+ 	 p2=377.898428MW 	 pl=0.377898MW
+daily water used 257972328.067993m^3 
+daily operating cost of thermal plant Rs204461.454825 
+\ No newline at end of file
diff --git a/416/CH11/EX11.3/exp11_3.sce b/416/CH11/EX11.3/exp11_3.sce
new file mode 100755
index 000000000..c01c08849
--- /dev/null
+++ b/416/CH11/EX11.3/exp11_3.sce
@@ -0,0 +1,23 @@
+clc
+clear
+disp("example 11.3")
+//given 
+p=250//load
+rt=14 //run time
+t=24//total time
+ac=5;bc=8;cc=0.05 //variables of cost equation
+bw=30;cw=0.05 //variables of water per power
+qw=500//quantity of water
+lam=bc+cc*2*p //lambda
+a=-qw*(10^6)/(3600*rt)
+inn=sqrt(bw^2-4*cw*a)
+phh1=(-bw+inn)/(2*cw)//solution of quadratic equation
+phh2=(-bw-inn)/(2*cw)
+if phh1>0 then
+     r=lam/(bw+cw*phh1) 
+    printf(" hydro plant power is %fMW \n the cost of water is %fRs.per hour/m^3/sec",phh1,r) 
+end
+if phh2>0 then
+     r=lam/(bw+cw*phh2) 
+    printf(" hydro plant power is %fMW \n the cost of water is %fRs.per hour/m^3/sec",phh2,r) 
+end
+\ No newline at end of file
diff --git a/416/CH11/EX11.3/exp11_3.txt b/416/CH11/EX11.3/exp11_3.txt
new file mode 100755
index 000000000..062eb7061
--- /dev/null
+++ b/416/CH11/EX11.3/exp11_3.txt
@@ -0,0 +1,4 @@
+
+ example 11.3   
+ hydro plant power is 237.040686MW 
+ the cost of water is 0.788492Rs.per hour/m^3/sec 
+\ No newline at end of file