[X, fs] = audioread('C:\Users\Dominik\Desktop\ISEL\SaU\Session1_SAu.wav');

a1 = 0.5;
a2 = 0.5;
a3 = 0.8;

L = 15;
Le1 = 5;
Le2 = 5;
Ls1 = 1;
Ls2 = 1;

headWidth = 0.18;

soundSpeed = 340.29;

%first refections
d1 = sqrt((abs(Le2-Ls2)*(Ls1/(Ls1+Le1)))^2 + Ls1^2);
d2 = sqrt((abs(Le2-Ls2)*(Le1/(Ls1+Le1)))^2 + Le1^2);
d = d1 + d2;
m1 = d/soundSpeed;
s1 = floor(m1*fs);

d1 = sqrt((abs(Le1-Ls1)*(Ls2/(Ls2+Le2)))^2 + (L-Ls2)^2);
d2 = sqrt((abs(Le1-Ls1)*(Le2/(Ls2+Le2)))^2 + (L-Le2)^2);
d = d1 + d2;
m2 = d/soundSpeed;
s2 = floor(m2*fs);

d1 = sqrt((abs(Le2-Ls2)*(Ls1/(Ls1+Le1)))^2 + Ls1^2);
d2 = sqrt((abs(Le2-Ls2)*(Le1/(Ls1+Le1)))^2 + Le1^2);
d = d1 + d2;
m3 = d/soundSpeed;
s3 = floor(m3*fs);

D = zeros([1, length(X)+max(s1, max(s2, s3))]);
D1 = zeros([1, length(X)+max(s1, max(s2, s3))]);
D2 = zeros([1, length(X)+max(s1, max(s2, s3))]);
D3 = zeros([1, length(X)+max(s1, max(s2, s3))]);
S = zeros([1, length(X)+max(s1, max(s2, s3))]);
for i = 1:length(X)
    D1(i+s1) = X(i)*(1-a1);
    D2(i+s2) = X(i)*(1-a2);
    D3(i+s3) = X(i)*(1-a3);
    D(i) = D1(i)+D2(i)+D3(i);
    S(i) = D(i) + X(i);
end