function mpc = case_skript_BEE
%zugrundegelegt wird die Schaltung BEE-Skript, Beispiel 7
%Admittanzen
y12=5*exp(-j*2*pi*(69.047/360));
y13=9*exp(-j*2*pi*(55/360));
y14=9*exp(-j*2*pi*(55/360));
y23=4*exp(-j*2*pi*(69.047/360));
y34=6*exp(-j*2*pi*(69.047/360));
y130=j*0.63;
y140=j*0.63;


%% MATPOWER Case Format : Version 2
mpc.version = '2';

%%-----  Power Flow Data  -----%%
%% system MVA base
mpc.baseMVA = 100;

%% bus data
%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin
mpc.bus = [
	1	1	120	59.3	0	0	1	1	0	110	1	1.1	0.9;
	2	1	60	20	0	0	1	1	0	110	1	1.1	0.9;
	3	2	0	0	0	0	1	1	0	110	1	1.0	1.0;
	4	3	0	0	0	0	1	1	0	110	1	1.0	1.0;
];

%% 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 = [
	3	100	0	300	-300	1	100	1	500	0	0	0	0	0	0	0	0	0	0	0	0 0;
	4	100	0	300	-300	1	100	1	500 0	0	0	0	0	0	0	0	0	0	0	0 0;
];

%% branch data
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax
mpc.branch = [
	1	2	real(1/y12)	imag(1/y12)	0 500	500	500	0	0	1	-360	360;
	1	3	real(1/y13)	imag(1/y13)	imag(y130)	500	500	500	0	0	1	-360	360;
	1	4	real(1/y14)	imag(1/y14)	imag(y140)	500	500	500	0	0	1	-360	360;
	3	4	real(1/y34)	imag(1/y34)	0	500	500	500	0	0	1	-360	360;
    2	3	real(1/y23)	imag(1/y23)	0	500	500	500	0	0	1	-360	360;
];


%%-----  OPF Data  -----%%
%% generator cost data
%	1	startup	shutdown	n	x1	y1	...	xn	yn
%	2	startup	shutdown	n	c(n-1)	...	c0
mpc.gencost = [
	2	0	0	3	1e-5 0.005 6;
	2	0	0	3	5e-6 0.006 5;
];

