Trim zeros and nan's off of an N-D array in matlab

trimnd.m

function trimmed = trimnd( trimmed )
% Trim an N-D matrix with extra zero columns.

% Delete zero slices for a more compact viz
for ii = 1 : 3
    b = all( trimmed == 0 | isnan( trimmed ) ,ii);
    for jj = setdiff( 1 : 3, ii )
        deletecol = setdiff( 1 : 3, [ii,jj] );
        b2 = all( b, jj );
        for qq = 1 : 3; evalindex{qq} = ':';end
        evalindex{deletecol} = 'b2';
        evalexp = sprintf( '%s(%s,%s,%s)= [];','trimmed', ...
                            evalindex{1}, evalindex{2}, evalindex{3} );
        eval( evalexp );                        
    end
end