function [X] = compute_ikeda_trajectory(u, x, y, N)
    X = zeros(N,2);
    
    X(1,:) = [x y];
    for n = 2:N
        t = 0.4 - 6/(1 + x^2 + y^2);
        x1 = 1 + u*(x*cos(t) - y*sin(t)) ;
        y1 = u*(x*sin(t) + y*cos(t)) ;
        x = x1;
        y = y1;
        X(n,:) = [x y];
    end
end