


  mismatch(0+1)=delta_t1(0+1)-delta(0+1)-.50000000000000001040834085586084e-3*delta_dot_t1(0+1)-.50000000000000001040834085586084e-3*delta_dot(0+1); % Depth 1, Step 1
  mismatch(1+1)=delta_t1(1+1)-delta(1+1)-.50000000000000001040834085586084e-3*delta_dot_t1(1+1)-.50000000000000001040834085586084e-3*delta_dot(1+1); % Depth 1, Step 2
  mismatch(2+1)=delta_t1(2+1)-delta(2+1)-.50000000000000001040834085586084e-3*delta_dot_t1(2+1)-.50000000000000001040834085586084e-3*delta_dot(2+1); % Depth 1, Step 3
  mismatch(3+1)=delta_t1(3+1)-delta(3+1)-.50000000000000001040834085586084e-3*delta_dot_t1(3+1)-.50000000000000001040834085586084e-3*delta_dot(3+1); % Depth 1, Step 4
  mismatch(4+1)=delta_t1(4+1)-delta(4+1)-.50000000000000001040834085586084e-3*delta_dot_t1(4+1)-.50000000000000001040834085586084e-3*delta_dot(4+1); % Depth 1, Step 5
  mismatch(5+1)=delta_t1(5+1)-delta(5+1)-.50000000000000001040834085586084e-3*delta_dot_t1(5+1)-.50000000000000001040834085586084e-3*delta_dot(5+1); % Depth 1, Step 6
  mismatch(6+1)=delta_t1(6+1)-delta(6+1)-.50000000000000001040834085586084e-3*delta_dot_t1(6+1)-.50000000000000001040834085586084e-3*delta_dot(6+1); % Depth 1, Step 7
  mismatch(7+1)=delta_t1(7+1)-delta(7+1)-.50000000000000001040834085586084e-3*delta_dot_t1(7+1)-.50000000000000001040834085586084e-3*delta_dot(7+1); % Depth 1, Step 8
  mismatch(8+1)=delta_t1(8+1)-delta(8+1)-.50000000000000001040834085586084e-3*delta_dot_t1(8+1)-.50000000000000001040834085586084e-3*delta_dot(8+1); % Depth 1, Step 9
  mismatch(9+1)=delta_t1(9+1)-delta(9+1)-.50000000000000001040834085586084e-3*delta_dot_t1(9+1)-.50000000000000001040834085586084e-3*delta_dot(9+1); % Depth 1, Step 10
  mismatch(10+1)=delta_t1(10+1)-delta(10+1)-.50000000000000001040834085586084e-3*delta_dot_t1(10+1)-.50000000000000001040834085586084e-3*delta_dot(10+1); % Depth 1, Step 11
  mismatch(11+1)=delta_t1(11+1)-delta(11+1)-.50000000000000001040834085586084e-3*delta_dot_t1(11+1)-.50000000000000001040834085586084e-3*delta_dot(11+1); % Depth 1, Step 12
  mismatch(12+1)=delta_t1(12+1)-delta(12+1)-.50000000000000001040834085586084e-3*delta_dot_t1(12+1)-.50000000000000001040834085586084e-3*delta_dot(12+1); % Depth 1, Step 13
  mismatch(13+1)=delta_t1(13+1)-delta(13+1)-.50000000000000001040834085586084e-3*delta_dot_t1(13+1)-.50000000000000001040834085586084e-3*delta_dot(13+1); % Depth 1, Step 14
  mismatch(14+1)=delta_t1(14+1)-delta(14+1)-.50000000000000001040834085586084e-3*delta_dot_t1(14+1)-.50000000000000001040834085586084e-3*delta_dot(14+1); % Depth 1, Step 15
  mismatch(15+1)=delta_t1(15+1)-delta(15+1)-.50000000000000001040834085586084e-3*delta_dot_t1(15+1)-.50000000000000001040834085586084e-3*delta_dot(15+1); % Depth 1, Step 16
  mismatch(16+1)=delta_t1(16+1)-delta(16+1)-.50000000000000001040834085586084e-3*delta_dot_t1(16+1)-.50000000000000001040834085586084e-3*delta_dot(16+1); % Depth 1, Step 17
  mismatch(17+1)=delta_t1(17+1)-delta(17+1)-.50000000000000001040834085586084e-3*delta_dot_t1(17+1)-.50000000000000001040834085586084e-3*delta_dot(17+1); % Depth 1, Step 18
  mismatch(18+1)=delta_t1(18+1)-delta(18+1)-.50000000000000001040834085586084e-3*delta_dot_t1(18+1)-.50000000000000001040834085586084e-3*delta_dot(18+1); % Depth 1, Step 19
  mismatch(19+1)=delta_t1(19+1)-delta(19+1)-.50000000000000001040834085586084e-3*delta_dot_t1(19+1)-.50000000000000001040834085586084e-3*delta_dot(19+1); % Depth 1, Step 20
  mismatch(20+1)=omega_t1(0+1)-omega(0+1)-.50000000000000001040834085586084e-3*omega_dot_t1(0+1)-.50000000000000001040834085586084e-3*omega_dot(0+1); % Depth 1, Step 21
  mismatch(21+1)=omega_t1(1+1)-omega(1+1)-.50000000000000001040834085586084e-3*omega_dot_t1(1+1)-.50000000000000001040834085586084e-3*omega_dot(1+1); % Depth 1, Step 22
  mismatch(22+1)=omega_t1(2+1)-omega(2+1)-.50000000000000001040834085586084e-3*omega_dot_t1(2+1)-.50000000000000001040834085586084e-3*omega_dot(2+1); % Depth 1, Step 23
  mismatch(23+1)=omega_t1(3+1)-omega(3+1)-.50000000000000001040834085586084e-3*omega_dot_t1(3+1)-.50000000000000001040834085586084e-3*omega_dot(3+1); % Depth 1, Step 24
  mismatch(24+1)=omega_t1(4+1)-omega(4+1)-.50000000000000001040834085586084e-3*omega_dot_t1(4+1)-.50000000000000001040834085586084e-3*omega_dot(4+1); % Depth 1, Step 25
  mismatch(25+1)=omega_t1(5+1)-omega(5+1)-.50000000000000001040834085586084e-3*omega_dot_t1(5+1)-.50000000000000001040834085586084e-3*omega_dot(5+1); % Depth 1, Step 26
  mismatch(26+1)=omega_t1(6+1)-omega(6+1)-.50000000000000001040834085586084e-3*omega_dot_t1(6+1)-.50000000000000001040834085586084e-3*omega_dot(6+1); % Depth 1, Step 27
  mismatch(27+1)=omega_t1(7+1)-omega(7+1)-.50000000000000001040834085586084e-3*omega_dot_t1(7+1)-.50000000000000001040834085586084e-3*omega_dot(7+1); % Depth 1, Step 28
  mismatch(28+1)=omega_t1(8+1)-omega(8+1)-.50000000000000001040834085586084e-3*omega_dot_t1(8+1)-.50000000000000001040834085586084e-3*omega_dot(8+1); % Depth 1, Step 29
  mismatch(29+1)=omega_t1(9+1)-omega(9+1)-.50000000000000001040834085586084e-3*omega_dot_t1(9+1)-.50000000000000001040834085586084e-3*omega_dot(9+1); % Depth 1, Step 30
  mismatch(30+1)=omega_t1(10+1)-omega(10+1)-.50000000000000001040834085586084e-3*omega_dot_t1(10+1)-.50000000000000001040834085586084e-3*omega_dot(10+1); % Depth 1, Step 31
  mismatch(31+1)=omega_t1(11+1)-omega(11+1)-.50000000000000001040834085586084e-3*omega_dot_t1(11+1)-.50000000000000001040834085586084e-3*omega_dot(11+1); % Depth 1, Step 32
  mismatch(32+1)=omega_t1(12+1)-omega(12+1)-.50000000000000001040834085586084e-3*omega_dot_t1(12+1)-.50000000000000001040834085586084e-3*omega_dot(12+1); % Depth 1, Step 33
  mismatch(33+1)=omega_t1(13+1)-omega(13+1)-.50000000000000001040834085586084e-3*omega_dot_t1(13+1)-.50000000000000001040834085586084e-3*omega_dot(13+1); % Depth 1, Step 34
  mismatch(34+1)=omega_t1(14+1)-omega(14+1)-.50000000000000001040834085586084e-3*omega_dot_t1(14+1)-.50000000000000001040834085586084e-3*omega_dot(14+1); % Depth 1, Step 35
  mismatch(35+1)=omega_t1(15+1)-omega(15+1)-.50000000000000001040834085586084e-3*omega_dot_t1(15+1)-.50000000000000001040834085586084e-3*omega_dot(15+1); % Depth 1, Step 36
  mismatch(36+1)=omega_t1(16+1)-omega(16+1)-.50000000000000001040834085586084e-3*omega_dot_t1(16+1)-.50000000000000001040834085586084e-3*omega_dot(16+1); % Depth 1, Step 37
  mismatch(37+1)=omega_t1(17+1)-omega(17+1)-.50000000000000001040834085586084e-3*omega_dot_t1(17+1)-.50000000000000001040834085586084e-3*omega_dot(17+1); % Depth 1, Step 38
  mismatch(38+1)=omega_t1(18+1)-omega(18+1)-.50000000000000001040834085586084e-3*omega_dot_t1(18+1)-.50000000000000001040834085586084e-3*omega_dot(18+1); % Depth 1, Step 39
  mismatch(39+1)=omega_t1(19+1)-omega(19+1)-.50000000000000001040834085586084e-3*omega_dot_t1(19+1)-.50000000000000001040834085586084e-3*omega_dot(19+1); % Depth 1, Step 40
  mismatch(40+1)=P_load(0+1)+V(0+1)*(Bus_Ymag(0+1)*cos(Bus_Yang(0+1))*V(0+1)+Bus_Ymag(3+1)*cos(-theta(0+1)+theta(1+1)+Bus_Yang(3+1))*V(1+1)+Bus_Ymag(6+1)*cos(theta(0+1)-theta(2+1)-Bus_Yang(6+1))*V(2+1)); % Depth 1, Step 41
  mismatch(41+1)=P_load(1+1)+V(1+1)*(Bus_Ymag(1+1)*cos(theta(1+1)-theta(0+1)-Bus_Yang(1+1))*V(0+1)+Bus_Ymag(4+1)*cos(Bus_Yang(4+1))*V(1+1)+Bus_Ymag(39+1)*cos(theta(1+1)-theta(11+1)-Bus_Yang(39+1))*V(11+1)); % Depth 1, Step 42
  mismatch(42+1)=P_load(2+1)+V(2+1)*(Bus_Ymag(2+1)*cos(-theta(2+1)+theta(0+1)+Bus_Yang(2+1))*V(0+1)+Bus_Ymag(7+1)*cos(Bus_Yang(7+1))*V(2+1)+Bus_Ymag(13+1)*cos(-theta(2+1)+theta(4+1)+Bus_Yang(13+1))*V(4+1)+Bus_Ymag(40+1)*cos(theta(2+1)-theta(11+1)-Bus_Yang(40+1))*V(11+1)); % Depth 1, Step 43
  mismatch(43+1)=P_load(3+1)+V(3+1)*(Bus_Ymag(10+1)*cos(Bus_Yang(10+1))*V(3+1)+Bus_Ymag(14+1)*cos(-theta(3+1)+theta(4+1)+Bus_Yang(14+1))*V(4+1)+Bus_Ymag(34+1)*cos(theta(3+1)-theta(10+1)-Bus_Yang(34+1))*V(10+1))-5.5555555555555553581825733999722*V(3+1)*E(0+1)*cos(theta(3+1)-delta_t1(0+1)+1.5707963267948965579989817342721); % Depth 1, Step 44
  mismatch(44+1)=P_load(4+1)+V(4+1)*(Bus_Ymag(8+1)*cos(theta(4+1)-theta(2+1)-Bus_Yang(8+1))*V(2+1)+Bus_Ymag(11+1)*cos(theta(4+1)-theta(3+1)-Bus_Yang(11+1))*V(3+1)+Bus_Ymag(15+1)*cos(Bus_Yang(15+1))*V(4+1)+Bus_Ymag(19+1)*cos(-theta(4+1)+theta(5+1)+Bus_Yang(19+1))*V(5+1)+Bus_Ymag(25+1)*cos(-theta(4+1)+theta(7+1)+Bus_Yang(25+1))*V(7+1)+Bus_Ymag(35+1)*cos(theta(4+1)-theta(10+1)-Bus_Yang(35+1))*V(10+1)); % Depth 1, Step 45
  mismatch(45+1)=P_load(5+1)+V(5+1)*(Bus_Ymag(16+1)*cos(theta(5+1)-theta(4+1)-Bus_Yang(16+1))*V(4+1)+Bus_Ymag(20+1)*cos(Bus_Yang(20+1))*V(5+1)+Bus_Ymag(22+1)*cos(-theta(5+1)+theta(6+1)+Bus_Yang(22+1))*V(6+1)); % Depth 1, Step 46
  mismatch(46+1)=P_load(6+1)+V(6+1)*(Bus_Ymag(21+1)*cos(theta(6+1)-theta(5+1)-Bus_Yang(21+1))*V(5+1)+Bus_Ymag(23+1)*cos(Bus_Yang(23+1))*V(6+1)+Bus_Ymag(41+1)*cos(theta(6+1)-theta(11+1)-Bus_Yang(41+1))*V(11+1)); % Depth 1, Step 47
  mismatch(47+1)=P_load(7+1)+V(7+1)*(Bus_Ymag(17+1)*cos(theta(7+1)-theta(4+1)-Bus_Yang(17+1))*V(4+1)+Bus_Ymag(26+1)*cos(Bus_Yang(26+1))*V(7+1)+Bus_Ymag(29+1)*cos(theta(7+1)-theta(8+1)-Bus_Yang(29+1))*V(8+1)+Bus_Ymag(113+1)*cos(theta(7+1)-theta(29+1)-Bus_Yang(113+1))*V(29+1)); % Depth 1, Step 48
  mismatch(48+1)=P_load(8+1)+V(8+1)*(Bus_Ymag(27+1)*cos(-theta(8+1)+theta(7+1)+Bus_Yang(27+1))*V(7+1)+Bus_Ymag(30+1)*cos(Bus_Yang(30+1))*V(8+1)+Bus_Ymag(32+1)*cos(-theta(8+1)+theta(9+1)+Bus_Yang(32+1))*V(9+1)); % Depth 1, Step 49
  mismatch(49+1)=P_load(9+1)+V(9+1)*(Bus_Ymag(31+1)*cos(theta(9+1)-theta(8+1)-Bus_Yang(31+1))*V(8+1)+Bus_Ymag(33+1)*cos(Bus_Yang(33+1))*V(9+1))-5.5555555555555553581825733999722*V(9+1)*E(1+1)*cos(theta(9+1)-delta_t1(1+1)+1.5707963267948965579989817342721); % Depth 1, Step 50
  mismatch(50+1)=P_load(10+1)+V(10+1)*(Bus_Ymag(12+1)*cos(-theta(10+1)+theta(3+1)+Bus_Yang(12+1))*V(3+1)+Bus_Ymag(18+1)*cos(-theta(10+1)+theta(4+1)+Bus_Yang(18+1))*V(4+1)+Bus_Ymag(36+1)*cos(Bus_Yang(36+1))*V(10+1)+Bus_Ymag(42+1)*cos(-theta(10+1)+theta(11+1)+Bus_Yang(42+1))*V(11+1)+Bus_Ymag(47+1)*cos(-theta(10+1)+theta(12+1)+Bus_Yang(47+1))*V(12+1)); % Depth 1, Step 51
  mismatch(51+1)=P_load(11+1)+V(11+1)*(Bus_Ymag(5+1)*cos(-theta(11+1)+theta(1+1)+Bus_Yang(5+1))*V(1+1)+Bus_Ymag(9+1)*cos(-theta(11+1)+theta(2+1)+Bus_Yang(9+1))*V(2+1)+Bus_Ymag(24+1)*cos(-theta(11+1)+theta(6+1)+Bus_Yang(24+1))*V(6+1)+Bus_Ymag(37+1)*cos(theta(11+1)-theta(10+1)-Bus_Yang(37+1))*V(10+1)+Bus_Ymag(43+1)*cos(Bus_Yang(43+1))*V(11+1)+Bus_Ymag(50+1)*cos(-theta(11+1)+theta(13+1)+Bus_Yang(50+1))*V(13+1)+Bus_Ymag(59+1)*cos(-theta(11+1)+theta(15+1)+Bus_Yang(59+1))*V(15+1)+Bus_Ymag(471+1)*cos(theta(11+1)-theta(116+1)-Bus_Yang(471+1))*V(116+1))-5.5555555555555553581825733999722*V(11+1)*E(2+1)*cos(-theta(11+1)+delta_t1(2+1)-1.5707963267948965579989817342721); % Depth 1, Step 52
  mismatch(52+1)=P_load(12+1)+V(12+1)*(Bus_Ymag(38+1)*cos(theta(12+1)-theta(10+1)-Bus_Yang(38+1))*V(10+1)+Bus_Ymag(48+1)*cos(Bus_Yang(48+1))*V(12+1)+Bus_Ymag(53+1)*cos(-theta(12+1)+theta(14+1)+Bus_Yang(53+1))*V(14+1)); % Depth 1, Step 53
  mismatch(53+1)=P_load(13+1)+V(13+1)*(Bus_Ymag(44+1)*cos(theta(13+1)-theta(11+1)-Bus_Yang(44+1))*V(11+1)+Bus_Ymag(51+1)*cos(Bus_Yang(51+1))*V(13+1)+Bus_Ymag(54+1)*cos(-theta(13+1)+theta(14+1)+Bus_Yang(54+1))*V(14+1)); % Depth 1, Step 54
  mismatch(54+1)=P_load(14+1)+V(14+1)*(Bus_Ymag(49+1)*cos(theta(14+1)-theta(12+1)-Bus_Yang(49+1))*V(12+1)+Bus_Ymag(52+1)*cos(theta(14+1)-theta(13+1)-Bus_Yang(52+1))*V(13+1)+Bus_Ymag(55+1)*cos(Bus_Yang(55+1))*V(14+1)+Bus_Ymag(62+1)*cos(theta(14+1)-theta(16+1)-Bus_Yang(62+1))*V(16+1)+Bus_Ymag(72+1)*cos(theta(14+1)-theta(18+1)-Bus_Yang(72+1))*V(18+1)+Bus_Ymag(128+1)*cos(theta(14+1)-theta(32+1)-Bus_Yang(128+1))*V(32+1)); % Depth 1, Step 55
  mismatch(55+1)=P_load(15+1)+V(15+1)*(Bus_Ymag(45+1)*cos(theta(15+1)-theta(11+1)-Bus_Yang(45+1))*V(11+1)+Bus_Ymag(60+1)*cos(Bus_Yang(60+1))*V(15+1)+Bus_Ymag(63+1)*cos(theta(15+1)-theta(16+1)-Bus_Yang(63+1))*V(16+1)); % Depth 1, Step 56
  mismatch(56+1)=P_load(16+1)+V(16+1)*(Bus_Ymag(56+1)*cos(-theta(16+1)+theta(14+1)+Bus_Yang(56+1))*V(14+1)+Bus_Ymag(61+1)*cos(-theta(16+1)+theta(15+1)+Bus_Yang(61+1))*V(15+1)+Bus_Ymag(64+1)*cos(Bus_Yang(64+1))*V(16+1)+Bus_Ymag(69+1)*cos(-theta(16+1)+theta(17+1)+Bus_Yang(69+1))*V(17+1)+Bus_Ymag(114+1)*cos(theta(16+1)-theta(29+1)-Bus_Yang(114+1))*V(29+1)+Bus_Ymag(118+1)*cos(theta(16+1)-theta(30+1)-Bus_Yang(118+1))*V(30+1)+Bus_Ymag(460+1)*cos(theta(16+1)-theta(112+1)-Bus_Yang(460+1))*V(112+1)); % Depth 1, Step 57
  mismatch(57+1)=P_load(17+1)+V(17+1)*(Bus_Ymag(65+1)*cos(theta(17+1)-theta(16+1)-Bus_Yang(65+1))*V(16+1)+Bus_Ymag(70+1)*cos(Bus_Yang(70+1))*V(17+1)+Bus_Ymag(73+1)*cos(-theta(17+1)+theta(18+1)+Bus_Yang(73+1))*V(18+1)); % Depth 1, Step 58
  mismatch(58+1)=P_load(18+1)+V(18+1)*(Bus_Ymag(57+1)*cos(-theta(18+1)+theta(14+1)+Bus_Yang(57+1))*V(14+1)+Bus_Ymag(71+1)*cos(theta(18+1)-theta(17+1)-Bus_Yang(71+1))*V(17+1)+Bus_Ymag(74+1)*cos(Bus_Yang(74+1))*V(18+1)+Bus_Ymag(77+1)*cos(-theta(18+1)+theta(19+1)+Bus_Yang(77+1))*V(19+1)+Bus_Ymag(131+1)*cos(theta(18+1)-theta(33+1)-Bus_Yang(131+1))*V(33+1)); % Depth 1, Step 59
  mismatch(59+1)=P_load(19+1)+V(19+1)*(Bus_Ymag(75+1)*cos(theta(19+1)-theta(18+1)-Bus_Yang(75+1))*V(18+1)+Bus_Ymag(78+1)*cos(Bus_Yang(78+1))*V(19+1)+Bus_Ymag(80+1)*cos(-theta(19+1)+theta(20+1)+Bus_Yang(80+1))*V(20+1)); % Depth 1, Step 60
  mismatch(60+1)=P_load(20+1)+V(20+1)*(Bus_Ymag(79+1)*cos(theta(20+1)-theta(19+1)-Bus_Yang(79+1))*V(19+1)+Bus_Ymag(81+1)*cos(Bus_Yang(81+1))*V(20+1)+Bus_Ymag(83+1)*cos(-theta(20+1)+theta(21+1)+Bus_Yang(83+1))*V(21+1)); % Depth 1, Step 61
  mismatch(61+1)=P_load(21+1)+V(21+1)*(Bus_Ymag(82+1)*cos(theta(21+1)-theta(20+1)-Bus_Yang(82+1))*V(20+1)+Bus_Ymag(84+1)*cos(Bus_Yang(84+1))*V(21+1)+Bus_Ymag(86+1)*cos(-theta(21+1)+theta(22+1)+Bus_Yang(86+1))*V(22+1)); % Depth 1, Step 62
  mismatch(62+1)=P_load(22+1)+V(22+1)*(Bus_Ymag(85+1)*cos(theta(22+1)-theta(21+1)-Bus_Yang(85+1))*V(21+1)+Bus_Ymag(87+1)*cos(Bus_Yang(87+1))*V(22+1)+Bus_Ymag(91+1)*cos(-theta(22+1)+theta(23+1)+Bus_Yang(91+1))*V(23+1)+Bus_Ymag(95+1)*cos(-theta(22+1)+theta(24+1)+Bus_Yang(95+1))*V(24+1)+Bus_Ymag(122+1)*cos(theta(22+1)-theta(31+1)-Bus_Yang(122+1))*V(31+1)); % Depth 1, Step 63
  mismatch(63+1)=P_load(23+1)+V(23+1)*(Bus_Ymag(88+1)*cos(theta(23+1)-theta(22+1)-Bus_Yang(88+1))*V(22+1)+Bus_Ymag(92+1)*cos(Bus_Yang(92+1))*V(23+1)+Bus_Ymag(287+1)*cos(-theta(23+1)+theta(69+1)+Bus_Yang(287+1))*V(69+1)+Bus_Ymag(297+1)*cos(-theta(23+1)+theta(71+1)+Bus_Yang(297+1))*V(71+1)); % Depth 1, Step 64
  mismatch(64+1)=P_load(24+1)+V(24+1)*(Bus_Ymag(89+1)*cos(theta(24+1)-theta(22+1)-Bus_Yang(89+1))*V(22+1)+Bus_Ymag(96+1)*cos(Bus_Yang(96+1))*V(24+1)+Bus_Ymag(99+1)*cos(-theta(24+1)+theta(25+1)+Bus_Yang(99+1))*V(25+1)+Bus_Ymag(102+1)*cos(-theta(24+1)+theta(26+1)+Bus_Yang(102+1))*V(26+1))-5.5555555555555553581825733999722*V(24+1)*E(3+1)*cos(-theta(24+1)+delta_t1(3+1)-1.5707963267948965579989817342721); % Depth 1, Step 65
  mismatch(65+1)=P_load(25+1)+V(25+1)*(Bus_Ymag(97+1)*cos(theta(25+1)-theta(24+1)-Bus_Yang(97+1))*V(24+1)+Bus_Ymag(100+1)*cos(Bus_Yang(100+1))*V(25+1)+Bus_Ymag(115+1)*cos(theta(25+1)-theta(29+1)-Bus_Yang(115+1))*V(29+1))-5.5555555555555553581825733999722*V(25+1)*E(4+1)*cos(-theta(25+1)+delta_t1(4+1)-1.5707963267948965579989817342721); % Depth 1, Step 66
  mismatch(66+1)=P_load(26+1)+V(26+1)*(Bus_Ymag(98+1)*cos(theta(26+1)-theta(24+1)-Bus_Yang(98+1))*V(24+1)+Bus_Ymag(103+1)*cos(Bus_Yang(103+1))*V(26+1)+Bus_Ymag(107+1)*cos(-theta(26+1)+theta(27+1)+Bus_Yang(107+1))*V(27+1)+Bus_Ymag(123+1)*cos(theta(26+1)-theta(31+1)-Bus_Yang(123+1))*V(31+1)+Bus_Ymag(466+1)*cos(theta(26+1)-theta(114+1)-Bus_Yang(466+1))*V(114+1))-5.5555555555555553581825733999722*V(26+1)*E(5+1)*cos(-theta(26+1)+delta_t1(5+1)-1.5707963267948965579989817342721); % Depth 1, Step 67
  mismatch(67+1)=P_load(27+1)+V(27+1)*(Bus_Ymag(104+1)*cos(theta(27+1)-theta(26+1)-Bus_Yang(104+1))*V(26+1)+Bus_Ymag(108+1)*cos(Bus_Yang(108+1))*V(27+1)+Bus_Ymag(110+1)*cos(-theta(27+1)+theta(28+1)+Bus_Yang(110+1))*V(28+1)); % Depth 1, Step 68
  mismatch(68+1)=P_load(28+1)+V(28+1)*(Bus_Ymag(109+1)*cos(theta(28+1)-theta(27+1)-Bus_Yang(109+1))*V(27+1)+Bus_Ymag(111+1)*cos(Bus_Yang(111+1))*V(28+1)+Bus_Ymag(119+1)*cos(theta(28+1)-theta(30+1)-Bus_Yang(119+1))*V(30+1)); % Depth 1, Step 69
  mismatch(69+1)=P_load(29+1)+V(29+1)*(Bus_Ymag(28+1)*cos(-theta(29+1)+theta(7+1)+Bus_Yang(28+1))*V(7+1)+Bus_Ymag(66+1)*cos(-theta(29+1)+theta(16+1)+Bus_Yang(66+1))*V(16+1)+Bus_Ymag(101+1)*cos(-theta(29+1)+theta(25+1)+Bus_Yang(101+1))*V(25+1)+Bus_Ymag(116+1)*cos(Bus_Yang(116+1))*V(29+1)+Bus_Ymag(149+1)*cos(-theta(29+1)+theta(37+1)+Bus_Yang(149+1))*V(37+1)); % Depth 1, Step 70
  mismatch(70+1)=P_load(30+1)+V(30+1)*(Bus_Ymag(67+1)*cos(-theta(30+1)+theta(16+1)+Bus_Yang(67+1))*V(16+1)+Bus_Ymag(112+1)*cos(-theta(30+1)+theta(28+1)+Bus_Yang(112+1))*V(28+1)+Bus_Ymag(120+1)*cos(Bus_Yang(120+1))*V(30+1)+Bus_Ymag(124+1)*cos(-theta(30+1)+theta(31+1)+Bus_Yang(124+1))*V(31+1))-5.5555555555555553581825733999722*V(30+1)*E(6+1)*cos(-theta(30+1)+delta_t1(6+1)-1.5707963267948965579989817342721); % Depth 1, Step 71
  mismatch(71+1)=P_load(31+1)+V(31+1)*(Bus_Ymag(90+1)*cos(-theta(31+1)+theta(22+1)+Bus_Yang(90+1))*V(22+1)+Bus_Ymag(105+1)*cos(-theta(31+1)+theta(26+1)+Bus_Yang(105+1))*V(26+1)+Bus_Ymag(121+1)*cos(theta(31+1)-theta(30+1)-Bus_Yang(121+1))*V(30+1)+Bus_Ymag(125+1)*cos(Bus_Yang(125+1))*V(31+1)+Bus_Ymag(461+1)*cos(theta(31+1)-theta(112+1)-Bus_Yang(461+1))*V(112+1)+Bus_Ymag(463+1)*cos(theta(31+1)-theta(113+1)-Bus_Yang(463+1))*V(113+1)); % Depth 1, Step 72
  mismatch(72+1)=P_load(32+1)+V(32+1)*(Bus_Ymag(58+1)*cos(-theta(32+1)+theta(14+1)+Bus_Yang(58+1))*V(14+1)+Bus_Ymag(129+1)*cos(Bus_Yang(129+1))*V(32+1)+Bus_Ymag(142+1)*cos(-theta(32+1)+theta(36+1)+Bus_Yang(142+1))*V(36+1)); % Depth 1, Step 73
  mismatch(73+1)=P_load(33+1)+V(33+1)*(Bus_Ymag(76+1)*cos(-theta(33+1)+theta(18+1)+Bus_Yang(76+1))*V(18+1)+Bus_Ymag(132+1)*cos(Bus_Yang(132+1))*V(33+1)+Bus_Ymag(139+1)*cos(-theta(33+1)+theta(35+1)+Bus_Yang(139+1))*V(35+1)+Bus_Ymag(143+1)*cos(-theta(33+1)+theta(36+1)+Bus_Yang(143+1))*V(36+1)+Bus_Ymag(168+1)*cos(theta(33+1)-theta(42+1)-Bus_Yang(168+1))*V(42+1)); % Depth 1, Step 74
  mismatch(74+1)=P_load(34+1)+V(34+1)*(Bus_Ymag(136+1)*cos(Bus_Yang(136+1))*V(34+1)+Bus_Ymag(140+1)*cos(-theta(34+1)+theta(35+1)+Bus_Yang(140+1))*V(35+1)+Bus_Ymag(144+1)*cos(-theta(34+1)+theta(36+1)+Bus_Yang(144+1))*V(36+1)); % Depth 1, Step 75
  mismatch(75+1)=P_load(35+1)+V(35+1)*(Bus_Ymag(133+1)*cos(theta(35+1)-theta(33+1)-Bus_Yang(133+1))*V(33+1)+Bus_Ymag(137+1)*cos(theta(35+1)-theta(34+1)-Bus_Yang(137+1))*V(34+1)+Bus_Ymag(141+1)*cos(Bus_Yang(141+1))*V(35+1)); % Depth 1, Step 76
  mismatch(76+1)=P_load(36+1)+V(36+1)*(Bus_Ymag(130+1)*cos(theta(36+1)-theta(32+1)-Bus_Yang(130+1))*V(32+1)+Bus_Ymag(134+1)*cos(theta(36+1)-theta(33+1)-Bus_Yang(134+1))*V(33+1)+Bus_Ymag(138+1)*cos(theta(36+1)-theta(34+1)-Bus_Yang(138+1))*V(34+1)+Bus_Ymag(145+1)*cos(Bus_Yang(145+1))*V(36+1)+Bus_Ymag(150+1)*cos(-theta(36+1)+theta(37+1)+Bus_Yang(150+1))*V(37+1)+Bus_Ymag(153+1)*cos(theta(36+1)-theta(38+1)-Bus_Yang(153+1))*V(38+1)+Bus_Ymag(156+1)*cos(theta(36+1)-theta(39+1)-Bus_Yang(156+1))*V(39+1)); % Depth 1, Step 77
  mismatch(77+1)=P_load(37+1)+V(37+1)*(Bus_Ymag(117+1)*cos(theta(37+1)-theta(29+1)-Bus_Yang(117+1))*V(29+1)+Bus_Ymag(146+1)*cos(theta(37+1)-theta(36+1)-Bus_Yang(146+1))*V(36+1)+Bus_Ymag(151+1)*cos(Bus_Yang(151+1))*V(37+1)+Bus_Ymag(262+1)*cos(theta(37+1)-theta(64+1)-Bus_Yang(262+1))*V(64+1)); % Depth 1, Step 78
  mismatch(78+1)=P_load(38+1)+V(38+1)*(Bus_Ymag(147+1)*cos(-theta(38+1)+theta(36+1)+Bus_Yang(147+1))*V(36+1)+Bus_Ymag(154+1)*cos(Bus_Yang(154+1))*V(38+1)+Bus_Ymag(157+1)*cos(-theta(38+1)+theta(39+1)+Bus_Yang(157+1))*V(39+1)); % Depth 1, Step 79
  mismatch(79+1)=P_load(39+1)+V(39+1)*(Bus_Ymag(148+1)*cos(-theta(39+1)+theta(36+1)+Bus_Yang(148+1))*V(36+1)+Bus_Ymag(155+1)*cos(theta(39+1)-theta(38+1)-Bus_Yang(155+1))*V(38+1)+Bus_Ymag(158+1)*cos(Bus_Yang(158+1))*V(39+1)+Bus_Ymag(161+1)*cos(-theta(39+1)+theta(40+1)+Bus_Yang(161+1))*V(40+1)+Bus_Ymag(164+1)*cos(-theta(39+1)+theta(41+1)+Bus_Yang(164+1))*V(41+1))-5.5555555555555553581825733999722*V(39+1)*E(7+1)*cos(theta(39+1)-delta_t1(7+1)+1.5707963267948965579989817342721); % Depth 1, Step 80
  mismatch(80+1)=P_load(40+1)+V(40+1)*(Bus_Ymag(159+1)*cos(theta(40+1)-theta(39+1)-Bus_Yang(159+1))*V(39+1)+Bus_Ymag(162+1)*cos(Bus_Yang(162+1))*V(40+1)+Bus_Ymag(165+1)*cos(-theta(40+1)+theta(41+1)+Bus_Yang(165+1))*V(41+1)); % Depth 1, Step 81
  mismatch(81+1)=P_load(41+1)+V(41+1)*(Bus_Ymag(160+1)*cos(theta(41+1)-theta(39+1)-Bus_Yang(160+1))*V(39+1)+Bus_Ymag(163+1)*cos(theta(41+1)-theta(40+1)-Bus_Yang(163+1))*V(40+1)+Bus_Ymag(166+1)*cos(Bus_Yang(166+1))*V(41+1)+Bus_Ymag(189+1)*cos(theta(41+1)-theta(48+1)-Bus_Yang(189+1))*V(48+1))-5.5555555555555553581825733999722*V(41+1)*E(8+1)*cos(theta(41+1)-delta_t1(8+1)+1.5707963267948965579989817342721); % Depth 1, Step 82
  mismatch(82+1)=P_load(42+1)+V(42+1)*(Bus_Ymag(135+1)*cos(-theta(42+1)+theta(33+1)+Bus_Yang(135+1))*V(33+1)+Bus_Ymag(169+1)*cos(Bus_Yang(169+1))*V(42+1)+Bus_Ymag(171+1)*cos(-theta(42+1)+theta(43+1)+Bus_Yang(171+1))*V(43+1)); % Depth 1, Step 83
  mismatch(83+1)=P_load(43+1)+V(43+1)*(Bus_Ymag(170+1)*cos(theta(43+1)-theta(42+1)-Bus_Yang(170+1))*V(42+1)+Bus_Ymag(172+1)*cos(Bus_Yang(172+1))*V(43+1)+Bus_Ymag(476+1)*cos(-theta(43+1)+theta(118+1)+Bus_Yang(476+1))*V(118+1)); % Depth 1, Step 84
  mismatch(84+1)=P_load(44+1)+V(44+1)*(Bus_Ymag(174+1)*cos(Bus_Yang(174+1))*V(44+1)+Bus_Ymag(178+1)*cos(-theta(44+1)+theta(45+1)+Bus_Yang(178+1))*V(45+1)+Bus_Ymag(190+1)*cos(theta(44+1)-theta(48+1)-Bus_Yang(190+1))*V(48+1)+Bus_Ymag(477+1)*cos(-theta(44+1)+theta(118+1)+Bus_Yang(477+1))*V(118+1)); % Depth 1, Step 85
  mismatch(85+1)=P_load(45+1)+V(45+1)*(Bus_Ymag(175+1)*cos(theta(45+1)-theta(44+1)-Bus_Yang(175+1))*V(44+1)+Bus_Ymag(179+1)*cos(Bus_Yang(179+1))*V(45+1)+Bus_Ymag(182+1)*cos(-theta(45+1)+theta(46+1)+Bus_Yang(182+1))*V(46+1)+Bus_Ymag(186+1)*cos(-theta(45+1)+theta(47+1)+Bus_Yang(186+1))*V(47+1)); % Depth 1, Step 86
  mismatch(86+1)=P_load(46+1)+V(46+1)*(Bus_Ymag(180+1)*cos(theta(46+1)-theta(45+1)-Bus_Yang(180+1))*V(45+1)+Bus_Ymag(183+1)*cos(Bus_Yang(183+1))*V(46+1)+Bus_Ymag(191+1)*cos(theta(46+1)-theta(48+1)-Bus_Yang(191+1))*V(48+1)+Bus_Ymag(280+1)*cos(-theta(46+1)+theta(68+1)+Bus_Yang(280+1))*V(68+1)); % Depth 1, Step 87
  mismatch(87+1)=P_load(47+1)+V(47+1)*(Bus_Ymag(181+1)*cos(theta(47+1)-theta(45+1)-Bus_Yang(181+1))*V(45+1)+Bus_Ymag(187+1)*cos(Bus_Yang(187+1))*V(47+1)+Bus_Ymag(192+1)*cos(theta(47+1)-theta(48+1)-Bus_Yang(192+1))*V(48+1)); % Depth 1, Step 88
  mismatch(88+1)=P_load(48+1)+V(48+1)*(Bus_Ymag(167+1)*cos(-theta(48+1)+theta(41+1)+Bus_Yang(167+1))*V(41+1)+Bus_Ymag(176+1)*cos(-theta(48+1)+theta(44+1)+Bus_Yang(176+1))*V(44+1)+Bus_Ymag(184+1)*cos(-theta(48+1)+theta(46+1)+Bus_Yang(184+1))*V(46+1)+Bus_Ymag(188+1)*cos(-theta(48+1)+theta(47+1)+Bus_Yang(188+1))*V(47+1)+Bus_Ymag(193+1)*cos(Bus_Yang(193+1))*V(48+1)+Bus_Ymag(199+1)*cos(-theta(48+1)+theta(49+1)+Bus_Yang(199+1))*V(49+1)+Bus_Ymag(202+1)*cos(-theta(48+1)+theta(50+1)+Bus_Yang(202+1))*V(50+1)+Bus_Ymag(212+1)*cos(theta(48+1)-theta(53+1)-Bus_Yang(212+1))*V(53+1)+Bus_Ymag(267+1)*cos(-theta(48+1)+theta(65+1)+Bus_Yang(267+1))*V(65+1)+Bus_Ymag(281+1)*cos(-theta(48+1)+theta(68+1)+Bus_Yang(281+1))*V(68+1))-5.5555555555555553581825733999722*V(48+1)*E(9+1)*cos(theta(48+1)-delta_t1(9+1)+1.5707963267948965579989817342721); % Depth 1, Step 89
  mismatch(89+1)=P_load(49+1)+V(49+1)*(Bus_Ymag(194+1)*cos(theta(49+1)-theta(48+1)-Bus_Yang(194+1))*V(48+1)+Bus_Ymag(200+1)*cos(Bus_Yang(200+1))*V(49+1)+Bus_Ymag(228+1)*cos(theta(49+1)-theta(56+1)-Bus_Yang(228+1))*V(56+1)); % Depth 1, Step 90
  mismatch(90+1)=P_load(50+1)+V(50+1)*(Bus_Ymag(195+1)*cos(theta(50+1)-theta(48+1)-Bus_Yang(195+1))*V(48+1)+Bus_Ymag(203+1)*cos(Bus_Yang(203+1))*V(50+1)+Bus_Ymag(206+1)*cos(-theta(50+1)+theta(51+1)+Bus_Yang(206+1))*V(51+1)+Bus_Ymag(231+1)*cos(theta(50+1)-theta(57+1)-Bus_Yang(231+1))*V(57+1)); % Depth 1, Step 91
  mismatch(91+1)=P_load(51+1)+V(51+1)*(Bus_Ymag(204+1)*cos(theta(51+1)-theta(50+1)-Bus_Yang(204+1))*V(50+1)+Bus_Ymag(207+1)*cos(Bus_Yang(207+1))*V(51+1)+Bus_Ymag(209+1)*cos(-theta(51+1)+theta(52+1)+Bus_Yang(209+1))*V(52+1)); % Depth 1, Step 92
  mismatch(92+1)=P_load(52+1)+V(52+1)*(Bus_Ymag(208+1)*cos(theta(52+1)-theta(51+1)-Bus_Yang(208+1))*V(51+1)+Bus_Ymag(210+1)*cos(Bus_Yang(210+1))*V(52+1)+Bus_Ymag(213+1)*cos(theta(52+1)-theta(53+1)-Bus_Yang(213+1))*V(53+1)); % Depth 1, Step 93
  mismatch(93+1)=P_load(53+1)+V(53+1)*(Bus_Ymag(196+1)*cos(-theta(53+1)+theta(48+1)+Bus_Yang(196+1))*V(48+1)+Bus_Ymag(211+1)*cos(-theta(53+1)+theta(52+1)+Bus_Yang(211+1))*V(52+1)+Bus_Ymag(214+1)*cos(Bus_Yang(214+1))*V(53+1)+Bus_Ymag(218+1)*cos(-theta(53+1)+theta(54+1)+Bus_Yang(218+1))*V(54+1)+Bus_Ymag(222+1)*cos(-theta(53+1)+theta(55+1)+Bus_Yang(222+1))*V(55+1)+Bus_Ymag(234+1)*cos(-theta(53+1)+theta(58+1)+Bus_Yang(234+1))*V(58+1)); % Depth 1, Step 94
  mismatch(94+1)=P_load(54+1)+V(54+1)*(Bus_Ymag(215+1)*cos(theta(54+1)-theta(53+1)-Bus_Yang(215+1))*V(53+1)+Bus_Ymag(219+1)*cos(Bus_Yang(219+1))*V(54+1)+Bus_Ymag(223+1)*cos(-theta(54+1)+theta(55+1)+Bus_Yang(223+1))*V(55+1)+Bus_Ymag(235+1)*cos(-theta(54+1)+theta(58+1)+Bus_Yang(235+1))*V(58+1)); % Depth 1, Step 95
  mismatch(95+1)=P_load(55+1)+V(55+1)*(Bus_Ymag(216+1)*cos(theta(55+1)-theta(53+1)-Bus_Yang(216+1))*V(53+1)+Bus_Ymag(220+1)*cos(theta(55+1)-theta(54+1)-Bus_Yang(220+1))*V(54+1)+Bus_Ymag(224+1)*cos(Bus_Yang(224+1))*V(55+1)+Bus_Ymag(229+1)*cos(-theta(55+1)+theta(56+1)+Bus_Yang(229+1))*V(56+1)+Bus_Ymag(232+1)*cos(-theta(55+1)+theta(57+1)+Bus_Yang(232+1))*V(57+1)+Bus_Ymag(236+1)*cos(-theta(55+1)+theta(58+1)+Bus_Yang(236+1))*V(58+1)); % Depth 1, Step 96
  mismatch(96+1)=P_load(56+1)+V(56+1)*(Bus_Ymag(201+1)*cos(-theta(56+1)+theta(49+1)+Bus_Yang(201+1))*V(49+1)+Bus_Ymag(225+1)*cos(theta(56+1)-theta(55+1)-Bus_Yang(225+1))*V(55+1)+Bus_Ymag(230+1)*cos(Bus_Yang(230+1))*V(56+1)); % Depth 1, Step 97
  mismatch(97+1)=P_load(57+1)+V(57+1)*(Bus_Ymag(205+1)*cos(-theta(57+1)+theta(50+1)+Bus_Yang(205+1))*V(50+1)+Bus_Ymag(226+1)*cos(theta(57+1)-theta(55+1)-Bus_Yang(226+1))*V(55+1)+Bus_Ymag(233+1)*cos(Bus_Yang(233+1))*V(57+1)); % Depth 1, Step 98
  mismatch(98+1)=P_load(58+1)+V(58+1)*(Bus_Ymag(217+1)*cos(theta(58+1)-theta(53+1)-Bus_Yang(217+1))*V(53+1)+Bus_Ymag(221+1)*cos(theta(58+1)-theta(54+1)-Bus_Yang(221+1))*V(54+1)+Bus_Ymag(227+1)*cos(theta(58+1)-theta(55+1)-Bus_Yang(227+1))*V(55+1)+Bus_Ymag(237+1)*cos(Bus_Yang(237+1))*V(58+1)+Bus_Ymag(241+1)*cos(theta(58+1)-theta(59+1)-Bus_Yang(241+1))*V(59+1)+Bus_Ymag(245+1)*cos(theta(58+1)-theta(60+1)-Bus_Yang(245+1))*V(60+1)+Bus_Ymag(255+1)*cos(-theta(58+1)+theta(62+1)+Bus_Yang(255+1))*V(62+1))-5.5555555555555553581825733999722*V(58+1)*E(10+1)*cos(theta(58+1)-delta_t1(10+1)+1.5707963267948965579989817342721); % Depth 1, Step 99
  mismatch(99+1)=P_load(59+1)+V(59+1)*(Bus_Ymag(238+1)*cos(-theta(59+1)+theta(58+1)+Bus_Yang(238+1))*V(58+1)+Bus_Ymag(242+1)*cos(Bus_Yang(242+1))*V(59+1)+Bus_Ymag(246+1)*cos(theta(59+1)-theta(60+1)-Bus_Yang(246+1))*V(60+1)+Bus_Ymag(250+1)*cos(theta(59+1)-theta(61+1)-Bus_Yang(250+1))*V(61+1)); % Depth 1, Step 100
  mismatch(100+1)=P_load(60+1)+V(60+1)*(Bus_Ymag(239+1)*cos(-theta(60+1)+theta(58+1)+Bus_Yang(239+1))*V(58+1)+Bus_Ymag(243+1)*cos(-theta(60+1)+theta(59+1)+Bus_Yang(243+1))*V(59+1)+Bus_Ymag(247+1)*cos(Bus_Yang(247+1))*V(60+1)+Bus_Ymag(251+1)*cos(theta(60+1)-theta(61+1)-Bus_Yang(251+1))*V(61+1)+Bus_Ymag(258+1)*cos(theta(60+1)-theta(63+1)-Bus_Yang(258+1))*V(63+1))-5.5555555555555553581825733999722*V(60+1)*E(11+1)*cos(theta(60+1)-delta_t1(11+1)+1.5707963267948965579989817342721); % Depth 1, Step 101
  mismatch(101+1)=P_load(61+1)+V(61+1)*(Bus_Ymag(244+1)*cos(-theta(61+1)+theta(59+1)+Bus_Yang(244+1))*V(59+1)+Bus_Ymag(248+1)*cos(-theta(61+1)+theta(60+1)+Bus_Yang(248+1))*V(60+1)+Bus_Ymag(252+1)*cos(Bus_Yang(252+1))*V(61+1)+Bus_Ymag(268+1)*cos(-theta(61+1)+theta(65+1)+Bus_Yang(268+1))*V(65+1)+Bus_Ymag(272+1)*cos(-theta(61+1)+theta(66+1)+Bus_Yang(272+1))*V(66+1)); % Depth 1, Step 102
  mismatch(102+1)=P_load(62+1)+V(62+1)*(Bus_Ymag(240+1)*cos(theta(62+1)-theta(58+1)-Bus_Yang(240+1))*V(58+1)+Bus_Ymag(256+1)*cos(Bus_Yang(256+1))*V(62+1)+Bus_Ymag(259+1)*cos(theta(62+1)-theta(63+1)-Bus_Yang(259+1))*V(63+1)); % Depth 1, Step 103
  mismatch(103+1)=P_load(63+1)+V(63+1)*(Bus_Ymag(249+1)*cos(-theta(63+1)+theta(60+1)+Bus_Yang(249+1))*V(60+1)+Bus_Ymag(257+1)*cos(-theta(63+1)+theta(62+1)+Bus_Yang(257+1))*V(62+1)+Bus_Ymag(260+1)*cos(Bus_Yang(260+1))*V(63+1)+Bus_Ymag(263+1)*cos(-theta(63+1)+theta(64+1)+Bus_Yang(263+1))*V(64+1)); % Depth 1, Step 104
  mismatch(104+1)=P_load(64+1)+V(64+1)*(Bus_Ymag(152+1)*cos(-theta(64+1)+theta(37+1)+Bus_Yang(152+1))*V(37+1)+Bus_Ymag(261+1)*cos(-theta(64+1)+theta(63+1)+Bus_Yang(261+1))*V(63+1)+Bus_Ymag(264+1)*cos(Bus_Yang(264+1))*V(64+1)+Bus_Ymag(269+1)*cos(-theta(64+1)+theta(65+1)+Bus_Yang(269+1))*V(65+1)+Bus_Ymag(275+1)*cos(-theta(64+1)+theta(67+1)+Bus_Yang(275+1))*V(67+1))-5.5555555555555553581825733999722*V(64+1)*E(12+1)*cos(-theta(64+1)+delta_t1(12+1)-1.5707963267948965579989817342721); % Depth 1, Step 105
  mismatch(105+1)=P_load(65+1)+V(65+1)*(Bus_Ymag(197+1)*cos(theta(65+1)-theta(48+1)-Bus_Yang(197+1))*V(48+1)+Bus_Ymag(253+1)*cos(theta(65+1)-theta(61+1)-Bus_Yang(253+1))*V(61+1)+Bus_Ymag(265+1)*cos(theta(65+1)-theta(64+1)-Bus_Yang(265+1))*V(64+1)+Bus_Ymag(270+1)*cos(Bus_Yang(270+1))*V(65+1)+Bus_Ymag(273+1)*cos(-theta(65+1)+theta(66+1)+Bus_Yang(273+1))*V(66+1))-5.5555555555555553581825733999722*V(65+1)*E(13+1)*cos(theta(65+1)-delta_t1(13+1)+1.5707963267948965579989817342721); % Depth 1, Step 106
  mismatch(106+1)=P_load(66+1)+V(66+1)*(Bus_Ymag(254+1)*cos(theta(66+1)-theta(61+1)-Bus_Yang(254+1))*V(61+1)+Bus_Ymag(271+1)*cos(theta(66+1)-theta(65+1)-Bus_Yang(271+1))*V(65+1)+Bus_Ymag(274+1)*cos(Bus_Yang(274+1))*V(66+1)); % Depth 1, Step 107
  mismatch(107+1)=P_load(67+1)+V(67+1)*(Bus_Ymag(266+1)*cos(theta(67+1)-theta(64+1)-Bus_Yang(266+1))*V(64+1)+Bus_Ymag(276+1)*cos(Bus_Yang(276+1))*V(67+1)+Bus_Ymag(282+1)*cos(-theta(67+1)+theta(68+1)+Bus_Yang(282+1))*V(68+1)+Bus_Ymag(335+1)*cos(theta(67+1)-theta(80+1)-Bus_Yang(335+1))*V(80+1)+Bus_Ymag(469+1)*cos(theta(67+1)-theta(115+1)-Bus_Yang(469+1))*V(115+1)); % Depth 1, Step 108
  mismatch(108+1)=P_load(68+1)+V(68+1)*(Bus_Ymag(185+1)*cos(theta(68+1)-theta(46+1)-Bus_Yang(185+1))*V(46+1)+Bus_Ymag(198+1)*cos(theta(68+1)-theta(48+1)-Bus_Yang(198+1))*V(48+1)+Bus_Ymag(277+1)*cos(theta(68+1)-theta(67+1)-Bus_Yang(277+1))*V(67+1)+Bus_Ymag(283+1)*cos(Bus_Yang(283+1))*V(68+1)+Bus_Ymag(288+1)*cos(-theta(68+1)+theta(69+1)+Bus_Yang(288+1))*V(69+1)+Bus_Ymag(305+1)*cos(theta(68+1)-theta(74+1)-Bus_Yang(305+1))*V(74+1)+Bus_Ymag(314+1)*cos(theta(68+1)-theta(76+1)-Bus_Yang(314+1))*V(76+1))-5.5555555555555553581825733999722*V(68+1)*E(14+1)*cos(theta(68+1)-delta_t1(14+1)+1.5707963267948965579989817342721); % Depth 1, Step 109
  mismatch(109+1)=P_load(69+1)+V(69+1)*(Bus_Ymag(93+1)*cos(theta(69+1)-theta(23+1)-Bus_Yang(93+1))*V(23+1)+Bus_Ymag(284+1)*cos(theta(69+1)-theta(68+1)-Bus_Yang(284+1))*V(68+1)+Bus_Ymag(289+1)*cos(Bus_Yang(289+1))*V(69+1)+Bus_Ymag(293+1)*cos(-theta(69+1)+theta(70+1)+Bus_Yang(293+1))*V(70+1)+Bus_Ymag(302+1)*cos(theta(69+1)-theta(73+1)-Bus_Yang(302+1))*V(73+1)+Bus_Ymag(306+1)*cos(theta(69+1)-theta(74+1)-Bus_Yang(306+1))*V(74+1)); % Depth 1, Step 110
  mismatch(110+1)=P_load(70+1)+V(70+1)*(Bus_Ymag(290+1)*cos(theta(70+1)-theta(69+1)-Bus_Yang(290+1))*V(69+1)+Bus_Ymag(294+1)*cos(Bus_Yang(294+1))*V(70+1)+Bus_Ymag(298+1)*cos(-theta(70+1)+theta(71+1)+Bus_Yang(298+1))*V(71+1)+Bus_Ymag(300+1)*cos(theta(70+1)-theta(72+1)-Bus_Yang(300+1))*V(72+1)); % Depth 1, Step 111
  mismatch(111+1)=P_load(71+1)+V(71+1)*(Bus_Ymag(94+1)*cos(theta(71+1)-theta(23+1)-Bus_Yang(94+1))*V(23+1)+Bus_Ymag(295+1)*cos(theta(71+1)-theta(70+1)-Bus_Yang(295+1))*V(70+1)+Bus_Ymag(299+1)*cos(Bus_Yang(299+1))*V(71+1)); % Depth 1, Step 112
  mismatch(112+1)=P_load(72+1)+V(72+1)*(Bus_Ymag(296+1)*cos(-theta(72+1)+theta(70+1)+Bus_Yang(296+1))*V(70+1)+Bus_Ymag(301+1)*cos(Bus_Yang(301+1))*V(72+1)); % Depth 1, Step 113
  mismatch(113+1)=P_load(73+1)+V(73+1)*(Bus_Ymag(291+1)*cos(-theta(73+1)+theta(69+1)+Bus_Yang(291+1))*V(69+1)+Bus_Ymag(303+1)*cos(Bus_Yang(303+1))*V(73+1)+Bus_Ymag(307+1)*cos(-theta(73+1)+theta(74+1)+Bus_Yang(307+1))*V(74+1)); % Depth 1, Step 114
  mismatch(114+1)=P_load(74+1)+V(74+1)*(Bus_Ymag(285+1)*cos(-theta(74+1)+theta(68+1)+Bus_Yang(285+1))*V(68+1)+Bus_Ymag(292+1)*cos(-theta(74+1)+theta(69+1)+Bus_Yang(292+1))*V(69+1)+Bus_Ymag(304+1)*cos(theta(74+1)-theta(73+1)-Bus_Yang(304+1))*V(73+1)+Bus_Ymag(308+1)*cos(Bus_Yang(308+1))*V(74+1)+Bus_Ymag(315+1)*cos(-theta(74+1)+theta(76+1)+Bus_Yang(315+1))*V(76+1)+Bus_Ymag(473+1)*cos(theta(74+1)-theta(117+1)-Bus_Yang(473+1))*V(117+1)); % Depth 1, Step 115
  mismatch(115+1)=P_load(75+1)+V(75+1)*(Bus_Ymag(311+1)*cos(Bus_Yang(311+1))*V(75+1)+Bus_Ymag(316+1)*cos(-theta(75+1)+theta(76+1)+Bus_Yang(316+1))*V(76+1)+Bus_Ymag(474+1)*cos(theta(75+1)-theta(117+1)-Bus_Yang(474+1))*V(117+1)); % Depth 1, Step 116
  mismatch(116+1)=P_load(76+1)+V(76+1)*(Bus_Ymag(286+1)*cos(-theta(76+1)+theta(68+1)+Bus_Yang(286+1))*V(68+1)+Bus_Ymag(309+1)*cos(theta(76+1)-theta(74+1)-Bus_Yang(309+1))*V(74+1)+Bus_Ymag(312+1)*cos(theta(76+1)-theta(75+1)-Bus_Yang(312+1))*V(75+1)+Bus_Ymag(317+1)*cos(Bus_Yang(317+1))*V(76+1)+Bus_Ymag(321+1)*cos(theta(76+1)-theta(77+1)-Bus_Yang(321+1))*V(77+1)+Bus_Ymag(327+1)*cos(-theta(76+1)+theta(79+1)+Bus_Yang(327+1))*V(79+1)+Bus_Ymag(338+1)*cos(theta(76+1)-theta(81+1)-Bus_Yang(338+1))*V(81+1)); % Depth 1, Step 117
  mismatch(117+1)=P_load(77+1)+V(77+1)*(Bus_Ymag(318+1)*cos(-theta(77+1)+theta(76+1)+Bus_Yang(318+1))*V(76+1)+Bus_Ymag(322+1)*cos(Bus_Yang(322+1))*V(77+1)+Bus_Ymag(324+1)*cos(-theta(77+1)+theta(78+1)+Bus_Yang(324+1))*V(78+1)); % Depth 1, Step 118
  mismatch(118+1)=P_load(78+1)+V(78+1)*(Bus_Ymag(323+1)*cos(theta(78+1)-theta(77+1)-Bus_Yang(323+1))*V(77+1)+Bus_Ymag(325+1)*cos(Bus_Yang(325+1))*V(78+1)+Bus_Ymag(328+1)*cos(-theta(78+1)+theta(79+1)+Bus_Yang(328+1))*V(79+1)); % Depth 1, Step 119
  mismatch(119+1)=P_load(79+1)+V(79+1)*(Bus_Ymag(319+1)*cos(-theta(79+1)+theta(76+1)+Bus_Yang(319+1))*V(76+1)+Bus_Ymag(326+1)*cos(-theta(79+1)+theta(78+1)+Bus_Yang(326+1))*V(78+1)+Bus_Ymag(329+1)*cos(Bus_Yang(329+1))*V(79+1)+Bus_Ymag(336+1)*cos(-theta(79+1)+theta(80+1)+Bus_Yang(336+1))*V(80+1)+Bus_Ymag(393+1)*cos(theta(79+1)-theta(95+1)-Bus_Yang(393+1))*V(95+1)+Bus_Ymag(399+1)*cos(theta(79+1)-theta(96+1)-Bus_Yang(399+1))*V(96+1)+Bus_Ymag(402+1)*cos(-theta(79+1)+theta(97+1)+Bus_Yang(402+1))*V(97+1)+Bus_Ymag(405+1)*cos(-theta(79+1)+theta(98+1)+Bus_Yang(405+1))*V(98+1))-5.5555555555555553581825733999722*V(79+1)*E(15+1)*cos(theta(79+1)-delta_t1(15+1)+1.5707963267948965579989817342721); % Depth 1, Step 120
  mismatch(120+1)=P_load(80+1)+V(80+1)*(Bus_Ymag(278+1)*cos(-theta(80+1)+theta(67+1)+Bus_Yang(278+1))*V(67+1)+Bus_Ymag(330+1)*cos(theta(80+1)-theta(79+1)-Bus_Yang(330+1))*V(79+1)+Bus_Ymag(337+1)*cos(Bus_Yang(337+1))*V(80+1)); % Depth 1, Step 121
  mismatch(121+1)=P_load(81+1)+V(81+1)*(Bus_Ymag(320+1)*cos(-theta(81+1)+theta(76+1)+Bus_Yang(320+1))*V(76+1)+Bus_Ymag(339+1)*cos(Bus_Yang(339+1))*V(81+1)+Bus_Ymag(342+1)*cos(-theta(81+1)+theta(82+1)+Bus_Yang(342+1))*V(82+1)+Bus_Ymag(394+1)*cos(theta(81+1)-theta(95+1)-Bus_Yang(394+1))*V(95+1)); % Depth 1, Step 122
  mismatch(122+1)=P_load(82+1)+V(82+1)*(Bus_Ymag(340+1)*cos(theta(82+1)-theta(81+1)-Bus_Yang(340+1))*V(81+1)+Bus_Ymag(343+1)*cos(Bus_Yang(343+1))*V(82+1)+Bus_Ymag(346+1)*cos(theta(82+1)-theta(83+1)-Bus_Yang(346+1))*V(83+1)+Bus_Ymag(349+1)*cos(theta(82+1)-theta(84+1)-Bus_Yang(349+1))*V(84+1)); % Depth 1, Step 123
  mismatch(123+1)=P_load(83+1)+V(83+1)*(Bus_Ymag(344+1)*cos(-theta(83+1)+theta(82+1)+Bus_Yang(344+1))*V(82+1)+Bus_Ymag(347+1)*cos(Bus_Yang(347+1))*V(83+1)+Bus_Ymag(350+1)*cos(-theta(83+1)+theta(84+1)+Bus_Yang(350+1))*V(84+1)); % Depth 1, Step 124
  mismatch(124+1)=P_load(84+1)+V(84+1)*(Bus_Ymag(345+1)*cos(-theta(84+1)+theta(82+1)+Bus_Yang(345+1))*V(82+1)+Bus_Ymag(348+1)*cos(theta(84+1)-theta(83+1)-Bus_Yang(348+1))*V(83+1)+Bus_Ymag(351+1)*cos(Bus_Yang(351+1))*V(84+1)+Bus_Ymag(355+1)*cos(-theta(84+1)+theta(85+1)+Bus_Yang(355+1))*V(85+1)+Bus_Ymag(360+1)*cos(-theta(84+1)+theta(87+1)+Bus_Yang(360+1))*V(87+1)+Bus_Ymag(363+1)*cos(-theta(84+1)+theta(88+1)+Bus_Yang(363+1))*V(88+1)); % Depth 1, Step 125
  mismatch(125+1)=P_load(85+1)+V(85+1)*(Bus_Ymag(352+1)*cos(theta(85+1)-theta(84+1)-Bus_Yang(352+1))*V(84+1)+Bus_Ymag(356+1)*cos(Bus_Yang(356+1))*V(85+1)+Bus_Ymag(358+1)*cos(-theta(85+1)+theta(86+1)+Bus_Yang(358+1))*V(86+1)); % Depth 1, Step 126
  mismatch(126+1)=P_load(86+1)+V(86+1)*(Bus_Ymag(357+1)*cos(theta(86+1)-theta(85+1)-Bus_Yang(357+1))*V(85+1)+Bus_Ymag(359+1)*cos(Bus_Yang(359+1))*V(86+1)); % Depth 1, Step 127
  mismatch(127+1)=P_load(87+1)+V(87+1)*(Bus_Ymag(353+1)*cos(theta(87+1)-theta(84+1)-Bus_Yang(353+1))*V(84+1)+Bus_Ymag(361+1)*cos(Bus_Yang(361+1))*V(87+1)+Bus_Ymag(364+1)*cos(-theta(87+1)+theta(88+1)+Bus_Yang(364+1))*V(88+1)); % Depth 1, Step 128
  mismatch(128+1)=P_load(88+1)+V(88+1)*(Bus_Ymag(354+1)*cos(theta(88+1)-theta(84+1)-Bus_Yang(354+1))*V(84+1)+Bus_Ymag(362+1)*cos(theta(88+1)-theta(87+1)-Bus_Yang(362+1))*V(87+1)+Bus_Ymag(365+1)*cos(Bus_Yang(365+1))*V(88+1)+Bus_Ymag(368+1)*cos(theta(88+1)-theta(89+1)-Bus_Yang(368+1))*V(89+1)+Bus_Ymag(374+1)*cos(theta(88+1)-theta(91+1)-Bus_Yang(374+1))*V(91+1))-5.5555555555555553581825733999722*V(88+1)*E(16+1)*cos(theta(88+1)-delta_t1(16+1)+1.5707963267948965579989817342721); % Depth 1, Step 129
  mismatch(129+1)=P_load(89+1)+V(89+1)*(Bus_Ymag(366+1)*cos(-theta(89+1)+theta(88+1)+Bus_Yang(366+1))*V(88+1)+Bus_Ymag(369+1)*cos(Bus_Yang(369+1))*V(89+1)+Bus_Ymag(371+1)*cos(-theta(89+1)+theta(90+1)+Bus_Yang(371+1))*V(90+1))-5.5555555555555553581825733999722*V(89+1)*E(17+1)*cos(theta(89+1)-delta_t1(17+1)+1.5707963267948965579989817342721); % Depth 1, Step 130
  mismatch(130+1)=P_load(90+1)+V(90+1)*(Bus_Ymag(370+1)*cos(theta(90+1)-theta(89+1)-Bus_Yang(370+1))*V(89+1)+Bus_Ymag(372+1)*cos(Bus_Yang(372+1))*V(90+1)+Bus_Ymag(375+1)*cos(theta(90+1)-theta(91+1)-Bus_Yang(375+1))*V(91+1)); % Depth 1, Step 131
  mismatch(131+1)=P_load(91+1)+V(91+1)*(Bus_Ymag(367+1)*cos(-theta(91+1)+theta(88+1)+Bus_Yang(367+1))*V(88+1)+Bus_Ymag(373+1)*cos(-theta(91+1)+theta(90+1)+Bus_Yang(373+1))*V(90+1)+Bus_Ymag(376+1)*cos(Bus_Yang(376+1))*V(91+1)+Bus_Ymag(381+1)*cos(-theta(91+1)+theta(92+1)+Bus_Yang(381+1))*V(92+1)+Bus_Ymag(384+1)*cos(theta(91+1)-theta(93+1)-Bus_Yang(384+1))*V(93+1)+Bus_Ymag(408+1)*cos(theta(91+1)-theta(99+1)-Bus_Yang(408+1))*V(99+1)+Bus_Ymag(420+1)*cos(-theta(91+1)+theta(101+1)+Bus_Yang(420+1))*V(101+1)); % Depth 1, Step 132
  mismatch(132+1)=P_load(92+1)+V(92+1)*(Bus_Ymag(377+1)*cos(theta(92+1)-theta(91+1)-Bus_Yang(377+1))*V(91+1)+Bus_Ymag(382+1)*cos(Bus_Yang(382+1))*V(92+1)+Bus_Ymag(385+1)*cos(theta(92+1)-theta(93+1)-Bus_Yang(385+1))*V(93+1)); % Depth 1, Step 133
  mismatch(133+1)=P_load(93+1)+V(93+1)*(Bus_Ymag(378+1)*cos(-theta(93+1)+theta(91+1)+Bus_Yang(378+1))*V(91+1)+Bus_Ymag(383+1)*cos(-theta(93+1)+theta(92+1)+Bus_Yang(383+1))*V(92+1)+Bus_Ymag(386+1)*cos(Bus_Yang(386+1))*V(93+1)+Bus_Ymag(390+1)*cos(theta(93+1)-theta(94+1)-Bus_Yang(390+1))*V(94+1)+Bus_Ymag(395+1)*cos(theta(93+1)-theta(95+1)-Bus_Yang(395+1))*V(95+1)+Bus_Ymag(409+1)*cos(-theta(93+1)+theta(99+1)+Bus_Yang(409+1))*V(99+1)); % Depth 1, Step 134
  mismatch(134+1)=P_load(94+1)+V(94+1)*(Bus_Ymag(387+1)*cos(-theta(94+1)+theta(93+1)+Bus_Yang(387+1))*V(93+1)+Bus_Ymag(391+1)*cos(Bus_Yang(391+1))*V(94+1)+Bus_Ymag(396+1)*cos(theta(94+1)-theta(95+1)-Bus_Yang(396+1))*V(95+1)); % Depth 1, Step 135
  mismatch(135+1)=P_load(95+1)+V(95+1)*(Bus_Ymag(331+1)*cos(-theta(95+1)+theta(79+1)+Bus_Yang(331+1))*V(79+1)+Bus_Ymag(341+1)*cos(-theta(95+1)+theta(81+1)+Bus_Yang(341+1))*V(81+1)+Bus_Ymag(388+1)*cos(-theta(95+1)+theta(93+1)+Bus_Yang(388+1))*V(93+1)+Bus_Ymag(392+1)*cos(-theta(95+1)+theta(94+1)+Bus_Yang(392+1))*V(94+1)+Bus_Ymag(397+1)*cos(Bus_Yang(397+1))*V(95+1)+Bus_Ymag(400+1)*cos(theta(95+1)-theta(96+1)-Bus_Yang(400+1))*V(96+1)); % Depth 1, Step 136
  mismatch(136+1)=P_load(96+1)+V(96+1)*(Bus_Ymag(332+1)*cos(-theta(96+1)+theta(79+1)+Bus_Yang(332+1))*V(79+1)+Bus_Ymag(398+1)*cos(-theta(96+1)+theta(95+1)+Bus_Yang(398+1))*V(95+1)+Bus_Ymag(401+1)*cos(Bus_Yang(401+1))*V(96+1)); % Depth 1, Step 137
  mismatch(137+1)=P_load(97+1)+V(97+1)*(Bus_Ymag(333+1)*cos(theta(97+1)-theta(79+1)-Bus_Yang(333+1))*V(79+1)+Bus_Ymag(403+1)*cos(Bus_Yang(403+1))*V(97+1)+Bus_Ymag(410+1)*cos(theta(97+1)-theta(99+1)-Bus_Yang(410+1))*V(99+1)); % Depth 1, Step 138
  mismatch(138+1)=P_load(98+1)+V(98+1)*(Bus_Ymag(334+1)*cos(theta(98+1)-theta(79+1)-Bus_Yang(334+1))*V(79+1)+Bus_Ymag(406+1)*cos(Bus_Yang(406+1))*V(98+1)+Bus_Ymag(411+1)*cos(theta(98+1)-theta(99+1)-Bus_Yang(411+1))*V(99+1)); % Depth 1, Step 139
  mismatch(139+1)=P_load(99+1)+V(99+1)*(Bus_Ymag(379+1)*cos(-theta(99+1)+theta(91+1)+Bus_Yang(379+1))*V(91+1)+Bus_Ymag(389+1)*cos(theta(99+1)-theta(93+1)-Bus_Yang(389+1))*V(93+1)+Bus_Ymag(404+1)*cos(-theta(99+1)+theta(97+1)+Bus_Yang(404+1))*V(97+1)+Bus_Ymag(407+1)*cos(-theta(99+1)+theta(98+1)+Bus_Yang(407+1))*V(98+1)+Bus_Ymag(412+1)*cos(Bus_Yang(412+1))*V(99+1)+Bus_Ymag(417+1)*cos(theta(99+1)-theta(100+1)-Bus_Yang(417+1))*V(100+1)+Bus_Ymag(423+1)*cos(-theta(99+1)+theta(102+1)+Bus_Yang(423+1))*V(102+1)+Bus_Ymag(428+1)*cos(-theta(99+1)+theta(103+1)+Bus_Yang(428+1))*V(103+1)+Bus_Ymag(438+1)*cos(-theta(99+1)+theta(105+1)+Bus_Yang(438+1))*V(105+1))-5.5555555555555553581825733999722*V(99+1)*E(18+1)*cos(theta(99+1)-delta_t1(18+1)+1.5707963267948965579989817342721); % Depth 1, Step 140
  mismatch(140+1)=P_load(100+1)+V(100+1)*(Bus_Ymag(413+1)*cos(-theta(100+1)+theta(99+1)+Bus_Yang(413+1))*V(99+1)+Bus_Ymag(418+1)*cos(Bus_Yang(418+1))*V(100+1)+Bus_Ymag(421+1)*cos(-theta(100+1)+theta(101+1)+Bus_Yang(421+1))*V(101+1)); % Depth 1, Step 141
  mismatch(141+1)=P_load(101+1)+V(101+1)*(Bus_Ymag(380+1)*cos(theta(101+1)-theta(91+1)-Bus_Yang(380+1))*V(91+1)+Bus_Ymag(419+1)*cos(theta(101+1)-theta(100+1)-Bus_Yang(419+1))*V(100+1)+Bus_Ymag(422+1)*cos(Bus_Yang(422+1))*V(101+1)); % Depth 1, Step 142
  mismatch(142+1)=P_load(102+1)+V(102+1)*(Bus_Ymag(414+1)*cos(theta(102+1)-theta(99+1)-Bus_Yang(414+1))*V(99+1)+Bus_Ymag(424+1)*cos(Bus_Yang(424+1))*V(102+1)+Bus_Ymag(429+1)*cos(-theta(102+1)+theta(103+1)+Bus_Yang(429+1))*V(103+1)+Bus_Ymag(432+1)*cos(theta(102+1)-theta(104+1)-Bus_Yang(432+1))*V(104+1)+Bus_Ymag(451+1)*cos(theta(102+1)-theta(109+1)-Bus_Yang(451+1))*V(109+1)); % Depth 1, Step 143
  mismatch(143+1)=P_load(103+1)+V(103+1)*(Bus_Ymag(415+1)*cos(theta(103+1)-theta(99+1)-Bus_Yang(415+1))*V(99+1)+Bus_Ymag(425+1)*cos(theta(103+1)-theta(102+1)-Bus_Yang(425+1))*V(102+1)+Bus_Ymag(430+1)*cos(Bus_Yang(430+1))*V(103+1)+Bus_Ymag(433+1)*cos(theta(103+1)-theta(104+1)-Bus_Yang(433+1))*V(104+1)); % Depth 1, Step 144
  mismatch(144+1)=P_load(104+1)+V(104+1)*(Bus_Ymag(426+1)*cos(-theta(104+1)+theta(102+1)+Bus_Yang(426+1))*V(102+1)+Bus_Ymag(431+1)*cos(-theta(104+1)+theta(103+1)+Bus_Yang(431+1))*V(103+1)+Bus_Ymag(434+1)*cos(Bus_Yang(434+1))*V(104+1)+Bus_Ymag(439+1)*cos(-theta(104+1)+theta(105+1)+Bus_Yang(439+1))*V(105+1)+Bus_Ymag(442+1)*cos(theta(104+1)-theta(106+1)-Bus_Yang(442+1))*V(106+1)+Bus_Ymag(445+1)*cos(theta(104+1)-theta(107+1)-Bus_Yang(445+1))*V(107+1)); % Depth 1, Step 145
  mismatch(145+1)=P_load(105+1)+V(105+1)*(Bus_Ymag(416+1)*cos(theta(105+1)-theta(99+1)-Bus_Yang(416+1))*V(99+1)+Bus_Ymag(435+1)*cos(theta(105+1)-theta(104+1)-Bus_Yang(435+1))*V(104+1)+Bus_Ymag(440+1)*cos(Bus_Yang(440+1))*V(105+1)+Bus_Ymag(443+1)*cos(theta(105+1)-theta(106+1)-Bus_Yang(443+1))*V(106+1)); % Depth 1, Step 146
  mismatch(146+1)=P_load(106+1)+V(106+1)*(Bus_Ymag(436+1)*cos(-theta(106+1)+theta(104+1)+Bus_Yang(436+1))*V(104+1)+Bus_Ymag(441+1)*cos(-theta(106+1)+theta(105+1)+Bus_Yang(441+1))*V(105+1)+Bus_Ymag(444+1)*cos(Bus_Yang(444+1))*V(106+1)); % Depth 1, Step 147
  mismatch(147+1)=P_load(107+1)+V(107+1)*(Bus_Ymag(437+1)*cos(-theta(107+1)+theta(104+1)+Bus_Yang(437+1))*V(104+1)+Bus_Ymag(446+1)*cos(Bus_Yang(446+1))*V(107+1)+Bus_Ymag(448+1)*cos(theta(107+1)-theta(108+1)-Bus_Yang(448+1))*V(108+1)); % Depth 1, Step 148
  mismatch(148+1)=P_load(108+1)+V(108+1)*(Bus_Ymag(447+1)*cos(-theta(108+1)+theta(107+1)+Bus_Yang(447+1))*V(107+1)+Bus_Ymag(449+1)*cos(Bus_Yang(449+1))*V(108+1)+Bus_Ymag(452+1)*cos(theta(108+1)-theta(109+1)-Bus_Yang(452+1))*V(109+1)); % Depth 1, Step 149
  mismatch(149+1)=P_load(109+1)+V(109+1)*(Bus_Ymag(427+1)*cos(-theta(109+1)+theta(102+1)+Bus_Yang(427+1))*V(102+1)+Bus_Ymag(450+1)*cos(-theta(109+1)+theta(108+1)+Bus_Yang(450+1))*V(108+1)+Bus_Ymag(453+1)*cos(Bus_Yang(453+1))*V(109+1)+Bus_Ymag(456+1)*cos(theta(109+1)-theta(110+1)-Bus_Yang(456+1))*V(110+1)+Bus_Ymag(458+1)*cos(theta(109+1)-theta(111+1)-Bus_Yang(458+1))*V(111+1)); % Depth 1, Step 150
  mismatch(150+1)=P_load(110+1)+V(110+1)*(Bus_Ymag(454+1)*cos(-theta(110+1)+theta(109+1)+Bus_Yang(454+1))*V(109+1)+Bus_Ymag(457+1)*cos(Bus_Yang(457+1))*V(110+1)); % Depth 1, Step 151
  mismatch(151+1)=P_load(111+1)+V(111+1)*(Bus_Ymag(455+1)*cos(-theta(111+1)+theta(109+1)+Bus_Yang(455+1))*V(109+1)+Bus_Ymag(459+1)*cos(Bus_Yang(459+1))*V(111+1))-5.5555555555555553581825733999722*V(111+1)*E(19+1)*cos(-theta(111+1)+delta_t1(19+1)-1.5707963267948965579989817342721); % Depth 1, Step 152
  mismatch(152+1)=P_load(112+1)+V(112+1)*(Bus_Ymag(68+1)*cos(-theta(112+1)+theta(16+1)+Bus_Yang(68+1))*V(16+1)+Bus_Ymag(126+1)*cos(-theta(112+1)+theta(31+1)+Bus_Yang(126+1))*V(31+1)+Bus_Ymag(462+1)*cos(Bus_Yang(462+1))*V(112+1)); % Depth 1, Step 153
  mismatch(153+1)=P_load(113+1)+V(113+1)*(Bus_Ymag(127+1)*cos(-theta(113+1)+theta(31+1)+Bus_Yang(127+1))*V(31+1)+Bus_Ymag(464+1)*cos(Bus_Yang(464+1))*V(113+1)+Bus_Ymag(467+1)*cos(-theta(113+1)+theta(114+1)+Bus_Yang(467+1))*V(114+1)); % Depth 1, Step 154
  mismatch(154+1)=P_load(114+1)+V(114+1)*(Bus_Ymag(106+1)*cos(-theta(114+1)+theta(26+1)+Bus_Yang(106+1))*V(26+1)+Bus_Ymag(465+1)*cos(-theta(114+1)+theta(113+1)+Bus_Yang(465+1))*V(113+1)+Bus_Ymag(468+1)*cos(Bus_Yang(468+1))*V(114+1)); % Depth 1, Step 155
  mismatch(155+1)=P_load(115+1)+V(115+1)*(Bus_Ymag(279+1)*cos(-theta(115+1)+theta(67+1)+Bus_Yang(279+1))*V(67+1)+Bus_Ymag(470+1)*cos(Bus_Yang(470+1))*V(115+1)); % Depth 1, Step 156
  mismatch(156+1)=P_load(116+1)+V(116+1)*(Bus_Ymag(46+1)*cos(-theta(116+1)+theta(11+1)+Bus_Yang(46+1))*V(11+1)+Bus_Ymag(472+1)*cos(Bus_Yang(472+1))*V(116+1)); % Depth 1, Step 157
  mismatch(157+1)=P_load(117+1)+V(117+1)*(Bus_Ymag(310+1)*cos(-theta(117+1)+theta(74+1)+Bus_Yang(310+1))*V(74+1)+Bus_Ymag(313+1)*cos(-theta(117+1)+theta(75+1)+Bus_Yang(313+1))*V(75+1)+Bus_Ymag(475+1)*cos(Bus_Yang(475+1))*V(117+1)); % Depth 1, Step 158
  mismatch(158+1)=P_load(118+1)+V(118+1)*(Bus_Ymag(173+1)*cos(theta(118+1)-theta(43+1)-Bus_Yang(173+1))*V(43+1)+Bus_Ymag(177+1)*cos(theta(118+1)-theta(44+1)-Bus_Yang(177+1))*V(44+1)+Bus_Ymag(478+1)*cos(Bus_Yang(478+1))*V(118+1)); % Depth 1, Step 159
  mismatch(159+1)=Q_load(0+1)+V(0+1)*(-Bus_Ymag(0+1)*sin(Bus_Yang(0+1))*V(0+1)-Bus_Ymag(3+1)*sin(-theta(0+1)+theta(1+1)+Bus_Yang(3+1))*V(1+1)+Bus_Ymag(6+1)*sin(theta(0+1)-theta(2+1)-Bus_Yang(6+1))*V(2+1)); % Depth 1, Step 160
  mismatch(160+1)=Q_load(1+1)+V(1+1)*(Bus_Ymag(1+1)*sin(theta(1+1)-theta(0+1)-Bus_Yang(1+1))*V(0+1)-Bus_Ymag(4+1)*sin(Bus_Yang(4+1))*V(1+1)+Bus_Ymag(39+1)*sin(theta(1+1)-theta(11+1)-Bus_Yang(39+1))*V(11+1)); % Depth 1, Step 161
  mismatch(161+1)=Q_load(2+1)+V(2+1)*(-Bus_Ymag(2+1)*sin(-theta(2+1)+theta(0+1)+Bus_Yang(2+1))*V(0+1)-Bus_Ymag(7+1)*sin(Bus_Yang(7+1))*V(2+1)-Bus_Ymag(13+1)*sin(-theta(2+1)+theta(4+1)+Bus_Yang(13+1))*V(4+1)+Bus_Ymag(40+1)*sin(theta(2+1)-theta(11+1)-Bus_Yang(40+1))*V(11+1)); % Depth 1, Step 162
  mismatch(162+1)=Q_load(3+1)+V(3+1)*(-Bus_Ymag(10+1)*sin(Bus_Yang(10+1))*V(3+1)-Bus_Ymag(14+1)*sin(-theta(3+1)+theta(4+1)+Bus_Yang(14+1))*V(4+1)+Bus_Ymag(34+1)*sin(theta(3+1)-theta(10+1)-Bus_Yang(34+1))*V(10+1))-5.5555555555555553581825733999722*V(3+1)*E(0+1)*sin(theta(3+1)-delta_t1(0+1)+1.5707963267948965579989817342721); % Depth 1, Step 163
  mismatch(163+1)=Q_load(4+1)+V(4+1)*(Bus_Ymag(8+1)*sin(theta(4+1)-theta(2+1)-Bus_Yang(8+1))*V(2+1)+Bus_Ymag(11+1)*sin(theta(4+1)-theta(3+1)-Bus_Yang(11+1))*V(3+1)-Bus_Ymag(15+1)*sin(Bus_Yang(15+1))*V(4+1)-Bus_Ymag(19+1)*sin(-theta(4+1)+theta(5+1)+Bus_Yang(19+1))*V(5+1)-Bus_Ymag(25+1)*sin(-theta(4+1)+theta(7+1)+Bus_Yang(25+1))*V(7+1)+Bus_Ymag(35+1)*sin(theta(4+1)-theta(10+1)-Bus_Yang(35+1))*V(10+1)); % Depth 1, Step 164
  mismatch(164+1)=Q_load(5+1)+V(5+1)*(Bus_Ymag(16+1)*sin(theta(5+1)-theta(4+1)-Bus_Yang(16+1))*V(4+1)-Bus_Ymag(20+1)*sin(Bus_Yang(20+1))*V(5+1)-Bus_Ymag(22+1)*sin(-theta(5+1)+theta(6+1)+Bus_Yang(22+1))*V(6+1)); % Depth 1, Step 165
  mismatch(165+1)=Q_load(6+1)+V(6+1)*(Bus_Ymag(21+1)*sin(theta(6+1)-theta(5+1)-Bus_Yang(21+1))*V(5+1)-Bus_Ymag(23+1)*sin(Bus_Yang(23+1))*V(6+1)+Bus_Ymag(41+1)*sin(theta(6+1)-theta(11+1)-Bus_Yang(41+1))*V(11+1)); % Depth 1, Step 166
  mismatch(166+1)=Q_load(7+1)+V(7+1)*(Bus_Ymag(17+1)*sin(theta(7+1)-theta(4+1)-Bus_Yang(17+1))*V(4+1)-Bus_Ymag(26+1)*sin(Bus_Yang(26+1))*V(7+1)+Bus_Ymag(29+1)*sin(theta(7+1)-theta(8+1)-Bus_Yang(29+1))*V(8+1)+Bus_Ymag(113+1)*sin(theta(7+1)-theta(29+1)-Bus_Yang(113+1))*V(29+1)); % Depth 1, Step 167
  mismatch(167+1)=Q_load(8+1)+V(8+1)*(-Bus_Ymag(27+1)*sin(-theta(8+1)+theta(7+1)+Bus_Yang(27+1))*V(7+1)-Bus_Ymag(30+1)*sin(Bus_Yang(30+1))*V(8+1)-Bus_Ymag(32+1)*sin(-theta(8+1)+theta(9+1)+Bus_Yang(32+1))*V(9+1)); % Depth 1, Step 168
  mismatch(168+1)=Q_load(9+1)+V(9+1)*(Bus_Ymag(31+1)*sin(theta(9+1)-theta(8+1)-Bus_Yang(31+1))*V(8+1)-Bus_Ymag(33+1)*sin(Bus_Yang(33+1))*V(9+1))-5.5555555555555553581825733999722*V(9+1)*E(1+1)*sin(theta(9+1)-delta_t1(1+1)+1.5707963267948965579989817342721); % Depth 1, Step 169
  mismatch(169+1)=Q_load(10+1)+V(10+1)*(-Bus_Ymag(12+1)*sin(-theta(10+1)+theta(3+1)+Bus_Yang(12+1))*V(3+1)-Bus_Ymag(18+1)*sin(-theta(10+1)+theta(4+1)+Bus_Yang(18+1))*V(4+1)-Bus_Ymag(36+1)*sin(Bus_Yang(36+1))*V(10+1)-Bus_Ymag(42+1)*sin(-theta(10+1)+theta(11+1)+Bus_Yang(42+1))*V(11+1)-Bus_Ymag(47+1)*sin(-theta(10+1)+theta(12+1)+Bus_Yang(47+1))*V(12+1)); % Depth 1, Step 170
  mismatch(170+1)=Q_load(11+1)+V(11+1)*(-Bus_Ymag(5+1)*sin(-theta(11+1)+theta(1+1)+Bus_Yang(5+1))*V(1+1)-Bus_Ymag(9+1)*sin(-theta(11+1)+theta(2+1)+Bus_Yang(9+1))*V(2+1)-Bus_Ymag(24+1)*sin(-theta(11+1)+theta(6+1)+Bus_Yang(24+1))*V(6+1)+Bus_Ymag(37+1)*sin(theta(11+1)-theta(10+1)-Bus_Yang(37+1))*V(10+1)-Bus_Ymag(43+1)*sin(Bus_Yang(43+1))*V(11+1)-Bus_Ymag(50+1)*sin(-theta(11+1)+theta(13+1)+Bus_Yang(50+1))*V(13+1)-Bus_Ymag(59+1)*sin(-theta(11+1)+theta(15+1)+Bus_Yang(59+1))*V(15+1)+Bus_Ymag(471+1)*sin(theta(11+1)-theta(116+1)-Bus_Yang(471+1))*V(116+1))+5.5555555555555553581825733999722*V(11+1)*E(2+1)*sin(-theta(11+1)+delta_t1(2+1)-1.5707963267948965579989817342721); % Depth 1, Step 171
  mismatch(171+1)=Q_load(12+1)+V(12+1)*(Bus_Ymag(38+1)*sin(theta(12+1)-theta(10+1)-Bus_Yang(38+1))*V(10+1)-Bus_Ymag(48+1)*sin(Bus_Yang(48+1))*V(12+1)-Bus_Ymag(53+1)*sin(-theta(12+1)+theta(14+1)+Bus_Yang(53+1))*V(14+1)); % Depth 1, Step 172
  mismatch(172+1)=Q_load(13+1)+V(13+1)*(Bus_Ymag(44+1)*sin(theta(13+1)-theta(11+1)-Bus_Yang(44+1))*V(11+1)-Bus_Ymag(51+1)*sin(Bus_Yang(51+1))*V(13+1)-Bus_Ymag(54+1)*sin(-theta(13+1)+theta(14+1)+Bus_Yang(54+1))*V(14+1)); % Depth 1, Step 173
  mismatch(173+1)=Q_load(14+1)+V(14+1)*(Bus_Ymag(49+1)*sin(theta(14+1)-theta(12+1)-Bus_Yang(49+1))*V(12+1)+Bus_Ymag(52+1)*sin(theta(14+1)-theta(13+1)-Bus_Yang(52+1))*V(13+1)-Bus_Ymag(55+1)*sin(Bus_Yang(55+1))*V(14+1)+Bus_Ymag(62+1)*sin(theta(14+1)-theta(16+1)-Bus_Yang(62+1))*V(16+1)+Bus_Ymag(72+1)*sin(theta(14+1)-theta(18+1)-Bus_Yang(72+1))*V(18+1)+Bus_Ymag(128+1)*sin(theta(14+1)-theta(32+1)-Bus_Yang(128+1))*V(32+1)); % Depth 1, Step 174
  mismatch(174+1)=Q_load(15+1)+V(15+1)*(Bus_Ymag(45+1)*sin(theta(15+1)-theta(11+1)-Bus_Yang(45+1))*V(11+1)-Bus_Ymag(60+1)*sin(Bus_Yang(60+1))*V(15+1)+Bus_Ymag(63+1)*sin(theta(15+1)-theta(16+1)-Bus_Yang(63+1))*V(16+1)); % Depth 1, Step 175
  mismatch(175+1)=Q_load(16+1)+V(16+1)*(-Bus_Ymag(56+1)*sin(-theta(16+1)+theta(14+1)+Bus_Yang(56+1))*V(14+1)-Bus_Ymag(61+1)*sin(-theta(16+1)+theta(15+1)+Bus_Yang(61+1))*V(15+1)-Bus_Ymag(64+1)*sin(Bus_Yang(64+1))*V(16+1)-Bus_Ymag(69+1)*sin(-theta(16+1)+theta(17+1)+Bus_Yang(69+1))*V(17+1)+Bus_Ymag(114+1)*sin(theta(16+1)-theta(29+1)-Bus_Yang(114+1))*V(29+1)+Bus_Ymag(118+1)*sin(theta(16+1)-theta(30+1)-Bus_Yang(118+1))*V(30+1)+Bus_Ymag(460+1)*sin(theta(16+1)-theta(112+1)-Bus_Yang(460+1))*V(112+1)); % Depth 1, Step 176
  mismatch(176+1)=Q_load(17+1)+V(17+1)*(Bus_Ymag(65+1)*sin(theta(17+1)-theta(16+1)-Bus_Yang(65+1))*V(16+1)-Bus_Ymag(70+1)*sin(Bus_Yang(70+1))*V(17+1)-Bus_Ymag(73+1)*sin(-theta(17+1)+theta(18+1)+Bus_Yang(73+1))*V(18+1)); % Depth 1, Step 177
  mismatch(177+1)=Q_load(18+1)+V(18+1)*(-Bus_Ymag(57+1)*sin(-theta(18+1)+theta(14+1)+Bus_Yang(57+1))*V(14+1)+Bus_Ymag(71+1)*sin(theta(18+1)-theta(17+1)-Bus_Yang(71+1))*V(17+1)-Bus_Ymag(74+1)*sin(Bus_Yang(74+1))*V(18+1)-Bus_Ymag(77+1)*sin(-theta(18+1)+theta(19+1)+Bus_Yang(77+1))*V(19+1)+Bus_Ymag(131+1)*sin(theta(18+1)-theta(33+1)-Bus_Yang(131+1))*V(33+1)); % Depth 1, Step 178
  mismatch(178+1)=Q_load(19+1)+V(19+1)*(Bus_Ymag(75+1)*sin(theta(19+1)-theta(18+1)-Bus_Yang(75+1))*V(18+1)-Bus_Ymag(78+1)*sin(Bus_Yang(78+1))*V(19+1)-Bus_Ymag(80+1)*sin(-theta(19+1)+theta(20+1)+Bus_Yang(80+1))*V(20+1)); % Depth 1, Step 179
  mismatch(179+1)=Q_load(20+1)+V(20+1)*(Bus_Ymag(79+1)*sin(theta(20+1)-theta(19+1)-Bus_Yang(79+1))*V(19+1)-Bus_Ymag(81+1)*sin(Bus_Yang(81+1))*V(20+1)-Bus_Ymag(83+1)*sin(-theta(20+1)+theta(21+1)+Bus_Yang(83+1))*V(21+1)); % Depth 1, Step 180
  mismatch(180+1)=Q_load(21+1)+V(21+1)*(Bus_Ymag(82+1)*sin(theta(21+1)-theta(20+1)-Bus_Yang(82+1))*V(20+1)-Bus_Ymag(84+1)*sin(Bus_Yang(84+1))*V(21+1)-Bus_Ymag(86+1)*sin(-theta(21+1)+theta(22+1)+Bus_Yang(86+1))*V(22+1)); % Depth 1, Step 181
  mismatch(181+1)=Q_load(22+1)+V(22+1)*(Bus_Ymag(85+1)*sin(theta(22+1)-theta(21+1)-Bus_Yang(85+1))*V(21+1)-Bus_Ymag(87+1)*sin(Bus_Yang(87+1))*V(22+1)-Bus_Ymag(91+1)*sin(-theta(22+1)+theta(23+1)+Bus_Yang(91+1))*V(23+1)-Bus_Ymag(95+1)*sin(-theta(22+1)+theta(24+1)+Bus_Yang(95+1))*V(24+1)+Bus_Ymag(122+1)*sin(theta(22+1)-theta(31+1)-Bus_Yang(122+1))*V(31+1)); % Depth 1, Step 182
  mismatch(182+1)=Q_load(23+1)+V(23+1)*(Bus_Ymag(88+1)*sin(theta(23+1)-theta(22+1)-Bus_Yang(88+1))*V(22+1)-Bus_Ymag(92+1)*sin(Bus_Yang(92+1))*V(23+1)-Bus_Ymag(287+1)*sin(-theta(23+1)+theta(69+1)+Bus_Yang(287+1))*V(69+1)-Bus_Ymag(297+1)*sin(-theta(23+1)+theta(71+1)+Bus_Yang(297+1))*V(71+1)); % Depth 1, Step 183
  mismatch(183+1)=Q_load(24+1)+V(24+1)*(Bus_Ymag(89+1)*sin(theta(24+1)-theta(22+1)-Bus_Yang(89+1))*V(22+1)-Bus_Ymag(96+1)*sin(Bus_Yang(96+1))*V(24+1)-Bus_Ymag(99+1)*sin(-theta(24+1)+theta(25+1)+Bus_Yang(99+1))*V(25+1)-Bus_Ymag(102+1)*sin(-theta(24+1)+theta(26+1)+Bus_Yang(102+1))*V(26+1))+5.5555555555555553581825733999722*V(24+1)*E(3+1)*sin(-theta(24+1)+delta_t1(3+1)-1.5707963267948965579989817342721); % Depth 1, Step 184
  mismatch(184+1)=Q_load(25+1)+V(25+1)*(Bus_Ymag(97+1)*sin(theta(25+1)-theta(24+1)-Bus_Yang(97+1))*V(24+1)-Bus_Ymag(100+1)*sin(Bus_Yang(100+1))*V(25+1)+Bus_Ymag(115+1)*sin(theta(25+1)-theta(29+1)-Bus_Yang(115+1))*V(29+1))+5.5555555555555553581825733999722*V(25+1)*E(4+1)*sin(-theta(25+1)+delta_t1(4+1)-1.5707963267948965579989817342721); % Depth 1, Step 185
  mismatch(185+1)=Q_load(26+1)+V(26+1)*(Bus_Ymag(98+1)*sin(theta(26+1)-theta(24+1)-Bus_Yang(98+1))*V(24+1)-Bus_Ymag(103+1)*sin(Bus_Yang(103+1))*V(26+1)-Bus_Ymag(107+1)*sin(-theta(26+1)+theta(27+1)+Bus_Yang(107+1))*V(27+1)+Bus_Ymag(123+1)*sin(theta(26+1)-theta(31+1)-Bus_Yang(123+1))*V(31+1)+Bus_Ymag(466+1)*sin(theta(26+1)-theta(114+1)-Bus_Yang(466+1))*V(114+1))+5.5555555555555553581825733999722*V(26+1)*E(5+1)*sin(-theta(26+1)+delta_t1(5+1)-1.5707963267948965579989817342721); % Depth 1, Step 186
  mismatch(186+1)=Q_load(27+1)+V(27+1)*(Bus_Ymag(104+1)*sin(theta(27+1)-theta(26+1)-Bus_Yang(104+1))*V(26+1)-Bus_Ymag(108+1)*sin(Bus_Yang(108+1))*V(27+1)-Bus_Ymag(110+1)*sin(-theta(27+1)+theta(28+1)+Bus_Yang(110+1))*V(28+1)); % Depth 1, Step 187
  mismatch(187+1)=Q_load(28+1)+V(28+1)*(Bus_Ymag(109+1)*sin(theta(28+1)-theta(27+1)-Bus_Yang(109+1))*V(27+1)-Bus_Ymag(111+1)*sin(Bus_Yang(111+1))*V(28+1)+Bus_Ymag(119+1)*sin(theta(28+1)-theta(30+1)-Bus_Yang(119+1))*V(30+1)); % Depth 1, Step 188
  mismatch(188+1)=Q_load(29+1)+V(29+1)*(-Bus_Ymag(28+1)*sin(-theta(29+1)+theta(7+1)+Bus_Yang(28+1))*V(7+1)-Bus_Ymag(66+1)*sin(-theta(29+1)+theta(16+1)+Bus_Yang(66+1))*V(16+1)-Bus_Ymag(101+1)*sin(-theta(29+1)+theta(25+1)+Bus_Yang(101+1))*V(25+1)-Bus_Ymag(116+1)*sin(Bus_Yang(116+1))*V(29+1)-Bus_Ymag(149+1)*sin(-theta(29+1)+theta(37+1)+Bus_Yang(149+1))*V(37+1)); % Depth 1, Step 189
  mismatch(189+1)=Q_load(30+1)+V(30+1)*(-Bus_Ymag(67+1)*sin(-theta(30+1)+theta(16+1)+Bus_Yang(67+1))*V(16+1)-Bus_Ymag(112+1)*sin(-theta(30+1)+theta(28+1)+Bus_Yang(112+1))*V(28+1)-Bus_Ymag(120+1)*sin(Bus_Yang(120+1))*V(30+1)-Bus_Ymag(124+1)*sin(-theta(30+1)+theta(31+1)+Bus_Yang(124+1))*V(31+1))+5.5555555555555553581825733999722*V(30+1)*E(6+1)*sin(-theta(30+1)+delta_t1(6+1)-1.5707963267948965579989817342721); % Depth 1, Step 190
  mismatch(190+1)=Q_load(31+1)+V(31+1)*(-Bus_Ymag(90+1)*sin(-theta(31+1)+theta(22+1)+Bus_Yang(90+1))*V(22+1)-Bus_Ymag(105+1)*sin(-theta(31+1)+theta(26+1)+Bus_Yang(105+1))*V(26+1)+Bus_Ymag(121+1)*sin(theta(31+1)-theta(30+1)-Bus_Yang(121+1))*V(30+1)-Bus_Ymag(125+1)*sin(Bus_Yang(125+1))*V(31+1)+Bus_Ymag(461+1)*sin(theta(31+1)-theta(112+1)-Bus_Yang(461+1))*V(112+1)+Bus_Ymag(463+1)*sin(theta(31+1)-theta(113+1)-Bus_Yang(463+1))*V(113+1)); % Depth 1, Step 191
  mismatch(191+1)=Q_load(32+1)+V(32+1)*(-Bus_Ymag(58+1)*sin(-theta(32+1)+theta(14+1)+Bus_Yang(58+1))*V(14+1)-Bus_Ymag(129+1)*sin(Bus_Yang(129+1))*V(32+1)-Bus_Ymag(142+1)*sin(-theta(32+1)+theta(36+1)+Bus_Yang(142+1))*V(36+1)); % Depth 1, Step 192
  mismatch(192+1)=Q_load(33+1)+V(33+1)*(-Bus_Ymag(76+1)*sin(-theta(33+1)+theta(18+1)+Bus_Yang(76+1))*V(18+1)-Bus_Ymag(132+1)*sin(Bus_Yang(132+1))*V(33+1)-Bus_Ymag(139+1)*sin(-theta(33+1)+theta(35+1)+Bus_Yang(139+1))*V(35+1)-Bus_Ymag(143+1)*sin(-theta(33+1)+theta(36+1)+Bus_Yang(143+1))*V(36+1)+Bus_Ymag(168+1)*sin(theta(33+1)-theta(42+1)-Bus_Yang(168+1))*V(42+1)); % Depth 1, Step 193
  mismatch(193+1)=Q_load(34+1)+V(34+1)*(-Bus_Ymag(136+1)*sin(Bus_Yang(136+1))*V(34+1)-Bus_Ymag(140+1)*sin(-theta(34+1)+theta(35+1)+Bus_Yang(140+1))*V(35+1)-Bus_Ymag(144+1)*sin(-theta(34+1)+theta(36+1)+Bus_Yang(144+1))*V(36+1)); % Depth 1, Step 194
  mismatch(194+1)=Q_load(35+1)+V(35+1)*(Bus_Ymag(133+1)*sin(theta(35+1)-theta(33+1)-Bus_Yang(133+1))*V(33+1)+Bus_Ymag(137+1)*sin(theta(35+1)-theta(34+1)-Bus_Yang(137+1))*V(34+1)-Bus_Ymag(141+1)*sin(Bus_Yang(141+1))*V(35+1)); % Depth 1, Step 195
  mismatch(195+1)=Q_load(36+1)+V(36+1)*(Bus_Ymag(130+1)*sin(theta(36+1)-theta(32+1)-Bus_Yang(130+1))*V(32+1)+Bus_Ymag(134+1)*sin(theta(36+1)-theta(33+1)-Bus_Yang(134+1))*V(33+1)+Bus_Ymag(138+1)*sin(theta(36+1)-theta(34+1)-Bus_Yang(138+1))*V(34+1)-Bus_Ymag(145+1)*sin(Bus_Yang(145+1))*V(36+1)-Bus_Ymag(150+1)*sin(-theta(36+1)+theta(37+1)+Bus_Yang(150+1))*V(37+1)+Bus_Ymag(153+1)*sin(theta(36+1)-theta(38+1)-Bus_Yang(153+1))*V(38+1)+Bus_Ymag(156+1)*sin(theta(36+1)-theta(39+1)-Bus_Yang(156+1))*V(39+1)); % Depth 1, Step 196
  mismatch(196+1)=Q_load(37+1)+V(37+1)*(Bus_Ymag(117+1)*sin(theta(37+1)-theta(29+1)-Bus_Yang(117+1))*V(29+1)+Bus_Ymag(146+1)*sin(theta(37+1)-theta(36+1)-Bus_Yang(146+1))*V(36+1)-Bus_Ymag(151+1)*sin(Bus_Yang(151+1))*V(37+1)+Bus_Ymag(262+1)*sin(theta(37+1)-theta(64+1)-Bus_Yang(262+1))*V(64+1)); % Depth 1, Step 197
  mismatch(197+1)=Q_load(38+1)+V(38+1)*(-Bus_Ymag(147+1)*sin(-theta(38+1)+theta(36+1)+Bus_Yang(147+1))*V(36+1)-Bus_Ymag(154+1)*sin(Bus_Yang(154+1))*V(38+1)-Bus_Ymag(157+1)*sin(-theta(38+1)+theta(39+1)+Bus_Yang(157+1))*V(39+1)); % Depth 1, Step 198
  mismatch(198+1)=Q_load(39+1)+V(39+1)*(-Bus_Ymag(148+1)*sin(-theta(39+1)+theta(36+1)+Bus_Yang(148+1))*V(36+1)+Bus_Ymag(155+1)*sin(theta(39+1)-theta(38+1)-Bus_Yang(155+1))*V(38+1)-Bus_Ymag(158+1)*sin(Bus_Yang(158+1))*V(39+1)-Bus_Ymag(161+1)*sin(-theta(39+1)+theta(40+1)+Bus_Yang(161+1))*V(40+1)-Bus_Ymag(164+1)*sin(-theta(39+1)+theta(41+1)+Bus_Yang(164+1))*V(41+1))-5.5555555555555553581825733999722*V(39+1)*E(7+1)*sin(theta(39+1)-delta_t1(7+1)+1.5707963267948965579989817342721); % Depth 1, Step 199
  mismatch(199+1)=Q_load(40+1)+V(40+1)*(Bus_Ymag(159+1)*sin(theta(40+1)-theta(39+1)-Bus_Yang(159+1))*V(39+1)-Bus_Ymag(162+1)*sin(Bus_Yang(162+1))*V(40+1)-Bus_Ymag(165+1)*sin(-theta(40+1)+theta(41+1)+Bus_Yang(165+1))*V(41+1)); % Depth 1, Step 200
  mismatch(200+1)=Q_load(41+1)+V(41+1)*(Bus_Ymag(160+1)*sin(theta(41+1)-theta(39+1)-Bus_Yang(160+1))*V(39+1)+Bus_Ymag(163+1)*sin(theta(41+1)-theta(40+1)-Bus_Yang(163+1))*V(40+1)-Bus_Ymag(166+1)*sin(Bus_Yang(166+1))*V(41+1)+Bus_Ymag(189+1)*sin(theta(41+1)-theta(48+1)-Bus_Yang(189+1))*V(48+1))-5.5555555555555553581825733999722*V(41+1)*E(8+1)*sin(theta(41+1)-delta_t1(8+1)+1.5707963267948965579989817342721); % Depth 1, Step 201
  mismatch(201+1)=Q_load(42+1)+V(42+1)*(-Bus_Ymag(135+1)*sin(-theta(42+1)+theta(33+1)+Bus_Yang(135+1))*V(33+1)-Bus_Ymag(169+1)*sin(Bus_Yang(169+1))*V(42+1)-Bus_Ymag(171+1)*sin(-theta(42+1)+theta(43+1)+Bus_Yang(171+1))*V(43+1)); % Depth 1, Step 202
  mismatch(202+1)=Q_load(43+1)+V(43+1)*(Bus_Ymag(170+1)*sin(theta(43+1)-theta(42+1)-Bus_Yang(170+1))*V(42+1)-Bus_Ymag(172+1)*sin(Bus_Yang(172+1))*V(43+1)-Bus_Ymag(476+1)*sin(-theta(43+1)+theta(118+1)+Bus_Yang(476+1))*V(118+1)); % Depth 1, Step 203
  mismatch(203+1)=Q_load(44+1)+V(44+1)*(-Bus_Ymag(174+1)*sin(Bus_Yang(174+1))*V(44+1)-Bus_Ymag(178+1)*sin(-theta(44+1)+theta(45+1)+Bus_Yang(178+1))*V(45+1)+Bus_Ymag(190+1)*sin(theta(44+1)-theta(48+1)-Bus_Yang(190+1))*V(48+1)-Bus_Ymag(477+1)*sin(-theta(44+1)+theta(118+1)+Bus_Yang(477+1))*V(118+1)); % Depth 1, Step 204
  mismatch(204+1)=Q_load(45+1)+V(45+1)*(Bus_Ymag(175+1)*sin(theta(45+1)-theta(44+1)-Bus_Yang(175+1))*V(44+1)-Bus_Ymag(179+1)*sin(Bus_Yang(179+1))*V(45+1)-Bus_Ymag(182+1)*sin(-theta(45+1)+theta(46+1)+Bus_Yang(182+1))*V(46+1)-Bus_Ymag(186+1)*sin(-theta(45+1)+theta(47+1)+Bus_Yang(186+1))*V(47+1)); % Depth 1, Step 205
  mismatch(205+1)=Q_load(46+1)+V(46+1)*(Bus_Ymag(180+1)*sin(theta(46+1)-theta(45+1)-Bus_Yang(180+1))*V(45+1)-Bus_Ymag(183+1)*sin(Bus_Yang(183+1))*V(46+1)+Bus_Ymag(191+1)*sin(theta(46+1)-theta(48+1)-Bus_Yang(191+1))*V(48+1)-Bus_Ymag(280+1)*sin(-theta(46+1)+theta(68+1)+Bus_Yang(280+1))*V(68+1)); % Depth 1, Step 206
  mismatch(206+1)=Q_load(47+1)+V(47+1)*(Bus_Ymag(181+1)*sin(theta(47+1)-theta(45+1)-Bus_Yang(181+1))*V(45+1)-Bus_Ymag(187+1)*sin(Bus_Yang(187+1))*V(47+1)+Bus_Ymag(192+1)*sin(theta(47+1)-theta(48+1)-Bus_Yang(192+1))*V(48+1)); % Depth 1, Step 207
  mismatch(207+1)=Q_load(48+1)+V(48+1)*(-Bus_Ymag(167+1)*sin(-theta(48+1)+theta(41+1)+Bus_Yang(167+1))*V(41+1)-Bus_Ymag(176+1)*sin(-theta(48+1)+theta(44+1)+Bus_Yang(176+1))*V(44+1)-Bus_Ymag(184+1)*sin(-theta(48+1)+theta(46+1)+Bus_Yang(184+1))*V(46+1)-Bus_Ymag(188+1)*sin(-theta(48+1)+theta(47+1)+Bus_Yang(188+1))*V(47+1)-Bus_Ymag(193+1)*sin(Bus_Yang(193+1))*V(48+1)-Bus_Ymag(199+1)*sin(-theta(48+1)+theta(49+1)+Bus_Yang(199+1))*V(49+1)-Bus_Ymag(202+1)*sin(-theta(48+1)+theta(50+1)+Bus_Yang(202+1))*V(50+1)+Bus_Ymag(212+1)*sin(theta(48+1)-theta(53+1)-Bus_Yang(212+1))*V(53+1)-Bus_Ymag(267+1)*sin(-theta(48+1)+theta(65+1)+Bus_Yang(267+1))*V(65+1)-Bus_Ymag(281+1)*sin(-theta(48+1)+theta(68+1)+Bus_Yang(281+1))*V(68+1))-5.5555555555555553581825733999722*V(48+1)*E(9+1)*sin(theta(48+1)-delta_t1(9+1)+1.5707963267948965579989817342721); % Depth 1, Step 208
  mismatch(208+1)=Q_load(49+1)+V(49+1)*(Bus_Ymag(194+1)*sin(theta(49+1)-theta(48+1)-Bus_Yang(194+1))*V(48+1)-Bus_Ymag(200+1)*sin(Bus_Yang(200+1))*V(49+1)+Bus_Ymag(228+1)*sin(theta(49+1)-theta(56+1)-Bus_Yang(228+1))*V(56+1)); % Depth 1, Step 209
  mismatch(209+1)=Q_load(50+1)+V(50+1)*(Bus_Ymag(195+1)*sin(theta(50+1)-theta(48+1)-Bus_Yang(195+1))*V(48+1)-Bus_Ymag(203+1)*sin(Bus_Yang(203+1))*V(50+1)-Bus_Ymag(206+1)*sin(-theta(50+1)+theta(51+1)+Bus_Yang(206+1))*V(51+1)+Bus_Ymag(231+1)*sin(theta(50+1)-theta(57+1)-Bus_Yang(231+1))*V(57+1)); % Depth 1, Step 210
  mismatch(210+1)=Q_load(51+1)+V(51+1)*(Bus_Ymag(204+1)*sin(theta(51+1)-theta(50+1)-Bus_Yang(204+1))*V(50+1)-Bus_Ymag(207+1)*sin(Bus_Yang(207+1))*V(51+1)-Bus_Ymag(209+1)*sin(-theta(51+1)+theta(52+1)+Bus_Yang(209+1))*V(52+1)); % Depth 1, Step 211
  mismatch(211+1)=Q_load(52+1)+V(52+1)*(Bus_Ymag(208+1)*sin(theta(52+1)-theta(51+1)-Bus_Yang(208+1))*V(51+1)-Bus_Ymag(210+1)*sin(Bus_Yang(210+1))*V(52+1)+Bus_Ymag(213+1)*sin(theta(52+1)-theta(53+1)-Bus_Yang(213+1))*V(53+1)); % Depth 1, Step 212
  mismatch(212+1)=Q_load(53+1)+V(53+1)*(-Bus_Ymag(196+1)*sin(-theta(53+1)+theta(48+1)+Bus_Yang(196+1))*V(48+1)-Bus_Ymag(211+1)*sin(-theta(53+1)+theta(52+1)+Bus_Yang(211+1))*V(52+1)-Bus_Ymag(214+1)*sin(Bus_Yang(214+1))*V(53+1)-Bus_Ymag(218+1)*sin(-theta(53+1)+theta(54+1)+Bus_Yang(218+1))*V(54+1)-Bus_Ymag(222+1)*sin(-theta(53+1)+theta(55+1)+Bus_Yang(222+1))*V(55+1)-Bus_Ymag(234+1)*sin(-theta(53+1)+theta(58+1)+Bus_Yang(234+1))*V(58+1)); % Depth 1, Step 213
  mismatch(213+1)=Q_load(54+1)+V(54+1)*(Bus_Ymag(215+1)*sin(theta(54+1)-theta(53+1)-Bus_Yang(215+1))*V(53+1)-Bus_Ymag(219+1)*sin(Bus_Yang(219+1))*V(54+1)-Bus_Ymag(223+1)*sin(-theta(54+1)+theta(55+1)+Bus_Yang(223+1))*V(55+1)-Bus_Ymag(235+1)*sin(-theta(54+1)+theta(58+1)+Bus_Yang(235+1))*V(58+1)); % Depth 1, Step 214
  mismatch(214+1)=Q_load(55+1)+V(55+1)*(Bus_Ymag(216+1)*sin(theta(55+1)-theta(53+1)-Bus_Yang(216+1))*V(53+1)+Bus_Ymag(220+1)*sin(theta(55+1)-theta(54+1)-Bus_Yang(220+1))*V(54+1)-Bus_Ymag(224+1)*sin(Bus_Yang(224+1))*V(55+1)-Bus_Ymag(229+1)*sin(-theta(55+1)+theta(56+1)+Bus_Yang(229+1))*V(56+1)-Bus_Ymag(232+1)*sin(-theta(55+1)+theta(57+1)+Bus_Yang(232+1))*V(57+1)-Bus_Ymag(236+1)*sin(-theta(55+1)+theta(58+1)+Bus_Yang(236+1))*V(58+1)); % Depth 1, Step 215
  mismatch(215+1)=Q_load(56+1)+V(56+1)*(-Bus_Ymag(201+1)*sin(-theta(56+1)+theta(49+1)+Bus_Yang(201+1))*V(49+1)+Bus_Ymag(225+1)*sin(theta(56+1)-theta(55+1)-Bus_Yang(225+1))*V(55+1)-Bus_Ymag(230+1)*sin(Bus_Yang(230+1))*V(56+1)); % Depth 1, Step 216
  mismatch(216+1)=Q_load(57+1)+V(57+1)*(-Bus_Ymag(205+1)*sin(-theta(57+1)+theta(50+1)+Bus_Yang(205+1))*V(50+1)+Bus_Ymag(226+1)*sin(theta(57+1)-theta(55+1)-Bus_Yang(226+1))*V(55+1)-Bus_Ymag(233+1)*sin(Bus_Yang(233+1))*V(57+1)); % Depth 1, Step 217
  mismatch(217+1)=Q_load(58+1)+V(58+1)*(Bus_Ymag(217+1)*sin(theta(58+1)-theta(53+1)-Bus_Yang(217+1))*V(53+1)+Bus_Ymag(221+1)*sin(theta(58+1)-theta(54+1)-Bus_Yang(221+1))*V(54+1)+Bus_Ymag(227+1)*sin(theta(58+1)-theta(55+1)-Bus_Yang(227+1))*V(55+1)-Bus_Ymag(237+1)*sin(Bus_Yang(237+1))*V(58+1)+Bus_Ymag(241+1)*sin(theta(58+1)-theta(59+1)-Bus_Yang(241+1))*V(59+1)+Bus_Ymag(245+1)*sin(theta(58+1)-theta(60+1)-Bus_Yang(245+1))*V(60+1)-Bus_Ymag(255+1)*sin(-theta(58+1)+theta(62+1)+Bus_Yang(255+1))*V(62+1))-5.5555555555555553581825733999722*V(58+1)*E(10+1)*sin(theta(58+1)-delta_t1(10+1)+1.5707963267948965579989817342721); % Depth 1, Step 218
  mismatch(218+1)=Q_load(59+1)+V(59+1)*(-Bus_Ymag(238+1)*sin(-theta(59+1)+theta(58+1)+Bus_Yang(238+1))*V(58+1)-Bus_Ymag(242+1)*sin(Bus_Yang(242+1))*V(59+1)+Bus_Ymag(246+1)*sin(theta(59+1)-theta(60+1)-Bus_Yang(246+1))*V(60+1)+Bus_Ymag(250+1)*sin(theta(59+1)-theta(61+1)-Bus_Yang(250+1))*V(61+1)); % Depth 1, Step 219
  mismatch(219+1)=Q_load(60+1)+V(60+1)*(-Bus_Ymag(239+1)*sin(-theta(60+1)+theta(58+1)+Bus_Yang(239+1))*V(58+1)-Bus_Ymag(243+1)*sin(-theta(60+1)+theta(59+1)+Bus_Yang(243+1))*V(59+1)-Bus_Ymag(247+1)*sin(Bus_Yang(247+1))*V(60+1)+Bus_Ymag(251+1)*sin(theta(60+1)-theta(61+1)-Bus_Yang(251+1))*V(61+1)+Bus_Ymag(258+1)*sin(theta(60+1)-theta(63+1)-Bus_Yang(258+1))*V(63+1))-5.5555555555555553581825733999722*V(60+1)*E(11+1)*sin(theta(60+1)-delta_t1(11+1)+1.5707963267948965579989817342721); % Depth 1, Step 220
  mismatch(220+1)=Q_load(61+1)+V(61+1)*(-Bus_Ymag(244+1)*sin(-theta(61+1)+theta(59+1)+Bus_Yang(244+1))*V(59+1)-Bus_Ymag(248+1)*sin(-theta(61+1)+theta(60+1)+Bus_Yang(248+1))*V(60+1)-Bus_Ymag(252+1)*sin(Bus_Yang(252+1))*V(61+1)-Bus_Ymag(268+1)*sin(-theta(61+1)+theta(65+1)+Bus_Yang(268+1))*V(65+1)-Bus_Ymag(272+1)*sin(-theta(61+1)+theta(66+1)+Bus_Yang(272+1))*V(66+1)); % Depth 1, Step 221
  mismatch(221+1)=Q_load(62+1)+V(62+1)*(Bus_Ymag(240+1)*sin(theta(62+1)-theta(58+1)-Bus_Yang(240+1))*V(58+1)-Bus_Ymag(256+1)*sin(Bus_Yang(256+1))*V(62+1)+Bus_Ymag(259+1)*sin(theta(62+1)-theta(63+1)-Bus_Yang(259+1))*V(63+1)); % Depth 1, Step 222
  mismatch(222+1)=Q_load(63+1)+V(63+1)*(-Bus_Ymag(249+1)*sin(-theta(63+1)+theta(60+1)+Bus_Yang(249+1))*V(60+1)-Bus_Ymag(257+1)*sin(-theta(63+1)+theta(62+1)+Bus_Yang(257+1))*V(62+1)-Bus_Ymag(260+1)*sin(Bus_Yang(260+1))*V(63+1)-Bus_Ymag(263+1)*sin(-theta(63+1)+theta(64+1)+Bus_Yang(263+1))*V(64+1)); % Depth 1, Step 223
  mismatch(223+1)=Q_load(64+1)+V(64+1)*(-Bus_Ymag(152+1)*sin(-theta(64+1)+theta(37+1)+Bus_Yang(152+1))*V(37+1)-Bus_Ymag(261+1)*sin(-theta(64+1)+theta(63+1)+Bus_Yang(261+1))*V(63+1)-Bus_Ymag(264+1)*sin(Bus_Yang(264+1))*V(64+1)-Bus_Ymag(269+1)*sin(-theta(64+1)+theta(65+1)+Bus_Yang(269+1))*V(65+1)-Bus_Ymag(275+1)*sin(-theta(64+1)+theta(67+1)+Bus_Yang(275+1))*V(67+1))+5.5555555555555553581825733999722*V(64+1)*E(12+1)*sin(-theta(64+1)+delta_t1(12+1)-1.5707963267948965579989817342721); % Depth 1, Step 224
  mismatch(224+1)=Q_load(65+1)+V(65+1)*(Bus_Ymag(197+1)*sin(theta(65+1)-theta(48+1)-Bus_Yang(197+1))*V(48+1)+Bus_Ymag(253+1)*sin(theta(65+1)-theta(61+1)-Bus_Yang(253+1))*V(61+1)+Bus_Ymag(265+1)*sin(theta(65+1)-theta(64+1)-Bus_Yang(265+1))*V(64+1)-Bus_Ymag(270+1)*sin(Bus_Yang(270+1))*V(65+1)-Bus_Ymag(273+1)*sin(-theta(65+1)+theta(66+1)+Bus_Yang(273+1))*V(66+1))-5.5555555555555553581825733999722*V(65+1)*E(13+1)*sin(theta(65+1)-delta_t1(13+1)+1.5707963267948965579989817342721); % Depth 1, Step 225
  mismatch(225+1)=Q_load(66+1)+V(66+1)*(Bus_Ymag(254+1)*sin(theta(66+1)-theta(61+1)-Bus_Yang(254+1))*V(61+1)+Bus_Ymag(271+1)*sin(theta(66+1)-theta(65+1)-Bus_Yang(271+1))*V(65+1)-Bus_Ymag(274+1)*sin(Bus_Yang(274+1))*V(66+1)); % Depth 1, Step 226
  mismatch(226+1)=Q_load(67+1)+V(67+1)*(Bus_Ymag(266+1)*sin(theta(67+1)-theta(64+1)-Bus_Yang(266+1))*V(64+1)-Bus_Ymag(276+1)*sin(Bus_Yang(276+1))*V(67+1)-Bus_Ymag(282+1)*sin(-theta(67+1)+theta(68+1)+Bus_Yang(282+1))*V(68+1)+Bus_Ymag(335+1)*sin(theta(67+1)-theta(80+1)-Bus_Yang(335+1))*V(80+1)+Bus_Ymag(469+1)*sin(theta(67+1)-theta(115+1)-Bus_Yang(469+1))*V(115+1)); % Depth 1, Step 227
  mismatch(227+1)=Q_load(68+1)+V(68+1)*(Bus_Ymag(185+1)*sin(theta(68+1)-theta(46+1)-Bus_Yang(185+1))*V(46+1)+Bus_Ymag(198+1)*sin(theta(68+1)-theta(48+1)-Bus_Yang(198+1))*V(48+1)+Bus_Ymag(277+1)*sin(theta(68+1)-theta(67+1)-Bus_Yang(277+1))*V(67+1)-Bus_Ymag(283+1)*sin(Bus_Yang(283+1))*V(68+1)-Bus_Ymag(288+1)*sin(-theta(68+1)+theta(69+1)+Bus_Yang(288+1))*V(69+1)+Bus_Ymag(305+1)*sin(theta(68+1)-theta(74+1)-Bus_Yang(305+1))*V(74+1)+Bus_Ymag(314+1)*sin(theta(68+1)-theta(76+1)-Bus_Yang(314+1))*V(76+1))-5.5555555555555553581825733999722*V(68+1)*E(14+1)*sin(theta(68+1)-delta_t1(14+1)+1.5707963267948965579989817342721); % Depth 1, Step 228
  mismatch(228+1)=Q_load(69+1)+V(69+1)*(Bus_Ymag(93+1)*sin(theta(69+1)-theta(23+1)-Bus_Yang(93+1))*V(23+1)+Bus_Ymag(284+1)*sin(theta(69+1)-theta(68+1)-Bus_Yang(284+1))*V(68+1)-Bus_Ymag(289+1)*sin(Bus_Yang(289+1))*V(69+1)-Bus_Ymag(293+1)*sin(-theta(69+1)+theta(70+1)+Bus_Yang(293+1))*V(70+1)+Bus_Ymag(302+1)*sin(theta(69+1)-theta(73+1)-Bus_Yang(302+1))*V(73+1)+Bus_Ymag(306+1)*sin(theta(69+1)-theta(74+1)-Bus_Yang(306+1))*V(74+1)); % Depth 1, Step 229
  mismatch(229+1)=Q_load(70+1)+V(70+1)*(Bus_Ymag(290+1)*sin(theta(70+1)-theta(69+1)-Bus_Yang(290+1))*V(69+1)-Bus_Ymag(294+1)*sin(Bus_Yang(294+1))*V(70+1)-Bus_Ymag(298+1)*sin(-theta(70+1)+theta(71+1)+Bus_Yang(298+1))*V(71+1)+Bus_Ymag(300+1)*sin(theta(70+1)-theta(72+1)-Bus_Yang(300+1))*V(72+1)); % Depth 1, Step 230
  mismatch(230+1)=Q_load(71+1)+V(71+1)*(Bus_Ymag(94+1)*sin(theta(71+1)-theta(23+1)-Bus_Yang(94+1))*V(23+1)+Bus_Ymag(295+1)*sin(theta(71+1)-theta(70+1)-Bus_Yang(295+1))*V(70+1)-Bus_Ymag(299+1)*sin(Bus_Yang(299+1))*V(71+1)); % Depth 1, Step 231
  mismatch(231+1)=Q_load(72+1)+V(72+1)*(-Bus_Ymag(296+1)*sin(-theta(72+1)+theta(70+1)+Bus_Yang(296+1))*V(70+1)-Bus_Ymag(301+1)*sin(Bus_Yang(301+1))*V(72+1)); % Depth 1, Step 232
  mismatch(232+1)=Q_load(73+1)+V(73+1)*(-Bus_Ymag(291+1)*sin(-theta(73+1)+theta(69+1)+Bus_Yang(291+1))*V(69+1)-Bus_Ymag(303+1)*sin(Bus_Yang(303+1))*V(73+1)-Bus_Ymag(307+1)*sin(-theta(73+1)+theta(74+1)+Bus_Yang(307+1))*V(74+1)); % Depth 1, Step 233
  mismatch(233+1)=Q_load(74+1)+V(74+1)*(-Bus_Ymag(285+1)*sin(-theta(74+1)+theta(68+1)+Bus_Yang(285+1))*V(68+1)-Bus_Ymag(292+1)*sin(-theta(74+1)+theta(69+1)+Bus_Yang(292+1))*V(69+1)+Bus_Ymag(304+1)*sin(theta(74+1)-theta(73+1)-Bus_Yang(304+1))*V(73+1)-Bus_Ymag(308+1)*sin(Bus_Yang(308+1))*V(74+1)-Bus_Ymag(315+1)*sin(-theta(74+1)+theta(76+1)+Bus_Yang(315+1))*V(76+1)+Bus_Ymag(473+1)*sin(theta(74+1)-theta(117+1)-Bus_Yang(473+1))*V(117+1)); % Depth 1, Step 234
  mismatch(234+1)=Q_load(75+1)+V(75+1)*(-Bus_Ymag(311+1)*sin(Bus_Yang(311+1))*V(75+1)-Bus_Ymag(316+1)*sin(-theta(75+1)+theta(76+1)+Bus_Yang(316+1))*V(76+1)+Bus_Ymag(474+1)*sin(theta(75+1)-theta(117+1)-Bus_Yang(474+1))*V(117+1)); % Depth 1, Step 235
  mismatch(235+1)=Q_load(76+1)+V(76+1)*(-Bus_Ymag(286+1)*sin(-theta(76+1)+theta(68+1)+Bus_Yang(286+1))*V(68+1)+Bus_Ymag(309+1)*sin(theta(76+1)-theta(74+1)-Bus_Yang(309+1))*V(74+1)+Bus_Ymag(312+1)*sin(theta(76+1)-theta(75+1)-Bus_Yang(312+1))*V(75+1)-Bus_Ymag(317+1)*sin(Bus_Yang(317+1))*V(76+1)+Bus_Ymag(321+1)*sin(theta(76+1)-theta(77+1)-Bus_Yang(321+1))*V(77+1)-Bus_Ymag(327+1)*sin(-theta(76+1)+theta(79+1)+Bus_Yang(327+1))*V(79+1)+Bus_Ymag(338+1)*sin(theta(76+1)-theta(81+1)-Bus_Yang(338+1))*V(81+1)); % Depth 1, Step 236
  mismatch(236+1)=Q_load(77+1)+V(77+1)*(-Bus_Ymag(318+1)*sin(-theta(77+1)+theta(76+1)+Bus_Yang(318+1))*V(76+1)-Bus_Ymag(322+1)*sin(Bus_Yang(322+1))*V(77+1)-Bus_Ymag(324+1)*sin(-theta(77+1)+theta(78+1)+Bus_Yang(324+1))*V(78+1)); % Depth 1, Step 237
  mismatch(237+1)=Q_load(78+1)+V(78+1)*(Bus_Ymag(323+1)*sin(theta(78+1)-theta(77+1)-Bus_Yang(323+1))*V(77+1)-Bus_Ymag(325+1)*sin(Bus_Yang(325+1))*V(78+1)-Bus_Ymag(328+1)*sin(-theta(78+1)+theta(79+1)+Bus_Yang(328+1))*V(79+1)); % Depth 1, Step 238
  mismatch(238+1)=Q_load(79+1)+V(79+1)*(-Bus_Ymag(319+1)*sin(-theta(79+1)+theta(76+1)+Bus_Yang(319+1))*V(76+1)-Bus_Ymag(326+1)*sin(-theta(79+1)+theta(78+1)+Bus_Yang(326+1))*V(78+1)-Bus_Ymag(329+1)*sin(Bus_Yang(329+1))*V(79+1)-Bus_Ymag(336+1)*sin(-theta(79+1)+theta(80+1)+Bus_Yang(336+1))*V(80+1)+Bus_Ymag(393+1)*sin(theta(79+1)-theta(95+1)-Bus_Yang(393+1))*V(95+1)+Bus_Ymag(399+1)*sin(theta(79+1)-theta(96+1)-Bus_Yang(399+1))*V(96+1)-Bus_Ymag(402+1)*sin(-theta(79+1)+theta(97+1)+Bus_Yang(402+1))*V(97+1)-Bus_Ymag(405+1)*sin(-theta(79+1)+theta(98+1)+Bus_Yang(405+1))*V(98+1))-5.5555555555555553581825733999722*V(79+1)*E(15+1)*sin(theta(79+1)-delta_t1(15+1)+1.5707963267948965579989817342721); % Depth 1, Step 239
  mismatch(239+1)=Q_load(80+1)+V(80+1)*(-Bus_Ymag(278+1)*sin(-theta(80+1)+theta(67+1)+Bus_Yang(278+1))*V(67+1)+Bus_Ymag(330+1)*sin(theta(80+1)-theta(79+1)-Bus_Yang(330+1))*V(79+1)-Bus_Ymag(337+1)*sin(Bus_Yang(337+1))*V(80+1)); % Depth 1, Step 240
  mismatch(240+1)=Q_load(81+1)+V(81+1)*(-Bus_Ymag(320+1)*sin(-theta(81+1)+theta(76+1)+Bus_Yang(320+1))*V(76+1)-Bus_Ymag(339+1)*sin(Bus_Yang(339+1))*V(81+1)-Bus_Ymag(342+1)*sin(-theta(81+1)+theta(82+1)+Bus_Yang(342+1))*V(82+1)+Bus_Ymag(394+1)*sin(theta(81+1)-theta(95+1)-Bus_Yang(394+1))*V(95+1)); % Depth 1, Step 241
  mismatch(241+1)=Q_load(82+1)+V(82+1)*(Bus_Ymag(340+1)*sin(theta(82+1)-theta(81+1)-Bus_Yang(340+1))*V(81+1)-Bus_Ymag(343+1)*sin(Bus_Yang(343+1))*V(82+1)+Bus_Ymag(346+1)*sin(theta(82+1)-theta(83+1)-Bus_Yang(346+1))*V(83+1)+Bus_Ymag(349+1)*sin(theta(82+1)-theta(84+1)-Bus_Yang(349+1))*V(84+1)); % Depth 1, Step 242
  mismatch(242+1)=Q_load(83+1)+V(83+1)*(-Bus_Ymag(344+1)*sin(-theta(83+1)+theta(82+1)+Bus_Yang(344+1))*V(82+1)-Bus_Ymag(347+1)*sin(Bus_Yang(347+1))*V(83+1)-Bus_Ymag(350+1)*sin(-theta(83+1)+theta(84+1)+Bus_Yang(350+1))*V(84+1)); % Depth 1, Step 243
  mismatch(243+1)=Q_load(84+1)+V(84+1)*(-Bus_Ymag(345+1)*sin(-theta(84+1)+theta(82+1)+Bus_Yang(345+1))*V(82+1)+Bus_Ymag(348+1)*sin(theta(84+1)-theta(83+1)-Bus_Yang(348+1))*V(83+1)-Bus_Ymag(351+1)*sin(Bus_Yang(351+1))*V(84+1)-Bus_Ymag(355+1)*sin(-theta(84+1)+theta(85+1)+Bus_Yang(355+1))*V(85+1)-Bus_Ymag(360+1)*sin(-theta(84+1)+theta(87+1)+Bus_Yang(360+1))*V(87+1)-Bus_Ymag(363+1)*sin(-theta(84+1)+theta(88+1)+Bus_Yang(363+1))*V(88+1)); % Depth 1, Step 244
  mismatch(244+1)=Q_load(85+1)+V(85+1)*(Bus_Ymag(352+1)*sin(theta(85+1)-theta(84+1)-Bus_Yang(352+1))*V(84+1)-Bus_Ymag(356+1)*sin(Bus_Yang(356+1))*V(85+1)-Bus_Ymag(358+1)*sin(-theta(85+1)+theta(86+1)+Bus_Yang(358+1))*V(86+1)); % Depth 1, Step 245
  mismatch(245+1)=Q_load(86+1)+V(86+1)*(Bus_Ymag(357+1)*sin(theta(86+1)-theta(85+1)-Bus_Yang(357+1))*V(85+1)-Bus_Ymag(359+1)*sin(Bus_Yang(359+1))*V(86+1)); % Depth 1, Step 246
  mismatch(246+1)=Q_load(87+1)+V(87+1)*(Bus_Ymag(353+1)*sin(theta(87+1)-theta(84+1)-Bus_Yang(353+1))*V(84+1)-Bus_Ymag(361+1)*sin(Bus_Yang(361+1))*V(87+1)-Bus_Ymag(364+1)*sin(-theta(87+1)+theta(88+1)+Bus_Yang(364+1))*V(88+1)); % Depth 1, Step 247
  mismatch(247+1)=Q_load(88+1)+V(88+1)*(Bus_Ymag(354+1)*sin(theta(88+1)-theta(84+1)-Bus_Yang(354+1))*V(84+1)+Bus_Ymag(362+1)*sin(theta(88+1)-theta(87+1)-Bus_Yang(362+1))*V(87+1)-Bus_Ymag(365+1)*sin(Bus_Yang(365+1))*V(88+1)+Bus_Ymag(368+1)*sin(theta(88+1)-theta(89+1)-Bus_Yang(368+1))*V(89+1)+Bus_Ymag(374+1)*sin(theta(88+1)-theta(91+1)-Bus_Yang(374+1))*V(91+1))-5.5555555555555553581825733999722*V(88+1)*E(16+1)*sin(theta(88+1)-delta_t1(16+1)+1.5707963267948965579989817342721); % Depth 1, Step 248
  mismatch(248+1)=Q_load(89+1)+V(89+1)*(-Bus_Ymag(366+1)*sin(-theta(89+1)+theta(88+1)+Bus_Yang(366+1))*V(88+1)-Bus_Ymag(369+1)*sin(Bus_Yang(369+1))*V(89+1)-Bus_Ymag(371+1)*sin(-theta(89+1)+theta(90+1)+Bus_Yang(371+1))*V(90+1))-5.5555555555555553581825733999722*V(89+1)*E(17+1)*sin(theta(89+1)-delta_t1(17+1)+1.5707963267948965579989817342721); % Depth 1, Step 249
  mismatch(249+1)=Q_load(90+1)+V(90+1)*(Bus_Ymag(370+1)*sin(theta(90+1)-theta(89+1)-Bus_Yang(370+1))*V(89+1)-Bus_Ymag(372+1)*sin(Bus_Yang(372+1))*V(90+1)+Bus_Ymag(375+1)*sin(theta(90+1)-theta(91+1)-Bus_Yang(375+1))*V(91+1)); % Depth 1, Step 250
  mismatch(250+1)=Q_load(91+1)+V(91+1)*(-Bus_Ymag(367+1)*sin(-theta(91+1)+theta(88+1)+Bus_Yang(367+1))*V(88+1)-Bus_Ymag(373+1)*sin(-theta(91+1)+theta(90+1)+Bus_Yang(373+1))*V(90+1)-Bus_Ymag(376+1)*sin(Bus_Yang(376+1))*V(91+1)-Bus_Ymag(381+1)*sin(-theta(91+1)+theta(92+1)+Bus_Yang(381+1))*V(92+1)+Bus_Ymag(384+1)*sin(theta(91+1)-theta(93+1)-Bus_Yang(384+1))*V(93+1)+Bus_Ymag(408+1)*sin(theta(91+1)-theta(99+1)-Bus_Yang(408+1))*V(99+1)-Bus_Ymag(420+1)*sin(-theta(91+1)+theta(101+1)+Bus_Yang(420+1))*V(101+1)); % Depth 1, Step 251
  mismatch(251+1)=Q_load(92+1)+V(92+1)*(Bus_Ymag(377+1)*sin(theta(92+1)-theta(91+1)-Bus_Yang(377+1))*V(91+1)-Bus_Ymag(382+1)*sin(Bus_Yang(382+1))*V(92+1)+Bus_Ymag(385+1)*sin(theta(92+1)-theta(93+1)-Bus_Yang(385+1))*V(93+1)); % Depth 1, Step 252
  mismatch(252+1)=Q_load(93+1)+V(93+1)*(-Bus_Ymag(378+1)*sin(-theta(93+1)+theta(91+1)+Bus_Yang(378+1))*V(91+1)-Bus_Ymag(383+1)*sin(-theta(93+1)+theta(92+1)+Bus_Yang(383+1))*V(92+1)-Bus_Ymag(386+1)*sin(Bus_Yang(386+1))*V(93+1)+Bus_Ymag(390+1)*sin(theta(93+1)-theta(94+1)-Bus_Yang(390+1))*V(94+1)+Bus_Ymag(395+1)*sin(theta(93+1)-theta(95+1)-Bus_Yang(395+1))*V(95+1)-Bus_Ymag(409+1)*sin(-theta(93+1)+theta(99+1)+Bus_Yang(409+1))*V(99+1)); % Depth 1, Step 253
  mismatch(253+1)=Q_load(94+1)+V(94+1)*(-Bus_Ymag(387+1)*sin(-theta(94+1)+theta(93+1)+Bus_Yang(387+1))*V(93+1)-Bus_Ymag(391+1)*sin(Bus_Yang(391+1))*V(94+1)+Bus_Ymag(396+1)*sin(theta(94+1)-theta(95+1)-Bus_Yang(396+1))*V(95+1)); % Depth 1, Step 254
  mismatch(254+1)=Q_load(95+1)+V(95+1)*(-Bus_Ymag(331+1)*sin(-theta(95+1)+theta(79+1)+Bus_Yang(331+1))*V(79+1)-Bus_Ymag(341+1)*sin(-theta(95+1)+theta(81+1)+Bus_Yang(341+1))*V(81+1)-Bus_Ymag(388+1)*sin(-theta(95+1)+theta(93+1)+Bus_Yang(388+1))*V(93+1)-Bus_Ymag(392+1)*sin(-theta(95+1)+theta(94+1)+Bus_Yang(392+1))*V(94+1)-Bus_Ymag(397+1)*sin(Bus_Yang(397+1))*V(95+1)+Bus_Ymag(400+1)*sin(theta(95+1)-theta(96+1)-Bus_Yang(400+1))*V(96+1)); % Depth 1, Step 255
  mismatch(255+1)=Q_load(96+1)+V(96+1)*(-Bus_Ymag(332+1)*sin(-theta(96+1)+theta(79+1)+Bus_Yang(332+1))*V(79+1)-Bus_Ymag(398+1)*sin(-theta(96+1)+theta(95+1)+Bus_Yang(398+1))*V(95+1)-Bus_Ymag(401+1)*sin(Bus_Yang(401+1))*V(96+1)); % Depth 1, Step 256
  mismatch(256+1)=Q_load(97+1)+V(97+1)*(Bus_Ymag(333+1)*sin(theta(97+1)-theta(79+1)-Bus_Yang(333+1))*V(79+1)-Bus_Ymag(403+1)*sin(Bus_Yang(403+1))*V(97+1)+Bus_Ymag(410+1)*sin(theta(97+1)-theta(99+1)-Bus_Yang(410+1))*V(99+1)); % Depth 1, Step 257
  mismatch(257+1)=Q_load(98+1)+V(98+1)*(Bus_Ymag(334+1)*sin(theta(98+1)-theta(79+1)-Bus_Yang(334+1))*V(79+1)-Bus_Ymag(406+1)*sin(Bus_Yang(406+1))*V(98+1)+Bus_Ymag(411+1)*sin(theta(98+1)-theta(99+1)-Bus_Yang(411+1))*V(99+1)); % Depth 1, Step 258
  mismatch(258+1)=Q_load(99+1)+V(99+1)*(-Bus_Ymag(379+1)*sin(-theta(99+1)+theta(91+1)+Bus_Yang(379+1))*V(91+1)+Bus_Ymag(389+1)*sin(theta(99+1)-theta(93+1)-Bus_Yang(389+1))*V(93+1)-Bus_Ymag(404+1)*sin(-theta(99+1)+theta(97+1)+Bus_Yang(404+1))*V(97+1)-Bus_Ymag(407+1)*sin(-theta(99+1)+theta(98+1)+Bus_Yang(407+1))*V(98+1)-Bus_Ymag(412+1)*sin(Bus_Yang(412+1))*V(99+1)+Bus_Ymag(417+1)*sin(theta(99+1)-theta(100+1)-Bus_Yang(417+1))*V(100+1)-Bus_Ymag(423+1)*sin(-theta(99+1)+theta(102+1)+Bus_Yang(423+1))*V(102+1)-Bus_Ymag(428+1)*sin(-theta(99+1)+theta(103+1)+Bus_Yang(428+1))*V(103+1)-Bus_Ymag(438+1)*sin(-theta(99+1)+theta(105+1)+Bus_Yang(438+1))*V(105+1))-5.5555555555555553581825733999722*V(99+1)*E(18+1)*sin(theta(99+1)-delta_t1(18+1)+1.5707963267948965579989817342721); % Depth 1, Step 259
  mismatch(259+1)=Q_load(100+1)+V(100+1)*(-Bus_Ymag(413+1)*sin(-theta(100+1)+theta(99+1)+Bus_Yang(413+1))*V(99+1)-Bus_Ymag(418+1)*sin(Bus_Yang(418+1))*V(100+1)-Bus_Ymag(421+1)*sin(-theta(100+1)+theta(101+1)+Bus_Yang(421+1))*V(101+1)); % Depth 1, Step 260
  mismatch(260+1)=Q_load(101+1)+V(101+1)*(Bus_Ymag(380+1)*sin(theta(101+1)-theta(91+1)-Bus_Yang(380+1))*V(91+1)+Bus_Ymag(419+1)*sin(theta(101+1)-theta(100+1)-Bus_Yang(419+1))*V(100+1)-Bus_Ymag(422+1)*sin(Bus_Yang(422+1))*V(101+1)); % Depth 1, Step 261
  mismatch(261+1)=Q_load(102+1)+V(102+1)*(Bus_Ymag(414+1)*sin(theta(102+1)-theta(99+1)-Bus_Yang(414+1))*V(99+1)-Bus_Ymag(424+1)*sin(Bus_Yang(424+1))*V(102+1)-Bus_Ymag(429+1)*sin(-theta(102+1)+theta(103+1)+Bus_Yang(429+1))*V(103+1)+Bus_Ymag(432+1)*sin(theta(102+1)-theta(104+1)-Bus_Yang(432+1))*V(104+1)+Bus_Ymag(451+1)*sin(theta(102+1)-theta(109+1)-Bus_Yang(451+1))*V(109+1)); % Depth 1, Step 262
  mismatch(262+1)=Q_load(103+1)+V(103+1)*(Bus_Ymag(415+1)*sin(theta(103+1)-theta(99+1)-Bus_Yang(415+1))*V(99+1)+Bus_Ymag(425+1)*sin(theta(103+1)-theta(102+1)-Bus_Yang(425+1))*V(102+1)-Bus_Ymag(430+1)*sin(Bus_Yang(430+1))*V(103+1)+Bus_Ymag(433+1)*sin(theta(103+1)-theta(104+1)-Bus_Yang(433+1))*V(104+1)); % Depth 1, Step 263
  mismatch(263+1)=Q_load(104+1)+V(104+1)*(-Bus_Ymag(426+1)*sin(-theta(104+1)+theta(102+1)+Bus_Yang(426+1))*V(102+1)-Bus_Ymag(431+1)*sin(-theta(104+1)+theta(103+1)+Bus_Yang(431+1))*V(103+1)-Bus_Ymag(434+1)*sin(Bus_Yang(434+1))*V(104+1)-Bus_Ymag(439+1)*sin(-theta(104+1)+theta(105+1)+Bus_Yang(439+1))*V(105+1)+Bus_Ymag(442+1)*sin(theta(104+1)-theta(106+1)-Bus_Yang(442+1))*V(106+1)+Bus_Ymag(445+1)*sin(theta(104+1)-theta(107+1)-Bus_Yang(445+1))*V(107+1)); % Depth 1, Step 264
  mismatch(264+1)=Q_load(105+1)+V(105+1)*(Bus_Ymag(416+1)*sin(theta(105+1)-theta(99+1)-Bus_Yang(416+1))*V(99+1)+Bus_Ymag(435+1)*sin(theta(105+1)-theta(104+1)-Bus_Yang(435+1))*V(104+1)-Bus_Ymag(440+1)*sin(Bus_Yang(440+1))*V(105+1)+Bus_Ymag(443+1)*sin(theta(105+1)-theta(106+1)-Bus_Yang(443+1))*V(106+1)); % Depth 1, Step 265
  mismatch(265+1)=Q_load(106+1)+V(106+1)*(-Bus_Ymag(436+1)*sin(-theta(106+1)+theta(104+1)+Bus_Yang(436+1))*V(104+1)-Bus_Ymag(441+1)*sin(-theta(106+1)+theta(105+1)+Bus_Yang(441+1))*V(105+1)-Bus_Ymag(444+1)*sin(Bus_Yang(444+1))*V(106+1)); % Depth 1, Step 266
  mismatch(266+1)=Q_load(107+1)+V(107+1)*(-Bus_Ymag(437+1)*sin(-theta(107+1)+theta(104+1)+Bus_Yang(437+1))*V(104+1)-Bus_Ymag(446+1)*sin(Bus_Yang(446+1))*V(107+1)+Bus_Ymag(448+1)*sin(theta(107+1)-theta(108+1)-Bus_Yang(448+1))*V(108+1)); % Depth 1, Step 267
  mismatch(267+1)=Q_load(108+1)+V(108+1)*(-Bus_Ymag(447+1)*sin(-theta(108+1)+theta(107+1)+Bus_Yang(447+1))*V(107+1)-Bus_Ymag(449+1)*sin(Bus_Yang(449+1))*V(108+1)+Bus_Ymag(452+1)*sin(theta(108+1)-theta(109+1)-Bus_Yang(452+1))*V(109+1)); % Depth 1, Step 268
  mismatch(268+1)=Q_load(109+1)+V(109+1)*(-Bus_Ymag(427+1)*sin(-theta(109+1)+theta(102+1)+Bus_Yang(427+1))*V(102+1)-Bus_Ymag(450+1)*sin(-theta(109+1)+theta(108+1)+Bus_Yang(450+1))*V(108+1)-Bus_Ymag(453+1)*sin(Bus_Yang(453+1))*V(109+1)+Bus_Ymag(456+1)*sin(theta(109+1)-theta(110+1)-Bus_Yang(456+1))*V(110+1)+Bus_Ymag(458+1)*sin(theta(109+1)-theta(111+1)-Bus_Yang(458+1))*V(111+1)); % Depth 1, Step 269
  mismatch(269+1)=Q_load(110+1)+V(110+1)*(-Bus_Ymag(454+1)*sin(-theta(110+1)+theta(109+1)+Bus_Yang(454+1))*V(109+1)-Bus_Ymag(457+1)*sin(Bus_Yang(457+1))*V(110+1)); % Depth 1, Step 270
  mismatch(270+1)=Q_load(111+1)+V(111+1)*(-Bus_Ymag(455+1)*sin(-theta(111+1)+theta(109+1)+Bus_Yang(455+1))*V(109+1)-Bus_Ymag(459+1)*sin(Bus_Yang(459+1))*V(111+1))+5.5555555555555553581825733999722*V(111+1)*E(19+1)*sin(-theta(111+1)+delta_t1(19+1)-1.5707963267948965579989817342721); % Depth 1, Step 271
  mismatch(271+1)=Q_load(112+1)+V(112+1)*(-Bus_Ymag(68+1)*sin(-theta(112+1)+theta(16+1)+Bus_Yang(68+1))*V(16+1)-Bus_Ymag(126+1)*sin(-theta(112+1)+theta(31+1)+Bus_Yang(126+1))*V(31+1)-Bus_Ymag(462+1)*sin(Bus_Yang(462+1))*V(112+1)); % Depth 1, Step 272
  mismatch(272+1)=Q_load(113+1)+V(113+1)*(-Bus_Ymag(127+1)*sin(-theta(113+1)+theta(31+1)+Bus_Yang(127+1))*V(31+1)-Bus_Ymag(464+1)*sin(Bus_Yang(464+1))*V(113+1)-Bus_Ymag(467+1)*sin(-theta(113+1)+theta(114+1)+Bus_Yang(467+1))*V(114+1)); % Depth 1, Step 273
  mismatch(273+1)=Q_load(114+1)+V(114+1)*(-Bus_Ymag(106+1)*sin(-theta(114+1)+theta(26+1)+Bus_Yang(106+1))*V(26+1)-Bus_Ymag(465+1)*sin(-theta(114+1)+theta(113+1)+Bus_Yang(465+1))*V(113+1)-Bus_Ymag(468+1)*sin(Bus_Yang(468+1))*V(114+1)); % Depth 1, Step 274
  mismatch(274+1)=Q_load(115+1)+V(115+1)*(-Bus_Ymag(279+1)*sin(-theta(115+1)+theta(67+1)+Bus_Yang(279+1))*V(67+1)-Bus_Ymag(470+1)*sin(Bus_Yang(470+1))*V(115+1)); % Depth 1, Step 275
  mismatch(275+1)=Q_load(116+1)+V(116+1)*(-Bus_Ymag(46+1)*sin(-theta(116+1)+theta(11+1)+Bus_Yang(46+1))*V(11+1)-Bus_Ymag(472+1)*sin(Bus_Yang(472+1))*V(116+1)); % Depth 1, Step 276
  mismatch(276+1)=Q_load(117+1)+V(117+1)*(-Bus_Ymag(310+1)*sin(-theta(117+1)+theta(74+1)+Bus_Yang(310+1))*V(74+1)-Bus_Ymag(313+1)*sin(-theta(117+1)+theta(75+1)+Bus_Yang(313+1))*V(75+1)-Bus_Ymag(475+1)*sin(Bus_Yang(475+1))*V(117+1)); % Depth 1, Step 277
  mismatch(277+1)=Q_load(118+1)+V(118+1)*(Bus_Ymag(173+1)*sin(theta(118+1)-theta(43+1)-Bus_Yang(173+1))*V(43+1)+Bus_Ymag(177+1)*sin(theta(118+1)-theta(44+1)-Bus_Yang(177+1))*V(44+1)-Bus_Ymag(478+1)*sin(Bus_Yang(478+1))*V(118+1)); % Depth 1, Step 278


