function img_target = diffracted_image_at_target(img_source, xpos_s, ypos_s, image_target, xpos_t, ypos_t, z_t, lambda)

[Ny,Nx]=size(image_target);

xim_t=linspace(xpos_t(1), xpos_t(end), Nx);
yim_t=linspace(ypos_t(1), ypos_t(end), Ny);

img_target=zeros(Ny, Nx);

for i=1:Nx
    for k=1:Ny
        intensity = beam_intensity_at_point_from_image(xim_t(i), yim_t(k), z_t, img_source, xpos_s, ypos_s, lambda);
        img_target(k,i) = intensity;
    end
end