Please go through and let me know if my calculations are cor

Solution

There is a problem with calculating Z bus

Ybus=

>> a=[-j*13.33 j*10 0 0; j*100 0.1098-j*11.23 -0.1089+j*1.2452 0; 0 -0.1089+j*1.2452 0.1098-j*11.23 j*10;0 0 j*10 -j*13.33 ]

a =

1.0e+002 *

0 - 0.1333i 0 + 0.1000i 0 0
0 + 1.0000i 0.0011 - 0.1123i -0.0011 + 0.0125i 0
0 -0.0011 + 0.0125i 0.0011 - 0.1123i 0 + 0.1000i
0 0 0 + 0.1000i 0 - 0.1333i

>> inv(a)

ans =

0.0001 - 0.0126i 0.0000 - 0.0117i 0.0002 - 0.0039i 0.0002 - 0.0029i
0.0001 - 0.1168i 0.0000 - 0.0156i 0.0003 - 0.0052i 0.0002 - 0.0039i
0.0023 - 0.0391i 0.0003 - 0.0052i 0.0081 + 0.2663i 0.0061 + 0.1997i
0.0017 - 0.0293i 0.0002 - 0.0039i 0.0061 + 0.1997i 0.0046 + 0.2249i

so your Zbus is not matching with YBus