BE2A = 2*0.3048;
BE2B = 2*0.3048;
BETA = 0;
for i = 1:50
        for j = 1:i
            for k = 1:(i+1-j)
                % 1st quadrant
                BEXC(i,j,k,1) = i*0+j*0+(0.5+(k-1))*BE2A*cos(BETA/180*pi);
                BEYC(i,j,k,1) = i*0+j*0+(0.5+(k-1))*BE2A*sin(BETA/180*pi);
                BEZC(i,j,k,1) = i*0+(0.5+(j-1))*BE2B+k*0;
                BE2A(i,j,k,1) = BE2A+i*0+j*0+k*0;
                BE2B(i,j,k,1) = BE2B+i*0+j*0+k*0;
                ANGLE(i,j,k,1) = BETA+i*0+j*0+k*0;
                % 2nd quadrant
                BEXC(i,j,k,2) = i*0+j*0-(0.5+(k-1))*BE2A*cos(BETA/180*pi);
                BEYC(i,j,k,2) = i*0+j*0-(0.5+(k-1))*BE2A*sin(BETA/180*pi);
                BEZC(i,j,k,2) = i*0+(0.5+(j-1))*BE2B+k*0;         
                BE2A(i,j,k,2) = BE2A+i*0+j*0+k*0;
                BE2B(i,j,k,2) = BE2B+i*0+j*0+k*0;
                ANGLE(i,j,k,2) = BETA+i*0+j*0+k*0; 
            end
        end
    end
When I run this code, it appears that: "Assignment has more non-singleton rhs dimensions than non-singleton subscripts". The errors start from the line
- "BEXC(i,j,k,1) = i*0+j*0+(0.5+(k-1))*BE2A*cos(BETA/180*pi)"
 
     
    