function mpc = FourAreaNetwork																					
mpc.version = '2';																					
mpc.baseMVA = 100;																					

%% Bus Data																					
%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin								
mpc.bus = [																					
    1   	2   	0	0	0	0	1	1	0	22	1	1.05	0.95;
    2   	1   	2800	0	0	0	1	1	0	400	1	1.05	0.95;
    3   	1   	1700	0	0	0	1	1	0	400	1	1.05	0.95;
    4   	2   	0	0	0	0	1	1	0	22	1	1.05	0.95;
    5   	2   	0	0	0	0	1	1	0	13.8	1	1.05	0.95;
    6   	1   	1300	0	0	0	1	1	0	400	1	1.05	0.95;
    7   	1   	1700	0	0	0	1	1	0	400	1	1.05	0.95;
    8   	3   	0	0	0	0	1	1	0	13.8	1	1.05	0.95;								
];																					
%% Generator Data																					
%	bus	Pg	Qg	Qmax	Qmin	Vg	mBase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	ramp_agc	ramp_10	ramp_30	ramp_q	apf
mpc.gen = [																					
    1	1740	0	1263	-1263	1	100	1	2608.2	0	0	0	0	0	0	0	0	0	0	0	0;
    4	1800	0	1263	-1263	1	100	1	2608.2	0	0	0	0	0	0	0	0	0	0	0	0;
    5	1960	0	1142.03	-1142.03	1	100	1	2358	0	0	0	0	0	0	0	0	0	0	0	0;
    8	2019	0	1256.23	-1256.23	1	100	1	2593.8	0	0	0	0	0	0	0	0	0	0	0	0;
];																					
%% Transmission Lines Data																					
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax								
mpc.branch = [																					
 
    2	3	0.001875	0.0265625	0	400	400	400	0	0	1	-360	360.00;
    2	6	0.0075	0.10625	0	400	400	400	0	0	1	-360	360.00;
    2	7	0.0075	0.10625	0	400	400	400	0	0	1	-360	360.00;
    
    3	6	0.001875	0.0265625	0	400	400	400	0	0	1	-360	360.00;
    3	7	0.001875	0.0265625	0	400	400	400	0	0	1	-360	360.00;
   
    6	7	0.0075	0.10625	0	400	400	400	0	0	1	-360	360.00;
   
    
    1	2	0	0.001	0	3000	3000	3000	0	0	1	-360	360.00;
    3	4	0	0.001	0	3000	3000	3000	0	0	1	-360	360.00;
    5	6	0	0.001	0	3000	3000	3000	0	0	1	-360	360.00;
    7	8	0	0.001	0	3000	3000	3000	0	0	1	-360	360.00;
];																					
mpc.areas =[																					
1 1;																					
];																					
%% Generator Cost Data (Made up)																					
%	1	startup	shutdown	n	x1	y1	...	xn	yn												
%	2	startup	shutdown	n	c(n-1)	...	c0														
mpc.gencost = [																					
2	1500	0	3	0	0.1	0.05;															
2	1500	0	3	0	0.1	0.01;															
2	1500	0	3	0	0.1	0.5;															
2	1500	0	3	0	0.1	0.02;																												
];																					
