# Fraunhoferova difrakce na kruhovem otvoru
function difrakce2
  [x,y]=meshgrid(-50:1:50);
  pupila=double(sqrt(x.^2+y.^2)<5); # uvnitr kruhu 1 jinak 0
  obrazec=abs(fftshift(fft2(ifftshift(pupila))));
  mesh(x,y,obrazec);
  print -djpg difrakce2.jpg
endfunction

