summaryrefslogtreecommitdiff
path: root/2441/CH1
diff options
context:
space:
mode:
Diffstat (limited to '2441/CH1')
-rwxr-xr-x2441/CH1/EX1.1/Ex1_1.sce28
-rwxr-xr-x2441/CH1/EX1.10/Ex1_10.sce12
-rwxr-xr-x2441/CH1/EX1.11/Ex1_11.sce21
-rwxr-xr-x2441/CH1/EX1.12/Ex1_12.sce28
-rwxr-xr-x2441/CH1/EX1.13/Ex1_13.sce8
-rwxr-xr-x2441/CH1/EX1.14/Ex1_14.sce30
-rwxr-xr-x2441/CH1/EX1.16/Ex1_16.sce16
-rwxr-xr-x2441/CH1/EX1.17/Ex1_17.sce9
-rwxr-xr-x2441/CH1/EX1.18/Ex1_18.sce8
-rwxr-xr-x2441/CH1/EX1.19/Ex1_19.sce23
-rwxr-xr-x2441/CH1/EX1.2/Ex1_2.sce46
-rwxr-xr-x2441/CH1/EX1.3/Ex1_3.sce46
-rwxr-xr-x2441/CH1/EX1.4/Ex1_4.sce23
-rwxr-xr-x2441/CH1/EX1.5/Ex1_5.sce27
-rwxr-xr-x2441/CH1/EX1.6/Ex1_6.sce47
-rwxr-xr-x2441/CH1/EX1.7/Ex1_7.sce14
-rwxr-xr-x2441/CH1/EX1.8/Ex1_8.sce28
-rwxr-xr-x2441/CH1/EX1.9/Ex1_9.sce51
18 files changed, 465 insertions, 0 deletions
diff --git a/2441/CH1/EX1.1/Ex1_1.sce b/2441/CH1/EX1.1/Ex1_1.sce
new file mode 100755
index 000000000..5f2292d9c
--- /dev/null
+++ b/2441/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,28 @@
+//exa 1.1
+clc;clear;close;
+format('v',6);
+B=100;//W(8Bulb)
+F=60;//W(2Fan)
+L=100;//W(2Light)
+LoadConnected=8*B+2*F+2*L;//W
+disp(LoadConnected,"(a) Connected Load (W)")
+//12 midnight to 5am
+demand1=1*F;//W
+//5am to 7am
+demand2=2*F+1*L;//W
+//7am to 9am
+demand3=0;//W
+//9am to 6pm
+demand4=2*F;//W
+//6pm to midnight
+demand5=2*F+4*B;//W
+DEMAND=[demand1 demand2 demand3 demand4 demand5]
+max_demand=max(DEMAND);
+disp(max_demand,"(b) Maximum demand (W)");
+df=max_demand/LoadConnected;//demand factor
+disp(df,"(c) Demand factor");
+E=demand1*5+demand2*2+demand3*2+demand4*9+demand5*6;//Wh
+E=E/1000;//kWh
+disp(E,"(d) Energy consumed during 24 hours(kWh)")
+Edash=LoadConnected*24/1000;//kWh
+ disp(Edash,"(e) Energy consumed during 24 hours if all devices are used(kWh)")
diff --git a/2441/CH1/EX1.10/Ex1_10.sce b/2441/CH1/EX1.10/Ex1_10.sce
new file mode 100755
index 000000000..0022b9595
--- /dev/null
+++ b/2441/CH1/EX1.10/Ex1_10.sce
@@ -0,0 +1,12 @@
+//exa 1.10
+clc;clear;close;
+format('v',6);
+E=438*10^4;//kWh
+LF=20;//% annual
+CF=15;//%//Capacity Factor
+Lmax=E/(LF/100)/24/365;//kW
+Lmax=Lmax/1000;//MW
+C=Lmax/CF*LF;//MW//Plant Capacity
+disp(C,"Plant Capacity(MW): ");
+RC=C-Lmax;//MW//Reserve Capacity
+disp(RC,"Reserve Capacity(MW) : ");
diff --git a/2441/CH1/EX1.11/Ex1_11.sce b/2441/CH1/EX1.11/Ex1_11.sce
new file mode 100755
index 000000000..1c315543f
--- /dev/null
+++ b/2441/CH1/EX1.11/Ex1_11.sce
@@ -0,0 +1,21 @@
+//exa 1.11
+clc;clear;close;
+format('v',7);
+L1=10000;//kW
+L2=6000;//kW
+L3=8000;//kW
+L4=7000;//kW
+df=1.5;//diversity factor
+LF=65;//%//Load Factor
+Dinc=60;//%//Increase in maximum demand
+L=L1+L2+L3+L4;//kW//Sum
+L=L/1000;//MW
+Dmax=L/df;//MW
+disp(Dmax,"Maximum demand on station(MWh)");
+E=Dmax*365*24*LF/100;//MWh//Annual Energy
+format('v',9);
+disp(E,"Annual Energy Supplied(MWh)");
+Din_max=Dinc/100*Dmax;//MW
+format('v',7);
+C=Dmax+Din_max;//MW
+disp(C,"Installed Capacity(MW)")
diff --git a/2441/CH1/EX1.12/Ex1_12.sce b/2441/CH1/EX1.12/Ex1_12.sce
new file mode 100755
index 000000000..fe5bbf220
--- /dev/null
+++ b/2441/CH1/EX1.12/Ex1_12.sce
@@ -0,0 +1,28 @@
+//exa 1.12
+clc;clear;close;
+format('v',5);
+//Arranging data for Load Duration Curve
+//week days 5-9pm load
+L1=350;//MW
+t1=4*5;//hours
+//week days 8-12am & 1-5pm load
+L2=250;//MW
+t2=t1+8*5;//hours
+//saturday & sunday 5-9pm load
+L3=200;//MW
+t3=t2+4*2;//hours
+//All days 150MW load
+L4=150;//MW
+t4=t3+6*5+15*2;//hours
+//All days 100MW load
+L5=100;//MW
+t5=t4+6*5+5*2;//hours
+A=31600;//Total Load Curve Area
+LF=A/L1/24/7*100;//%//Weekly load factor
+disp(LF,"Weekly Load factor(%)");
+disp("Load Duration Curve is shown in figure.");
+//Load Duration Curve
+L=[L1 L2 L3 L4 L5];//MW
+T=[t1 t2 t3 t4 t5];//hours
+plot2d2(T,L);
+xtitle('Load Duration Curve','Time(Hours)','Load(MW)');
diff --git a/2441/CH1/EX1.13/Ex1_13.sce b/2441/CH1/EX1.13/Ex1_13.sce
new file mode 100755
index 000000000..b73a10cb6
--- /dev/null
+++ b/2441/CH1/EX1.13/Ex1_13.sce
@@ -0,0 +1,8 @@
+//exa 1.13
+clc;clear;close;
+format('v',7);
+LF=0.825;//Daily Load Factor
+ratio1=0.87;//daily peak load to monthly peak load
+ratio2=0.78;//monthly peak load to annually peak load
+LF_annual=LF*ratio1*ratio2;//Annual Load Factor
+disp(LF_annual,"Annual Load Factor : ");
diff --git a/2441/CH1/EX1.14/Ex1_14.sce b/2441/CH1/EX1.14/Ex1_14.sce
new file mode 100755
index 000000000..928950255
--- /dev/null
+++ b/2441/CH1/EX1.14/Ex1_14.sce
@@ -0,0 +1,30 @@
+//exa 1.14
+clc;clear;close;
+format('v',5);
+//Transformer1
+Lm=300;//kW
+df_m=0.6;//demand factor
+Lc=100;//kW//Commercial Load
+df_c=0.5;//demand factor
+//Transformer2
+Lr2=500;//kW//Residential Load
+df_Lr2=0.4;//demand factor
+//Transformer3
+Lr3=400;//kW
+df_Lr3=0.5;//demand factor
+//Diversity factors
+df1=2.3;
+df2=2.5;
+df3=2;
+DF=1.4;//Diversity factor between transformers
+//Solution :
+disp("Part(a)");
+Lp1=(Lm*df_m+Lc*df_c)/df1;//kW//Peak load on Transformer1
+disp(Lp1,"Peak load on Transformer1(kW)");
+Lp2=Lr2*df_Lr2/df2;//kW//Peak load on Transformer2
+disp(Lp2,"Peak load on Transformer2(kW)");
+Lp3=Lr3*df_Lr3/df3;//kW//Peak load on Transformer3
+disp(Lp3,"Peak load on Transformer3(kW)");
+disp("Part(b)");
+LpF=(Lp1+Lp2+Lp3)/DF;//Peak load on feeder
+disp(LpF,"Peak load on feeder(kW)");
diff --git a/2441/CH1/EX1.16/Ex1_16.sce b/2441/CH1/EX1.16/Ex1_16.sce
new file mode 100755
index 000000000..061f416d8
--- /dev/null
+++ b/2441/CH1/EX1.16/Ex1_16.sce
@@ -0,0 +1,16 @@
+//exa 1.16
+clc;clear;close;
+format('v',8);
+L=[20 25 30 25 35 20];//MW
+T=[6 4 2 4 4 4];//Hours
+Lmax=max(L);//MW
+disp(Lmax,"(a) Maximum demand(MW)");
+E=L(1)*sum(T)+(L(2)-L(1))*T(2)+(L(3)-L(1))*T(3)+(L(4)-L(1))*T(4)+(L(5)-L(1))*T(5)+(L(6)-L(1))*T(6);//MWh
+E=E*1000;//kWh
+disp(E,"(b) Units generated per day(kWh)");
+Lavg=E/sum(T);///kWh
+Lavg=Lavg/1000;///MW
+disp(Lavg,"(c) Average Load(MW)");
+format('v',6);
+LF=Lavg/Lmax*100;//%
+disp(LF,"(d) Load Factor(%)");
diff --git a/2441/CH1/EX1.17/Ex1_17.sce b/2441/CH1/EX1.17/Ex1_17.sce
new file mode 100755
index 000000000..a9369a035
--- /dev/null
+++ b/2441/CH1/EX1.17/Ex1_17.sce
@@ -0,0 +1,9 @@
+//exa 1.17
+clc;clear;close;
+format('v',8);
+pf=0.8;//power factor
+delf=1;//%//drop in frequency(delf/f)
+//delP=-2*(sind(theta))^2*delf
+theta=acosd(pf);//degree
+delP_BY_delf=-2*sind(theta)^2;//increase in load wrt frequency
+disp(-delP_BY_delf,"1% drop in frequency, Increased in Load(%)");
diff --git a/2441/CH1/EX1.18/Ex1_18.sce b/2441/CH1/EX1.18/Ex1_18.sce
new file mode 100755
index 000000000..f26b26a91
--- /dev/null
+++ b/2441/CH1/EX1.18/Ex1_18.sce
@@ -0,0 +1,8 @@
+//exa 1.18
+clc;clear;close;
+format('v',8);
+Lmax=100;//MW
+LF=40;//%//Load Factor
+Lavg=Lmax*LF/100;//MW
+E=Lavg*24*365;//MWh
+disp(E,"Energy generated in a year(MWh)");
diff --git a/2441/CH1/EX1.19/Ex1_19.sce b/2441/CH1/EX1.19/Ex1_19.sce
new file mode 100755
index 000000000..32ac5cd0f
--- /dev/null
+++ b/2441/CH1/EX1.19/Ex1_19.sce
@@ -0,0 +1,23 @@
+//exa 1.19 page 25
+clc;clear;close;
+format('v',5);
+V=400;//V
+s1=0.03;//initial slip
+delV=1;//%///Voltage Drop
+R1=0.290;//ohm/phase
+R2=0.15;//ohm/phase
+X=0.7;//ohm/phase(X1+X2)
+//V1^2*s1=V2^2*s2 for speed independent torque
+//taking for calculating s2
+V1=1;//V
+V2=V1-V1*delV/100;//V
+s2=V1^2/V2^2*s1;//slip
+I2ByI1=sqrt([R1+R2/s1]^2+X^2)/sqrt([R1+R2/s2]^2+X^2)*(V2/V1)
+delI=(I2ByI1-1)*100;//%//Current Increase
+disp(delI,"1% drop in Voltage increases current by(%)");
+//P=(R1+R2/s)*I^2
+P2ByP1=(R1+R2/s2)/(R1+R2/s1)*I2ByI1^2;//ratio
+delP=(1-P2ByP1)*100;//%//Power Decrease
+format('v',4);
+disp(delP,"1% drop in Voltage decreases power input by(%)");
+//Answer in the textbook is not accurate.
diff --git a/2441/CH1/EX1.2/Ex1_2.sce b/2441/CH1/EX1.2/Ex1_2.sce
new file mode 100755
index 000000000..0b32cf684
--- /dev/null
+++ b/2441/CH1/EX1.2/Ex1_2.sce
@@ -0,0 +1,46 @@
+//exa 1.2
+clc;clear;close;
+format('v',6);
+LoadA=2.5*1000;//W
+//12 midnight to 5am
+d1A=100;//W
+//5am to 6am
+d2A=1.1*1000;//W
+//6am to 8am
+d3A=200;//W
+//8am to 5pm
+d4A=0;//W
+//5pm to 12 midnight
+d5A=500;//W
+LoadB=3*1000;//W
+//11 pm to 7am
+d1B=0;//W
+//7 am to 8 am
+d2B=300;//W
+//8 am to 10 am
+d3B=1*1000;//W
+//10 am to 6 pm
+d4B=200;//W
+//6 pm to 11 pm
+d5B=600;//W
+DEMAND_A=[d1A d2A d3A d4A d5A];//W
+DEMAND_B=[d1B d2B d3B d4B d5B];//W
+max_demand_A=max(DEMAND_A);//W
+max_demand_B=max(DEMAND_B);//W
+df_A=max_demand_A/LoadA;//demand factor
+df_B=max_demand_B/LoadB;//demand factor
+disp(df_B,df_A,"Demand factor of consumer A & B are");
+gd_factor=(max_demand_A+max_demand_B)/max_demand_A;
+disp(gd_factor,"Group diversity factor")
+E_A=d1A*5+d2A*1+d3A*2+d4A*9+d5A*7;//Wh
+E_B=d1B*8+d2B*1+d3B*2+d4B*8+d5B*5;//Wh
+E_A=E_A/1000;//kWh
+E_B=E_B/1000;//kWh
+disp(E_B,E_A,"Energy consumed by A & B during 24 hours(kWh)")
+Emax_A=max_demand_A*24/1000;//kWh
+Emax_B=max_demand_B*24/1000;//kWh
+ disp(Emax_B,Emax_A,"Maximum energy consumer A & B can consume during 24 hours(kWh)")
+ratio_A=E_A/Emax_A;
+format('v',7);
+ratio_B=E_B/Emax_B;
+disp(ratio_B,ratio_A,"Ratio of actual energy to maximum energy of consumer A & B")
diff --git a/2441/CH1/EX1.3/Ex1_3.sce b/2441/CH1/EX1.3/Ex1_3.sce
new file mode 100755
index 000000000..df7fff87e
--- /dev/null
+++ b/2441/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,46 @@
+//exa 1.3
+clc;clear;close;
+format('v',6);
+n1=600;//No. of apartments
+L1=5;//kW//Each Apartment Load
+n2=20;//No. of general purpose shops
+L2=2;//kW//Each Shop Load
+df=0.8;//demand factor
+//1 Floor mill
+L3=10;//kW//Load
+df3=0.7;//demand factor
+//1 Saw mill
+L4=5;//kW//Load
+df4=0.8;//demand factor
+//1 Laundry
+L5=20;//kW//Load
+df5=0.65;//demand factor
+//1 Cinema
+L6=80;//kW//Load
+df6=0.5;//demand factor
+//Street lights
+n7=200;//no. of tube lights
+L7=40;//W//Load of each light
+//Residential Load
+df8=0.5;//demand factor
+gdf_r=3;//group diversity factor
+pdf_r=1.25;//peak diversity factor
+//Commertial Load
+gdf_c=2;//group diversity factor
+pdf_c=1.6;//peak diversity factor
+//Solution :
+//Maximum demand of each apartment
+dmax_1a=L1*df8;//kW
+//Maximum demand of 600 apartment
+dmax_a=n1*dmax_1a/gdf_r;//kW
+//demand of apartments at system peak time
+d_a_sp=dmax_a/pdf_r;//kW
+//Maximum Commercial demand
+dmax_c=(n2*L2*df+L3*df3+L4*df4+L5*df5+L6*df6)/gdf_c;//kW
+//Commercial demand at system peak time
+d_c_sp=dmax_c/pdf_c;//kW
+//demand of street light at system peak time
+d_sl_sp=n7*L7/1000;//kW
+//Increase in system peak demand
+DI=d_a_sp+d_c_sp+d_sl_sp;//kW
+disp(DI,"Increase in system peak demand(kW)");
diff --git a/2441/CH1/EX1.4/Ex1_4.sce b/2441/CH1/EX1.4/Ex1_4.sce
new file mode 100755
index 000000000..c43a13918
--- /dev/null
+++ b/2441/CH1/EX1.4/Ex1_4.sce
@@ -0,0 +1,23 @@
+//exa 1.4
+clc;clear;close;
+format('v',6);
+//12 to 5 am
+L1=20;//MW
+t1=5;//hours
+//5 to 9 am
+L2=40;//MW
+t2=4;//hours
+//9 to 6 pm
+L3=80;//MW
+t3=9;//hours
+//6 to 10 pm
+L4=100;//MW
+t4=4;//hours
+//10 to 12 am
+L5=20;//MW
+t5=2;//hours
+//Energy Poduced in 24 hours
+E=L1*t1+L2*t2+L3*t3+L4*t4+L5*t5;//MWh
+disp(E,"Energy Supplied by the plant in 24 hours(MWh) :");
+LF=E/24;//%//Load Factor
+disp(LF,"Load Factor(%)");
diff --git a/2441/CH1/EX1.5/Ex1_5.sce b/2441/CH1/EX1.5/Ex1_5.sce
new file mode 100755
index 000000000..27ca0d311
--- /dev/null
+++ b/2441/CH1/EX1.5/Ex1_5.sce
@@ -0,0 +1,27 @@
+//exa 1.5
+clc;clear;close;
+format('v',6);
+C=125;//MW//Installed Capacity
+//12 to 5 am
+L1=20;//MW
+t1=5;//hours
+//5 to 9 am
+L2=40;//MW
+t2=4;//hours
+//9 to 6 pm
+L3=80;//MW
+t3=9;//hours
+//6 to 10 pm
+L4=100;//MW
+t4=4;//hours
+//10 to 12 am
+L5=20;//MW
+t5=2;//hours
+//Energy Poduced in 24 hours
+E=L1*t1+L2*t2+L3*t3+L4*t4+L5*t5;//MWh
+LF=E/24;//%//Load Factor
+CF=LF/C;//%//Capacity Factor
+disp(CF,"Capacity Factor(%) : ");
+UF=100/C;//%//Utilisation Factor
+disp(UF,"Utilisation Factor(%) : ");
+
diff --git a/2441/CH1/EX1.6/Ex1_6.sce b/2441/CH1/EX1.6/Ex1_6.sce
new file mode 100755
index 000000000..2672452e1
--- /dev/null
+++ b/2441/CH1/EX1.6/Ex1_6.sce
@@ -0,0 +1,47 @@
+//exa 1.6
+clc;clear;close;
+format('v',6);
+//12 to 5 am & 10 to 12 am
+L1=20;//MW
+E1=L1*24;//MWh
+//5 to 9 am
+L2=40;//MW
+E2=E1+(L2-L1)*17;//MWh
+//9 to 6 pm
+L3=80;//MW
+E3=E2+(L3-L2)*13;//MWh
+//6 to 10 pm
+L4=100;//MW
+E4=E3+(L4-L3)*4;//MWh
+//Plotting Energy load curve
+L=[0,L1,L2,L3,L4];//MW
+E=[0,E1,E2,E3,E4];//Mwh
+subplot(2,1,1)
+plot(E,L)
+xlabel('Energy(MWh)');
+ylabel('Load(MW)');
+title('Energy Load Curve');
+//Energy Supplied
+//Upto 5am
+t1=5;//hours
+E1=L1*t1;//MWh
+//Upto 9am
+t2=4;//hours
+E2=E1+L2*t2;//MWh
+//Upto 6pm
+t3=9;//hours
+E3=E2+L3*t3;//MWh
+//Upto 10pm
+t4=4;//hours
+E4=E3+L4*t4;//MWh
+//Upto 12pm
+t4=2;//hours
+E4=E3+L4*t4;//MWh
+//Plotting Mass curve
+T=[0,1,2,3,4];//MW
+E=[0,E1,E2,E3,E4];//Mwh
+subplot(2,1,2)
+plot(T,E)
+ylabel('Energy(MWh)');
+xlabel('0-1: 12-5am 1-2: 5-9am 2-3: 9-6pm 3-4: 6-10pm above4: 10-12pm');
+title('Mass Curve');
diff --git a/2441/CH1/EX1.7/Ex1_7.sce b/2441/CH1/EX1.7/Ex1_7.sce
new file mode 100755
index 000000000..fa7fe2be9
--- /dev/null
+++ b/2441/CH1/EX1.7/Ex1_7.sce
@@ -0,0 +1,14 @@
+//exa 1.7
+clc;clear;close;
+format('v',9);
+dmax=40;//MW//Maximum demand
+CF=0.5;//Capacity Factor
+UF=0.8;//Utilisation Factor
+LF=CF/UF;///Load Factor
+disp(LF,"(a) Load Factor : ");
+C=dmax/UF;//MW//Plant Capacity
+disp(C,"(b) Plant Capacity(MW) : ");
+RC=C-dmax;//MW//Reserve Capacity
+disp(RC,"(c) Reserve Capacity(MW) : ");
+p=dmax*LF*24*365;//MWh//Annual Energy Production
+disp(p,"(d) Annual Energy Production(MWh) : ");
diff --git a/2441/CH1/EX1.8/Ex1_8.sce b/2441/CH1/EX1.8/Ex1_8.sce
new file mode 100755
index 000000000..76edbe592
--- /dev/null
+++ b/2441/CH1/EX1.8/Ex1_8.sce
@@ -0,0 +1,28 @@
+//exa 1.8
+clc;clear;close;
+format('v',6);
+L1=50;//MW//Initial
+t1=5;//hours
+L2=50;//MW//5am
+t2=4;//hours
+L3=100;//MW//9am
+t3=9;//hours
+L4=100;//MW//6pm
+t4=2;//hours
+L5=150;//MW//8pm
+t5=2;//hours
+L6=80;//MW//10pm
+t6=2;//hours
+L7=50;//MW
+//Energy Required in 24 hours
+E=L1*t1+(L2+L3)/2*t2+(L3+L4)/2*t3+(L4+L5)/2*t4+(L5+L6)/2*t5+(L6+L1)/2*t6;//MWh
+disp(E,"Energy required in one day(MWh)");
+DLF=E/L5/24*100;//%//Daily Load Factor
+disp(DLF,"Daily Load Factor(%)");
+//Plotting load curve
+T=[0,1,2,3,4,5,6];//Slots
+L=[L1,L2,L3,L4,L5,L6,L7];//MW
+plot(T,L)
+ylabel('Load(MW)');
+xlabel('0-1: 12-5am 1-2: 5-9am 2-3: 9-6pm 3-4: 6-8pm 4-5:8-10pm 5-6 :10-12pm');
+title('Chronological Load Curve');
diff --git a/2441/CH1/EX1.9/Ex1_9.sce b/2441/CH1/EX1.9/Ex1_9.sce
new file mode 100755
index 000000000..d9ab81093
--- /dev/null
+++ b/2441/CH1/EX1.9/Ex1_9.sce
@@ -0,0 +1,51 @@
+//exa 1.9
+clc;clear;close;
+format('v',6);
+L1=50;//MW//Initial
+t1=5;//hours
+L2=50;//MW//5am
+t2=4;//hours
+L3=100;//MW//9am
+t3=9;//hours
+L4=100;//MW//6pm
+t4=2;//hours
+L5=150;//MW//8pm
+t5=2;//hours
+L6=80;//MW//10pm
+t6=2;//hours
+L7=50;//MW
+//Load Duration Curve
+l1=L5;//Mw
+l2=L4;//MW
+l3=L1;//MW
+L=[l1,l2,l2,l3,l3]
+T=0:6:24;//Duration in hours
+subplot(2,1,1)
+plot(T,L)
+ylabel('Load(MW)');
+xlabel('Hours');
+title('Load Duration Curve');
+//Energy Consumed
+//Upto 5am
+t1=5;//hours
+E1=L1*t1;//MWh
+//Upto 9am
+t2=4;//hours
+E2=E1+L2*t2;//MWh
+//Upto 6pm
+t3=9;//hours
+E3=E2+L3*t3;//MWh
+//Upto 10pm
+t4=4;//hours
+E4=E3+L4*t4;//MWh
+//Upto 12pm
+t4=2;//hours
+E4=E3+L4*t4;//MWh
+//Plotting Mass curve
+T=[0,1,2,3,4];//MW
+E=[0,E1,E2,E3,E4];//Mwh
+subplot(2,1,2)
+plot(T,E)
+ylabel('Energy(MWh)');
+xlabel('0-1: 12-5am 1-2: 5-9am 2-3: 9-6pm 3-4: 6-10pm above4: 10-12pm');
+title('Mass Curve');