% Session p. 259/260

h = 1/128; % N = 2559 x 255
p = poisson([0,0],[10,1],0,[1,1,0,0],'Cholesky',h)
p = poisson([0,0],[10,1],0,[1,1,0,0],'FFT',h)