function ind = realmutate(ind, domains, rate)%REALMUTATE Summary of this function goes here%   Detailed explanation goes here   % double rnd, delta1, delta2, mut_pow, deltaq;   % double y, yl, yu, val, xy;   % double eta_m = id_mu;   eta_m=20;   numVariables = size(domains,1);   if (isstruct(ind))       a = ind.parameter;   else       a = ind;   end   for j = 1:numVariables        if (rand() <= rate)             y = a(j);                        yl = domains(j,1);            yu = domains(j,2);            delta1 = (y - yl) / (yu - yl);            delta2 = (yu - y) / (yu - yl);            rnd = rand();            mut_pow = 1.0 / (eta_m + 1.0);            if (rnd <= 0.5)                 xy = 1.0 - delta1;                val = 2.0 * rnd + (1.0 - 2.0 * rnd) * (xy^(eta_m + 1.0));                deltaq = (val^mut_pow) - 1.0;            else                 xy = 1.0 - delta2;                val = 2.0 * (1.0 - rnd) + 2.0 * (rnd - 0.5) * (xy^ (eta_m + 1.0));                deltaq = 1.0 - (val^mut_pow);            end                        y = y + deltaq * (yu - yl);            if (y < yl)                y = yl;            end            if (y > yu)                y = yu;            end            a(j) = y;                end   end   if isstruct(ind)       ind.parameter = a;   else       ind = a;   endend