Maximum recursion limit of 500 reached. Use set(0,'Rec​ursionLimi​t',N) to change the limit. Be aware that exceeding your available stack space can crash MATLAB and/or your computer.

Hi, please how i can solve this problem:

Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N) to change the limit. Be aware that exceeding your
available stack space can crash MATLAB and/or your computer.

My code:

function [Pmacro Pfemto] = UAQP_power_control_method(handles)
Pmacro_dbm = str2double(get(handles.edit1,'String'));
Pmacro = (10.^(Pmacro_dbm / 10)) / 1000;
Pfemto_max_dbm = str2double(get(handles.edit2,'String'));
Pfemto_max = (10.^(Pfemto_max_dbm / 10)) / 1000;
Pfemto_def_dbm = str2double(get(handles.edit4,'String'));
Pfemto_def = (10.^(Pfemto_def_dbm / 10)) / 1000;
target_SINR = str2double(get(handles.edit5,'String'));
coordVectors = get(handles.figure1, 'UserData');
numOfFemtocells = str2num(get(handles.edit20,'String'));
Pfemto_init = ones(1, numOfFemtocells)*Pfemto_def;
pow_default = str2num(get(handles.edit13,'String'));
coordVectors.powFemto = (ones(1, numOfFemtocells))*(10^(pow_default / 10)) / 1000;
Pfemto_last = coordVectors.powFemto;
Pfemto=zeros(1,numOfFemtocells);
N21 = str2double(get(handles.edit15,'String'));
N31 = str2double(get(handles.edit35,'String'));
N1=str2num(get(handles.edit47,'String'));
N2=str2num(get(handles.edit18,'String'));
N3=str2num(get(handles.edit19,'String'));
coordVectors = get(handles.figure1, 'UserData');
numOfFemtousers = handles.edit11;
numOfMacrocells = 1;
P0=10;
o2=0.4;
o1=0.9;
o3=0.1;
Q21=0.7;
Q31=0.3;
Q22=0.6;
Q32=0.2;
D=0.7;
Pfemto = zeros(1,numOfFemtocells);
a = o2*Q21*ones(N2,1);
b = o2*Q22*ones(N2,1 );
c = o2*Q31*ones(N3,1 );
d = o2*Q32*ones(N3,1);
SINR1=zeros(1,N1);
SINR2=zeros(1,N2);
SINR3=zeros(1,N3);
R = str2num(get(handles.edit11,'String'));
currentSINR=zeros(1,numOfFemtocells);
target_SINR = str2double(get(handles.edit5,'String'));
coordVectors = get(handles.figure1, 'UserData');
P1=Pfemto;
P2=P0+(N1*o1)+sum(a(1:N2))+(sum(c(1:N3))*D);
for i=1:numOfFemtocells
    [row col]=size(i);
    SINRp11= zeros(row,col);
    SINRp12=zeros(row,col);
    SINRp21 = zeros(row,col);
    SINRp22 = zeros(row,col);
    SINRp31 = zeros(row,col);
    SINRp32=zeros(row,col);
    SINR11=zeros(row,col);
    SINR12=zeros(row,col);
    SINR21=zeros(row,col);
    SINR22=zeros(row,col);
    SINR31=zeros(row,col);
    SINR32=zeros(row,col);
    for j=1:N1
        [row col] = size(j);
        x1 = coordVectors.femt_x(j);
        y1 = coordVectors.femt_y(j);
        x = [x1 x1 max(x1-R,0) min(x1+R,496)];
        y = [min(y1+R,428) max(y1-R,0) y1 y1];

          if (Pfemto == P1)
              SINR11(j)=SINRfemto(x1, y1,j, handles);
              %SINR11(j)=findRangeSINR(handles);
          else
              SINR12(j) = SINRfemto(x1, y1,j, handles);
              %SINR12(j)=findRangeSINR(handles);
          end
      end
      SINRp11=sum( SINR11(j)/N1);
      SINRp12=sum(SINR12(j)/N1);
      for k=1:N2
          if (Pfemto == P1);
              %set(handles.figure1, 'UserData', coordVectors);
              SINR21(k)=SINRfemto(x, y,k, handles);
          else
              %coordVectors.powFemto = P2;
              % set(handles.figure1, 'UserData', coordVectors);
              SINR22(k) =SINRfemto(x, y, k, handles);
          end
          SINRp21=sum(SINR21(k)/N2);
          SINRp22=sum(SINR22(k)/N2);
      end
      for l=1:N3
          if (Pfemto == P1)
              %set(handles.figure1, 'UserData', coordVectors);
              SINR31(l)=SINRfemto(x, y, l, handles);
          else
              %set(handles.figure1, 'UserData', coordVectors);
              SINR32(l) = SINRfemto(x, y, l, handles);
          end
      end
      SINRp31=sum(SINR31(l)/N3);
      SINRp32=sum(SINR32(l)/N3);
  end
  for i=1:numOfFemtocells
      if (SINRp11> SINRp12 && SINRp21> SINRp22&& SINRp31> SINRp32)
          Pfemto=P2;
      else if (SINRp11> SINRp12 && SINRp21> SINRp22);
              Pfemto(i)=P0+(N1*o1)+sum(a(1:N2))*D;
          else if (SINRp11> SINRp12)
                  Pfemto(i)=P0+(N1*o1)+sum(a(1:N21))*D;
              else
                  Pfemto=P1;
              end
          end
      end
  end
  end

