% initial operating point
% voltage magnitudes in pu
% voltage phase angles in radian
% powers in MW / Mvar; positive if injected into the network

opt(  1).name= '1       ' ;
opt(  1).magv=   1.0200000     ;
opt(  1).phav=   0.0000000     ;
opt(  1).injp=  -295.18842     ;
opt(  1).injq=   38.516266     ;

opt(  2).name= '2       ' ;
opt(  2).magv=   1.0037710     ;
opt(  2).phav=  0.32197397E-01 ;
opt(  2).injp=  -150.00000     ;
opt(  2).injq=  -30.000002     ;

opt(  3).name= '3       ' ;
opt(  3).magv=   1.0153710     ;
opt(  3).phav=  0.11973300     ;
opt(  3).injp=   0.0000000     ;
opt(  3).injq=   0.0000000     ;

opt(  4).name= '4       ' ;
opt(  4).magv=   1.0183720     ;
opt(  4).phav=  0.14718790     ;
opt(  4).injp=   0.0000000     ;
opt(  4).injq=   0.0000000     ;

opt(  5).name= '5       ' ;
opt(  5).magv=   1.0000000     ;
opt(  5).phav=  0.28415780     ;
opt(  5).injp=   450.00000     ;
opt(  5).injq=   68.497604     ;

% Y0 nodal admittance matrix
% base power : 100 MVA
% base voltages : 225 kV, 20 kV
% Contributions of lines and transformers only

Y0=zeros(  5) ;

Y0(  1,  3)=Y0(  1,  3)+  -1.2048982     +   12.048983     *j ;
Y0(  3,  1)=Y0(  3,  1)+  -1.2048982     +   12.048983     *j ;
Y0(  1,  1)=Y0(  1,  1)+   1.2048982     +  -11.950264     *j ;
Y0(  3,  3)=Y0(  3,  3)+   1.2048982     +  -11.950264     *j ;
Y0(  1,  3)=Y0(  1,  3)+  -1.2048982     +   12.048983     *j ;
Y0(  3,  1)=Y0(  3,  1)+  -1.2048982     +   12.048983     *j ;
Y0(  1,  1)=Y0(  1,  1)+   1.2048982     +  -11.950264     *j ;
Y0(  3,  3)=Y0(  3,  3)+   1.2048982     +  -11.950264     *j ;
Y0(  3,  4)=Y0(  3,  4)+  -7.8318372     +   78.318375     *j ;
Y0(  4,  3)=Y0(  4,  3)+  -7.8318372     +   78.318375     *j ;
Y0(  3,  3)=Y0(  3,  3)+   7.8318372     +  -78.303185     *j ;
Y0(  4,  4)=Y0(  4,  4)+   7.8318372     +  -78.303185     *j ;
Y0(  3,  4)=Y0(  3,  4)+  -7.8318372     +   78.318375     *j ;
Y0(  4,  3)=Y0(  4,  3)+  -7.8318372     +   78.318375     *j ;
Y0(  3,  3)=Y0(  3,  3)+   7.8318372     +  -78.303185     *j ;
Y0(  4,  4)=Y0(  4,  4)+   7.8318372     +  -78.303185     *j ;
Y0(  5,  4)=Y0(  5,  4)+   0.0000000     +   32.362461     *j ;
Y0(  4,  5)=Y0(  4,  5)+   0.0000000     +   32.362461     *j ;
Y0(  5,  5)=Y0(  5,  5)+   0.0000000     +  -33.333332     *j ;
Y0(  4,  4)=Y0(  4,  4)+   0.0000000     +  -31.419863     *j ;
Y0(  2,  3)=Y0(  2,  3)+   0.0000000     +   16.835016     *j ;
Y0(  3,  2)=Y0(  3,  2)+   0.0000000     +   16.835016     *j ;
Y0(  2,  2)=Y0(  2,  2)+   0.0000000     +  -16.666666     *j ;
Y0(  3,  3)=Y0(  3,  3)+   0.0000000     +  -17.005066     *j ;

