W1 = [0.999, 0.001];                                  %wlamanie
U1 = [0.998, 0.002];                                  %uderzenie pioruna
SA = [0.99, 0.05; 0.01, 0.95];                    %stefan i alarm
BA = [0.98, 0.35; 0.02, 0.65];                    %basia i alarm

%wlaczenie alarmu, wlamanie i uderzenie pioruna
AWU = zeros(2,2,2);
AWU(2, :, :) = [0.001, 0.29; 0.94, 0.95];   %alarm sie wlaczyl
AWU(1, :, :) = [0.999, 0.71; 0.06, 0.05];   %alarm sie nie wlaczyl

suma = 0;
result = zeros(2,2,2,2,2);
for A=1:2
    for B=1:2
        for S=1:2
            for U=1:2
                for W=1:2
                    result(A,B,S,U,W) = SA(S,A)*BA(B,A)*AWU(A,W,U)*W1(W)*U1(U);
                    suma = suma + result(A,B,S,U,W);
                end
            end
        end
    end
end