 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 removing line 4-5   
 
  - 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     5.2099462 - 35.876658i    0                       
    0             7.1428571i  - 7.0391553 + 41.355037i    0                         7.5283342 - 48.556432i  






 
 Elements of post faulted bus admittance matrices   
 
    0.5004560 - 7.7897217i    0                       - 0.2215766 + 7.6290933i  
    0                         0.1590864 - 6.1167778i  - 0.0901191 + 6.0974738i  
  - 0.2215766 + 7.6290933i  - 0.0901191 + 6.0974738i    1.3927745 - 13.872781i  

 Pe1 = 0.6057 + 8.3965sin(delta_1 - 1.664) 

 Pe2 = 0.1804 + 6.4932sin(delta_2 - 0.847) 

 Swing equation for machine 1 on fault in elec deg/square sec
 d^2(delta_2)/dt^2 = (180f/11.2) * (2.8943 - 8.3965sin(delta_1 - 1.664)

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

