%Hunt 1978 3-D dispersion solution Eqn.10.

clear all
close('all')
[x y z] = meshgrid(-1:0.05:3,-1:0.05:1,-1:0.05:1);
M3=1
Dxx=.001
Dyy=.0005
Dzz=.0005
sigma=.5
V=0.04


for t=1:10:51


data = M3*exp(-(x-V*t).^2/(4*Dxx*t)-y.^2/(4*Dyy*t)-z.^2/(4*Dzz*t))/(8*sigma*sqrt(pi^3*t^3*Dxx*Dyy*Dzz));
    
data(1,1,1)=30;%color scaling trick?


SLICE(x,y,z,data,[0],[0],[0])

box on

clear data

view(150,30); daspect([1 1 1]);axis([-1,3,-1,1,-1,1])
a=gcf;%get current figure
saveas(a,sprintf('%04d%s',t,'.bmp'),'bmp');
close

end