 Elements of prefault bus admittance matrix   
 
  - 11.235955i    0             0                         11.235955i                0                       
    0           - 7.1428571i    0                         0                         7.1428571i              
    0             0             11.284152 - 65.473163i  - 4.244997 + 24.257126i   - 7.0391553 + 41.355037i  
    11.235955i    0           - 4.244997 + 24.257126i     6.6587549 - 44.61749i   - 1.4488088 + 8.8538313i  
    0             7.1428571i  - 7.0391553 + 41.355037i  - 1.4488088 + 8.8538313i    8.9771429 - 57.297263i  



 
 After bus 4 is shorted the matrix becomes   
 
  - 11.235955i    0             0                         0                       
    0           - 7.1428571i    0                         7.1428571i              
    0             0             11.284152 - 65.473163i  - 7.0391553 + 41.355037i  
    0             7.1428571i  - 7.0391553 + 41.355037i    8.9771429 - 57.297263i  



 
 Elements of faulted bus admittance matrices   
 
  - 11.235955i    0                         0                       
    0             0.1361702 - 6.2737408i  - 0.0681128 + 5.1661208i  
    0           - 0.0681128 + 5.1661208i    5.7986816 - 35.62997i   

 Pe1 = 0 
 Pe2 = 0.1544 + 5.5012sin(delta_2 - 0.755) 

 Swing Equation for machine 1  on fault in elec deg/square sec 
 d^2(delta_1)/dt^2 = (180f/11.2) * (3.5) 

 Swing equation for machine 2 on fault in elec deg/square sec 
 d^2(delta_2)/dt^2 = (180f/8.0) * (1.6956 - 5.5012sin(delta_2 - 0.755)

