% Forty-Units Test System with Valve-Point Loading Effects & Emission-Rates
% This data is entered by: Ali R. Al-Roomi in the 1st of May, 2016
% Refer to http://www.al-roomi.org for more details about the system

% System Data 
n=40; % number of units
Load=10500; % Load Demand (in MW)
Pmin=[36;36;60;80;47;68;110;135;135;130;94;94;125;125;125;125;220;220;242;242;254;254;254;254;254;254;10;10;10;47;60;60;60;90;90;90;25;25;25;242]; % for P1, P2, ..., Pn (in MW)
Pmax=[114;114;120;190;97;140;300;300;300;300;375;375;500;500;500;500;500;500;550;550;550;550;550;550;550;550;150;150;150;97;190;190;190;200;200;200;110;110;110;550]; % for P1, P2, ..., Pn (in MW)

% Fuel-Cost Curves Coefficients
a=[94.705;94.705;309.54;369.03;148.89;222.33;287.71;391.98;455.76;722.82;635.20;654.69;913.40;1760.4;1728.3;1728.3;647.85;649.69;647.83;647.81;785.96;785.96;794.53; ...
    794.53;801.32;801.32;1055.1;1055.1;1055.1;148.89;222.92;222.92;222.92;107.87;116.58;116.58;307.45;307.45;307.45;647.83]; % for unit 1, 2, ..., n
b=[6.73;6.73;7.07;8.18;5.35;8.05;8.03;6.99;6.60;12.9;12.9;12.8;12.5;8.84;9.15;9.15;7.97;7.95;7.97;7.97;6.63;6.63;6.66;6.66;7.10;7.10;3.33;3.33;3.33;5.35;6.43;6.43; ...
    6.43;8.95;8.62;8.62;5.88;5.88;5.88;7.97]; % for unit 1, 2, ..., n
c=[0.00690;0.00690;0.02028;0.00942;0.01140;0.01142;0.00357;0.00492;0.00573;0.00605;0.00515;0.00569;0.00421;0.00752;0.00708;0.00708;0.00313;0.00313;0.00313;0.00313; ...
    0.00298;0.00298;0.00284;0.00284;0.00277;0.00277;0.52124;0.52124;0.52124;0.01140;0.00160;0.00160;0.00160;0.00010;0.00010;0.00010;0.01610;0.01610;0.01610;0.00313]; % for unit 1, 2, ..., n

% Valve-Point Loading Effects Coefficients
d=[100;100;100;150;120;100;200;200;200;200;200;200;300;300;300;300;300;300;300;300;300;300;300;300;300;300;120;120;120;120;150;150;150;200;200;200;80;80;80;300]; % for unit 1, 2, ..., n
e=[0.084;0.084;0.084;0.063;0.077;0.084;0.042;0.042;0.042;0.042;0.042;0.042;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.035;0.077; ...
    0.077;0.077;0.077;0.063;0.063;0.063;0.042;0.042;0.042;0.098;0.098;0.098;0.035]; % for unit 1, 2, ..., n

% Emission-Rates Coefficients
Alpha=[60;60;100;120;50;80;100;130;150;280;220;225;300;520;510;510;220;222;220;220;290;285;295;295;310;310;360;360;360;50;80;80;80;150;130;150;100;100;100;220]; % for unit 1, 2, ..., n
Beta=-1*[2.2200;2.2200;2.3600;3.1400;1.8900;3.0800;3.0600;2.3200;2.1100;4.3400;4.3400;4.2800;4.1800;3.3400;3.5500;3.5500;2.6800;2.6600;2.6800;2.6800;2.2200; ...
          2.2200;2.2600;2.2600;2.4200;2.4200;1.1100;1.1100;1.1100;1.8900;2.0800;2.0800;2.0800;2.0600;2.3200;2.1100;1.9800;1.9800;1.9800;2.6800]; % for unit 1, 2, ..., n
Gamma=[0.0480;0.0480;0.0762;0.0540;0.0850;0.0854;0.0242;0.0310;0.0335;0.4250;0.0322;0.0338;0.0296;0.0512;0.0496;0.0496;0.0151;0.0151;0.0151;0.0151;0.0145;0.0145; ...
        0.0138;0.0138;0.0132;0.01321.84201.84201.8420;0.0850;0.0121;0.0121;0.0121;0.0542;0.0410;0.0435;0.0950;0.0950;0.0950;0.0151]; % for unit 1, 2, ..., n
Xi=[1.3100;1.3100;1.3100;0.9142;0.9936;1.3100;0.6550;0.6550;0.6550;0.6550;0.6550;0.6550;0.5035;0.5035;0.5035;0.5035;0.5035;0.5035;0.5035;0.5035;0.5035;0.5035; ...
     0.5035;0.5035;0.5035;0.5035;0.9936;0.9936;0.9936;0.9936;0.9142;0.9142;0.9142;0.6550;0.6550;0.6550;1.4200;1.4200;1.4200;0.5035]; % for unit 1, 2, ..., n
Delta=[0.0569;0.0569;0.0569;0.0454;0.0406;0.0569;0.0285;0.0285;0.0285;0.0285;0.0285;0.0285;0.0208;0.0208;0.0208;0.0208;0.0208;0.0208;0.0208;0.0208;0.0208; ...
        0.0208;0.0208;0.0208;0.0208;0.0208;0.0406;0.0406;0.0406;0.0406;0.0454;0.0454;0.0454;0.0285;0.0285;0.0285;0.0677;0.0677;0.0677;0.0208]; % for unit 1, 2, ..., n