function couple = TP1_3(I, x, y, d)
    couple=manhattan(I, x, y, d);
end

function couple = manhattan(I, x, y, d)
    couple=[[x y+d]; [x y-d]];
    for i=1:d
       couple=[couple; [x+i y+d-i]; [x-i y+d-i]; [x+i y-d+i]; [x-i y-d+i]];
    end
end

function manhattan2(d)
    I=zeros(d*2,d*2);
    for i=1:d
        
end