Index in position 2 exceeds array bounds (must not exceed 3).

When calculating C (bolded) I get the following error "Index in position 2 exceeds array bounds (must not exceed 3)."

CXYZ=[166 64 64 0 0 0; 64 166 64 0 0 0; 64 64 166 0 0 0;0 0 0 80 0 0; 0 0 0 0 80 0; 0 0 0 0 0 80];
l(1,1) = 1;
l(1,2) = 0 ;
l(1,3) = 0;
l(2,1) = 0;
l(2,2) =1/sqrt(2);
l(2,3) =1/sqrt(2);
l(3,1) =0 ;
l(3,2)= -1/sqrt(2);
l(3,3) =1/sqrt(2);

g= [1, 2, 3, 4]
a=1
for (i = g)
    for (j= g)
        for(k=g)
            for (p=g)
C(a) = l(1,i).*l(1,j).*l(1,k).*l(1,p).*CXYZ(1,1);    % bolded
a= a+1;
            end
        end
    end
end

[per isakson]

Try to fix your code with the help of: Debug a MATLAB Program and Examine Values While Debugging

Answer 1

[Walter Roberson]

Because g includes 4, at some point you are going to try to access l(1,4) . However you only initialize l to (3,3) so l(4,:) and l(:,4) do not exist