param :    B:       R          	X       	Imax :=
#      from to      (Ω)         (Ω)       	(A)   	Branches		Start and end buses, Phase resistance and reoctance, Maximum line current
		1	2		0.0922		0.0477		500		#1
		2	3		0.4930		0.2511		500		#2
		3	4		0.3660		0.1864		500		#3
		4	5		0.3811		0.1941		500		#4
		5	6		0.8190		0.7070		500		#5
		6	7		0.1872		0.6188		500		#6
		7	8		0.7114		0.2351		500		#7
		8	9		1.0300		0.7400		500		#8
		9	10		1.0440		0.7400		500		#9
		10	11		0.1966		0.0650		500		#10
		11	12		0.3744		0.1238		500		#11
		12	13		1.4680		1.1550		500		#12
		13	14		0.5416		0.7129		500		#13
		14	15		0.5910		0.5260		500		#14
		15	16		0.7463		0.5450		500		#15
		16	17		1.2890		1.7210		500		#16
		17	18		0.7320		0.5740		500		#17
		2	19		0.1640		0.1565		500		#18
		19	20		1.5042		1.3554		500		#19
		20	21		0.4095		0.4784		500		#20
		21	22		0.7089		0.9373		500		#21
		3	23		0.4512		0.3083		500		#22
		23	24		0.8980		0.7091		500		#23
		24	25		0.8960		0.7011		500		#24
		6	26		0.2030		0.1034		500		#25
		26	27		0.2842		0.1447		500		#26
		27	28		1.0590		0.9337		500		#27
		28	29		0.8042		0.7006		500		#28
		29	30		0.5075		0.2585		500		#29
		30	31		0.9744		0.9630		500		#30
		31	32		0.3105		0.3619		500		#31
		32	33		0.3410		0.5302		500		#32
#		Tie Lines
#      from to        (Ω)         (Ω)       (A)   Branches
		8	21		2.0000		2.0000		500		#33
		9	15		2.0000		2.0000		500		#34
		12	22		2.0000		2.0000		500		#35
		18	33		0.5000		0.5000		500		#36
		25	29		0.5000		0.5000		500;	#37

param : N:  Tb      PD      QD  PDNT  QDNT :=
#      Bus Type    (kW)   (kVAr)		Bus number, Bus type, 3-phase active and reactive demand
		1	1		0		0     0  0
		2	0		100		60    1  0.6
		3	0		90		40    0.9  0.4
		4	0		120		80    0  0
		5	0		60		30    0.6  0.3
		6	0		60		20    0  0
		7	0		200		100   0  0
		8	0		200		100   0  0
		9	0		60		20    0.6  0.2
		10	0		60		20    0  0
		11	0		45		30    0  0
		12	0		60		35    0  0
		13	0		60		35    0  0
		14	0		120		80    0  0
		15	0		60		10    0  0
		16	0		60		20    0  0
		17	0		60		20    0  0
		18	0		90		40    0  0
		19	0		90		40    0  0
		20	0		90		40    0  0
		21	0		90		40    0  0
		22	0		90		40    0  0
		23	0		90		50    0.9  0.5
		24	0		420		200   4.2  2
		25	0		420		200   0  0
		26	0		60		25    0  0
		27	0		60		25    0  0
		28	0		60		20    0  0
		29	0		120		70    0  0
		30	0		200		600   0  0
		31	0		150		70    0  0
		32	0		210		100   0  0
		33	0		60		40    0.6  0.4;

param Vnom  := 12.66;		#Nominal line voltage (kV)
param Vbase := 12.66;		#Base line voltage (kV)
param Sbase := 1000;		#3-Phase base power (kVA)=10000 or 100000