1 files changed, 32 insertions, 0 deletions
diff --git a/416/CH10/EX10.15/exp10_15pp.sce b/416/CH10/EX10.15/exp10_15pp.sce
new file mode 100755
index 000000000..add8435c7
--- /dev/null
+++ b/416/CH10/EX10.15/exp10_15pp.sce
@@ -0,0 +1,32 @@
+clc
+clear
+disp("example 10_15")
+c1=0.002;b1=0.86;a1=20
+c2=0.004;b2=1.08;a2=20
+c3=0.0028;b3=0.64;a3=36
+fc=500
+maxl=120
+minl=36
+tl=200
+d=[1 1 1;2*fc*c1 -fc*2*c2 0;0 -fc*2*c2 fc*2*c3]
+p=[tl;fc*(b2-b1);fc*(b2-b3)] 
+pp=inv(d)*p //matrix inversion method
+printf("\nloads on generaating station by economic creatirian method is %dMW,%dMW,%dMW",pp(1),pp(2),pp(3))
+for i=1:3
+    if pp(i)<minl then
+        pp(i)=minl
+        printf("\nload on generating station %d is less then minimum value %dMW \n so it is made equal to minimum value %dMW",i,minl,minl)
+        e=[1 1;d(2,1) -d(3,3)]
+        q=[(tl-pp(i));-p(i)]
+        qq=inv(e)*q //matrix inversion method
+printf("\nloads on generaating station by economic creatirian method is %.3fMW,%.3fMW",qq(1),qq(2))
+    end
+    if pp(i)>maxl then
+        pp(i)=maxl
+        printf("\nload on generating station %d is greater than maximum value %dMW \n so it is made equal to mmaximum value %dMW",i,maxl,maxl)
+        e=[1 1;d(2,1) -d(3,3)]
+        q=[(tl-pp(i));-p(i)]
+        qq=inv(e)*q //matrix inversion method
+printf("\nloads on generaating station by economic creatirian method is %.2fMW,%.2fMW",qq(1),qq(2))
+    end
+end