function Y = CreateY(branchData)
if nargin < 1
branchData = InputData();
end
busf = branchData(:,1);
bust = branchData(:,2);
z = branchData(:,3);
branchCount = length(busf);
busCount = max([busf',bust']);
Y = zeros(busCount);
for n = 1:branchCount
Y(busf(n),bust(n)) = Y(busf(n),bust(n))-1/z(n);
Y(bust(n),busf(n)) = Y(busf(n),bust(n));
Y(busf(n),busf(n)) = Y(busf(n),busf(n))+1/z(n);
Y(bust(n),bust(n)) = Y(bust(n),bust(n))+1/z(n);
end
end
1