thanks.

[Stephen23]

@brahmi ibtissem: your function does not appear to call itself from within, so this does not explain the recursion. How are you calling this function? It would help to use the debugger and have a look at which functions have been called recursively.

[brahmi ibtissem]

thanks. This is the code of my fuction.My object is to calling this function in my code in order to calculate the SINR of each femto.

function [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto(var1, var2, var3, var4)
nargchk(2, 4, nargin);
if nargin == 2
    [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto1(var1, var2);
else
    [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto2(var1, var2, var3, var4);
end

% function [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto1(userID, handles)
% coordVectors = get(handles.figure1, 'UserData');
% cellID = coordVectors.att_cell(userID);
% [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto2(coordVectors.fuser_x(userID), coordVectors.fuser_y(userID), cellID, handles);

function [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto2(x, y, cellID, handles)
PSD_dBm = -174; %dBm/Hz
awgn_N0 = (10^(PSD_dBm / 10)) / 1000;
spacing_Df = 15*10^3;
[row col] = size(x);
upper_thres = ones(row, col)*(100);

coordVectors = get(handles.figure1, 'UserData');
femto_fr_mode = get(handles.popupmenu7,'Value');
macro_fr_mode = get(handles.popupmenu9,'Value');

numOfFemtocells = min(size(coordVectors.femt_x, 2), size(coordVectors.femt_y, 2));
numOfMacrocells = 1;

st = dbstack;

if coordVectors.gain_calculated_f(cellID) == 0 || coordVectors.gain_calculated_f(cellID) == 1 || coordVectors.gain_calculated_f(cellID) == 3  || strcmp(st(3).file,'findRangeSINR.m')
    [G, PL, numOfWalls, dist, isIndoor] = Gain_femto(x, y, coordVectors.femt_x(cellID), coordVectors.femt_y(cellID), handles);
    if (~strcmp(st(3).file,'findRangeSINR.m'))
        coordVectors = get(handles.figure1, 'UserData');
        coordVectors.(genvarname(num2str(cellID))) = G;
        coordVectors.(genvarname(num2str(cellID+numOfFemtocells))) = PL;
        coordVectors.(genvarname(num2str(cellID+ 2*numOfFemtocells))) = numOfWalls;
        coordVectors.(genvarname(num2str(cellID+ 3*numOfFemtocells))) = dist;
        coordVectors.(genvarname(num2str(cellID+ 4*numOfFemtocells))) = isIndoor;

          coordVectors.gain_calculated_f(cellID) = 2;
          set(handles.figure1, 'UserData', coordVectors);
      end
  else

      G = coordVectors.(genvarname(num2str(cellID)));
      PL = coordVectors.(genvarname(num2str(cellID+numOfFemtocells)));
      numOfWalls = coordVectors.(genvarname(num2str(cellID+ 2*numOfFemtocells)));
      dist = coordVectors.(genvarname(num2str(cellID+ 3*numOfFemtocells)));
      isIndoor = coordVectors.(genvarname(num2str(cellID+ 4*numOfFemtocells)));
  end

if strcmp(st(3).file,'findRangeSINR.m')
        pow_default = str2num(get(handles.edit13,'String'));

      coordVectors.powFemto = (ones(1, numOfFemtocells))*(10^(pow_default / 10)) / 1000;
      Pfemto = coordVectors.powFemto;

      coordVectors.powFemto = (ones(1, numOfFemtocells))*(10^(pow_default / 10)) / 1000;
          Pmacro_dbm = str2num(get(handles.edit1,'String'));
      Pmacro = (10^(Pmacro_dbm / 10)) / 1000;
  else
      [Pmacro Pfemto] = power_scheme(handles);
  end

numer = Pfemto(cellID) * G;

sum1 = zeros(row,col);
sum2 = zeros(row,col);

coordVectors = get(handles.figure1, 'UserData');
for i=1:numOfFemtocells
    if i ~= cellID
        if coordVectors.gain_calculated_f(i) == 0 || coordVectors.gain_calculated_f(i) == 3 || strcmp(st(3).file,'findRangeSINR.m')
            G_f = Gain_femto(x, y, coordVectors.femt_x(i), coordVectors.femt_y(i), handles);
            if (~strcmp(st(3).file,'findRangeSINR.m'))
                coordVectors = get(handles.figure1, 'UserData');
                  coordVectors.(genvarname(num2str(i))) = G_f;
                  coordVectors.gain_calculated_f(i)=1;
                  set(handles.figure1, 'UserData', coordVectors);
            end
        else
            G_f = coordVectors.(genvarname(num2str(i)));
        end

          if femto_fr_mode == 3 || femto_fr_mode == 4
              sum1 = sum1 + 0;
          else
              sum1 = sum1 + Pfemto(i) * G_f;
          end

      end
  end

for i=1:numOfMacrocells
    if coordVectors.gain_calculated_m == 0 || strcmp(st(3).file,'findRangeSINR.m')
        G_m = Gain_macro(x, y, 250.0, 217.0, handles);
        if (~strcmp(st(3).file,'findRangeSINR.m'))
                    coordVectors = get(handles.figure1, 'UserData');
                    coordVectors.Gain_m = G_m;
                    coordVectors.gain_calculated_m = 1;
                    set(handles.figure1, 'UserData', coordVectors);
        end
    else
        G_m = coordVectors.Gain_m;
    end

      sumall = sum(sum(Pmacro * G_m))
      if femto_fr_mode == 2 ||femto_fr_mode == 3
          sum2 = sum2 + 0;
      else
          sum2 = sum2 + Pmacro * G_m;
      end
          %sum2 = sum2 + Pmacro * G_m;
  end

sum2 = interference_calc_macro(x, y, sum2, Pmacro, handles, 2);

denom = awgn_N0 * spacing_Df + sum1 + sum2;
temp = numer./denom;
if strcmp(st(3).file,'findRangeSINR.m')
    SINR = temp;
else
    SINR = min(temp, upper_thres);
end
%coordVectors = get(handles.figure1, 'UserData');
%coordVectors.SINR = SINR;

[Stephen23]

@brahmi ibtissem: and what is the function stack when you set a dbstop if error ? Which function/s are being called recursively?

[brahmi ibtissem]

ST = dbstack(___) returns the stack trace information in an m-by-1 structure, ST. dbstack displays the line numbers and file names of the function calls that led to the current pause condition, listed in the order in which they execute. The display starts with the currently executing functions and continues until it reaches the topmost function. Each line number is a hyperlink to that line in the Editor. The notation functionname>localfunctionname describes the location of a local function.

[brahmi ibtissem]

Hi, the fuction is findRangeSINR and this is the code of this function:

function currentSINR = findRangeSINR(handles)

R = str2num(get(handles.edit11,'String'));
coordVectors = get(handles.figure1, 'UserData');

numOfFemtocells = min(size(coordVectors.femt_x, 2), size(coordVectors.femt_y, 2));
currentSINR=zeros(1,numOfFemtocells);

for i=1:numOfFemtocells
    x1 = coordVectors.femt_x(i);
    y1 = coordVectors.femt_y(i);

      x = [x1 x1 max(x1-R,0) min(x1+R,496)];
      y = [min(y1+R,428) max(y1-R,0) y1 y1];
      [SINR, G, PL, numOfWalls, dist, isIndoor] = SINRfemto(x, y, i, handles);

      if histc(coordVectors.att_cell, i)==0
          NoF=1;
      else
          NoF = histc(coordVectors.att_cell, i);
      end
      currentSINR(i) = min(SINR)/NoF;
  end

[Stephen23]

@brahmi ibtissem: I still don't see anything obvious that would cause recursion of those functions. You should check if the filenames match the function names: this can cause recursion, and users can easily miss this. Double check the filenames against the function names!

If you still cannot find anything please upload your files in a new comment.