coma movie
N=256; D=16; [p x] = make_pupil(N,D); [xg yg] = meshgrid(x); w = (xg).*((xg.^2+yg.^2)-1); z = psf(p,w); logim(z,3);
movie1('coma.avi',p,w);
N=256; D=16; [p x] = make_pupil(N,D); [xg yg] = meshgrid(x); w = (xg).*((xg.^2+yg.^2)-1); z = psf(p,w); logim(z,3);
movie1('coma.avi',p,w);