Max - please include all of your code and describe what the above is intended to do. It appears to be some sort of factoring (?) algorithm.
A recursive function for the above for loop could be
function [r,v] = getrv(q,qn,r,v)
if qn<=q+1
if r>q
vn = floor(r/q);
rn = 10*(r-vn*q);
else
vn = 0;
rn = 10*r;
end
[rn,vn] = getrv(q,qn+1,rn,vn);
r = [r rn];
v = [v vn];
end
and it would be called from your main function as
[r,v] = getrv(q,2,10*p,0);
Note how the recursive call to getrv is concatenated with the r and v inputs.
