alpha           = 0;
beta            = 0.3;
gamma           = sqrt(1 - alpha*alpha - beta*beta);

wt_vec          = linspace(0,2*pi,10);

xvec            = linspace(-2,2,50);
yvec            = linspace(-2,2,50);
zvec            = linspace(-2,2,50);

[X Y Z]         = meshgrid(xvec,yvec,zvec);

Ebase           = exp(j*2*pi*(X*alpha+Y*beta+Z*gamma));

for wt = wt_vec
 
    figure(1)
    Ereal       = real(Ebase*exp(-j*wt));
    xslice      = max(xvec);
    yslice      = [0 max(yvec)];
    zslice      = min(zvec);
    
    slice(X,Y,Z,Ereal,xslice,yslice,zslice)
    colormap jet;shading interp;%colorbar
    xlabel('X axis');ylabel('Yaxis');zlabel('Z Axis')
    axis('equal')
    caxis([-1 1])
   
    %view([0 0])
    if wt ==0
        disp('use enter key to start')
        colorbar
        pause
    else
        pause(1)
    end
    
end


disp('done')