Question

I have this algorithm, but I am not too keen on the many if-statements.

Can someone see if this function can be written in a cleaner way?

rand('twister',101)

n = 10;
f = A.^(0:n)./factorial(0:n);
f = f/sum(f);
n = 10000;
Xi = 2;
X = zeros(1,n);

for i =1:n,
    dXi = (-1)^round(rand);
    Yi = Xi + dXi;
    if Yi > 0 & Yi <= length(f),
        if f(Yi) >= f(Xi),
            X(i) = Yi;
            Xi = Yi;
        else
            if rand <= f(Yi)/f(Xi),
                X(i) = Yi;
                Xi = Yi;
            else
                X(i) = Xi;
            end
        end
    end
    X(i) = Xi;
end

Answer 1

I don’t know Matlab syntax, but generally something like this:

if (cond1) then
   mainAction
else if (cond2) then
        mainAction
     else
        otherAction

can be simplified as:

if (cond1 OR cond2) then
   mainAction
else
   otherAction

The OR would have to be short-circuiting for exact equivalence, but if cond2 has no side-effects then it doesn’t really matter.