magandangles3(x,y,z,n)

function [g] = magandangles3(x,y,z,n)
    mag = sqrt(x.^2 + y.^2 + z.^2);
    %angle = atan2(y,x,z);
    
    maxmag = max(max(max(max(mag(:,:,:,1)))));
    
    g = figure;
    
    %subplot(1,2,1);
    k1 = imshow(mag(:,:,n,1), [0 maxmag]);
    colorbar;
    
%     %subplot(1,2,2);
%     k2 = imshow(angle(:,:,n,1), [0 2*pi]);
%     colorbar;
%     colormap jet;
    
    for i = 1:size(x,4)
        set(k1,'cdata',mag(:,:,n,i)); pause;
        %set(k2,'cdata',angle(:,:,n,i)); pause;
    end
end