function [arena] = ConcentrationRandomMultipleOdourSourcesNEW(X,Y)
arena = zeros(X, Y);

nos_min = 3; % minimum number of food sources
nos_max = 10; % maximum number of food sources
nos = randi([nos_min nos_max],1,1); % the actual number we are going to use

dist_min = 0; % minimum distance between food sourses where larvae can still sense odour
dist_max = 27; % maximum distance between food sourses where larvae can still sense odour

sources_arr = zeros(nos, 2); % array of positions of food sourses

sources_arr(1, 1) = -650+(1300*rand(1,1)); % random X-coordinate of the first food source
sources_arr(1, 2) = -1000+(2000*rand(1,1)); % random Y-coordinate of the first food source

dist = randi([dist_min dist_max],1,1); % random distance from the first food source to the second
angle_mem = 0; % angle for the direction of the second food source
rx = 0;
ry = 0;

while check ~= 1
    angle_mem = randi([0 359], 1, 1);
    rx = sources_arr(1, 1) + dist * sind(angle_mem); % hypothetical X-coordinate of the 2nd food source
    ry = sources_arr(1, 2) + dist * cosd(angle_mem); % hypothetical Y-coordinate of the 2nd food source
    if (rx > 0) and (rx < 650) and (ry > 0) and (ry < 1000)
        sources_arr(2, 1) = rx;
        sources_arr(2, 2) = ry;
        check = 1;
    end
end


dist_mem = dist; % distance between 2 previous food sources
angle = 0; % angle for the direction of the next food source
dist_p = 0; % distance from the (i-2)nd food source to the current one
dist_unimp = 0;
dist_unimp2 = 0;
angle_add = 0; % angle for completeting the angle of a triangle

% cycle for filling the array with random locations of food sources
i = 3; % because I already defined the position of the first 2 food sources

for i = 3 : nos
    while check ~= 1
        dist = randi([dist_min dist_max],1,1); % random distance from the current food source to the next one
        angle = randi([0 359],1,1);

        rx = sources_arr(i, 1) + dist * sind(angle); % hypothetical X-coordinate of the 2nd food source
        ry = sources_arr(i, 2) + dist * cosd(angle); % hypothetical Y-coordinate of the 2nd food source

        if (rx > 0) and (rx < 650) and (ry > 0) and (ry < 1000)

            dist_unimp = dist_mem * sind(angle_mem);
            dist_unimp2 = sqrt(dist_unimp^2 + dist_mem^2);
            angle_add = 90 - asind(dist_unimp2 / dist_mem);

            dist_p = sqrt(dist_mem^2 + dist^2 - 2 * dist_mem * dist * cosd(angle+angle_add)); % dist from the (i-2)nd to i'th food source (pythagorian)

            if (dist_p <= dist_max) % check the distance between (i-2)nd and i'th sources
                sources_arr(i, 1) = rx;
                sources_arr(i, 2) = ry;
                check = 1;
            end
        end
    end
end

all_dist = zeros(nos, X, Y); % distances to all the food sources from each point in arena
k = 1;
p = 1;
q = 1;

for k = 1 : nos
    for p = 1 : X
        for q = 1 : Y
            all_dist(k, p, q) = sqrt((p - (X + sources_arr(k, 1))^2) + (q - (Y + sources_arr(k, 2)))^2);
        end
    end
end

% cycle for filling the concentrations "arena" matrix

j = 1;
f = 1;
l = 1;

for j = 1 : X
    for f = 1 : Y
        for l = 1 : nos
            arena(j, f) = arena(j, f) + exp(-(all_dist(l, j, f)^2)/40000))
        end
    end
end

im = zeros(Y, X);
v = 1;
w = 1;

for v = 1 : X
    for w = 1 : Y
        im(w, v) = arena(v, w);
    end
end

imagesc([0 650], [0 1000], im);
hold on;