initial commit / add all books

3 files changed, 58 insertions, 0 deletions

diff --git a/1055/CH19/EX19.1/ch19_1.sce b/1055/CH19/EX19.1/ch19_1.sce
new file mode 100755
index 000000000..ebcea34d7
--- /dev/null
+++ b/1055/CH19/EX19.1/ch19_1.sce
@@ -0,0 +1,24 @@
+// To Determine the economic operating schedule and the corresponding cost of generation.(b)Determine the savings obtained by loading the units.

+clear 

+clc;

+//dF1/dP1=.4*P1+40 per MWhr

+//dF2/dP2=.5*P1+30 per MWhr

+mprintf("two equations are :\n");

+mprintf("%.1f P1 %.1f P2 = %.1f\n",.4,-.5,-10);

+mprintf("%.1f P1+ %.1fP2 = %.1f\n",1,1,180);

+A=[.4 -.5;1 1];

+B=[-10;180];

+P=(inv(A))*B;

+P1=P(1,1);

+P2=P(2,1);

+F1=.2*(P1)^2 +40*P1+120;

+F2=.25*(P2)^2+30*P2+150;

+Total=F1+F2;//Total cost

+mprintf("(a)Cost of Generation=Rs %.2f /hr\n",Total);

+P1=90;

+P2=90;

+F1=.2*(P1)^2 +40*P1+120;

+F2=.25*(P2)^2+30*P2+150;

+Total2=F1+F2;//Total cost

+savings=Total2-Total

+mprintf("(b)Savings=Rs %.2f /hr\n",savings)

diff --git a/1055/CH19/EX19.2/ch19_2.sce b/1055/CH19/EX19.2/ch19_2.sce
new file mode 100755
index 000000000..f456cc345
--- /dev/null
+++ b/1055/CH19/EX19.2/ch19_2.sce
@@ -0,0 +1,7 @@
+//Determine the incremental cost of recieved power and penalty factor of the plant

+clear 

+clc;

+pf=10/8;//penalty factor

+cost=(.1*10+3)*pf;//Cost of recieved power=dF1/dP1

+mprintf("Penalty Factor=%.1f\n",pf);

+mprintf("Cost of recieved Power=Rs %.1f /MWhr",cost);

diff --git a/1055/CH19/EX19.4/ch19_4.sce b/1055/CH19/EX19.4/ch19_4.sce
new file mode 100755
index 000000000..1a9734c65
--- /dev/null
+++ b/1055/CH19/EX19.4/ch19_4.sce
@@ -0,0 +1,27 @@
+//Determine the minimum cost of generation .

+clear 

+clc;

+//dF1/dP1=.048*P1+8

+//dF2/dP2=.08*P1+6

+mprintf("two equations are :\n");

+mprintf("%.3f P1 %.2f P2 = %.1f\n",.048,-.08,-2);

+mprintf("%.1f P1+ %.1fP2 = %.1f\n",1,1,50);

+A=[.048 -.08;1 1];

+B=[-2;50];

+P=(inv(A))*B;

+P1=P(1,1);

+P2=P(2,1);

+F1=(.024*(P1)^2 +8*P1+80)*(10^6);

+F2=(.04*(P2)^2+6*P2+120)*(10^6);

+mprintf("when load is 150MW , equations are: :\n");

+mprintf("%.3f P1 %.2f P2 = %.1f\n",.048,-.08,-2);

+mprintf("%.1f P1+ %.1fP2 = %.1f\n",1,1,150);

+A=[.048 -.08;1 1];

+B=[-2;150];

+P=(inv(A))*B;

+P1=P(1,1);

+P2=P(2,1);

+f1=(.024*(P1)^2 +8*P1+80)*(10^6);

+f2=(.04*(P2)^2+6*P2+120)*(10^6);

+Total=(F1+F2+f1+f2)*12*2/(10^6);

+mprintf("Total cost=Rs. %.2f",Total)