for i = 1:7
    Y(i) = input(['Podaj wartosc Y ', num2str(i),': ']);
end
for i = 1:4
    E(i) = input(['Podaj wartosc E ', num2str(i),': ']);
end
I(1) = input('Podaj wartosc I1: ');
%Y(1) = 1;Y(2) = 2;Y(3) = 3;Y(4) = 4;Y(5) = 5;Y(6) = 6;Y(7) = 7;E(1) = 10;E(2) = 20;E(3) = 30;E(4) = 40;I(1) = 50;
Y
E
I
A = [
    (Y(1) + Y(2) + Y(7)), (-Y(2)), (-Y(1));
    (-Y(2)), (Y(2) + Y(3) + Y(4)), (-Y(3));
    (-Y(1)), (-Y(3)), (Y(1) + Y(3) + Y(5) + Y(6));
    ]
B = [
    ((-E(1) .* Y(1)) - (E(2) .* Y(7)));
    I(1);
    ((-E(4) .* Y(6)) + (E(3) .* Y(5)) + (E(1) .* Y(1)));
    ]
Ainv = inv(A)
V = Ainv * B
A * V
I1 = ((V(1) - V(3) + E(1)) * Y(1)); 
I2 = ((V(2) - V(1)) * Y(2));
I3 = ((V(2) - V(3)) * Y(3));
I4 = V(2) * Y(4);
I5 = ((-V(3) + E(3)) * Y(5));
I6 = ((V(3) + E(4)) * Y(6));
I7 = ((V(1) + E(2)) * Y(7));

IV1 = [
    -I1;
    I2;
    - I7;
    ]
sum(IV1)

IV2 = [
    -I2;
    -I3;
    I(1);
    -I4;
    ]
sum(IV2)

IV3 = [
    -I6;
    I1;
    I3;
    I5;
    ]
sum(IV3)