Nested Summation Computation for M_K with MATLAB

[omer21]

<br /> M_{K}=\frac{1}{2^{k+1}-2}\sum_{i=0}^{L-1}\sum_{l=1}^{K}\binom{K}{l}h_{i}i^{l}M_{K-l}<br />

M_0=1 and the size of h_i is L.

I tried to compute this summation in matlab, my attempt is as following:

clear
h=[ (1+sqrt(3))/4 (3+sqrt(3))/4 (3-sqrt(3))/4 (1-sqrt(3))/4]'; 
% for simplicity i take K=4 and L=4 
K=4;
L=4;
k=1;
M=zeros(1,K+1);
M(1)=1;

for l=1:K
    for i=1:L
        for j=1:k
            M(l+1)=M(l+1)+nchoosek(k,j)*h(i)*((i-1)^l)*M(l);
        end
    end
    k=k+1;
    M(l+1)=(1/(2^(K+1)-2))*M(l+1);
end

But I'm not sure whether the codes are correct or not. If the codes are not correct how can i fix them?

 

[omer21]

problem is solved.