def project_data(X, U, k):
    """Computes reduced data representation (projected data)

    Args:
        X (numpy.array): Normalized features' dataset
        U (numpy.array): eigenvectors of covariance matrix
        k (int): Number of features in reduced data representation

    Returns:
        numpy.array: Reduced data representation (projection)
    """
    U_reduce = U[:, 0:k]
    Z = X.dot(U_reduce)
    return Z