Thank you for your code. In the function of "function_GenerateFresnelPropagationStack", you give the code that "HStack = exp(1i*lambda*pi*(us.^2+vs.^2)* z/2);". I do't know why there must be "z/2". It looks like "z/2" should be "z" for the Fresnel propagation.