//Example 5.6
clc;clear;close;
format('v',6);
f=50;//Hz
R=2;//Hz/pu MW
Pr=10000;//MW(Rated Capacity)
P=Pr/2;//MW(Operating Power)
delP=2;//%(Load Increase)
del_f=f*1/100;//Hz(1% change in frequency)
del_PL=P*1/100;//MW(1% change in load)
//Rate of change of load with frequency :
D=del_PL/del_f;//MW/Hz
D=D/Pr;//p.u. MW/Hz
//Frequency response characteristic : 
Beta=D+1/R;//p.u. MW/Hz
M=delP/100*P;//MW
M=M/Pr;//p.u. MW
del_fo=-M/Beta;//Hz
disp("Frequency drop contribution to increase in load(MW) : ");
delP_fo=-del_fo*(D*Pr);//MW
disp(delP_fo);
disp("Increase in generation to meet the increase load(MW) ");
delP_gen=-del_fo/R*Pr;//MW
disp(delP_gen);