matlab 极坐标计算 代码实现

想请教一下，这种极坐标形式怎么用matlab代码实现，要随n变化而改变

你好，同学假设m=1, n=10

m = 1;
n = 10;
q = zeros(1,n)
for i = 1:n
if(i>=2)
q(1,i) = prod(cos(psi(1,i:n-1)))*sin(psi(1,i-1));
else
q(1,i) = prod(cos(psi(1,i:n-1)));
end
end