Espiral arquimediana em C ++

Estou tentando traçar as posições x e y de uma espiral arquimediana em C ++.

Até agora eu tenho tentado algo assim, mas sem sorte:

int dx = 0;
int dy = 0;
int x = 0;
int y = 0;

for (int i = 0; i < maxPoints; i++)
{
    dx = sin(i * PI / 2);
    dy = cos(-i * PI / 2);
    x += dx;
    y += dy;

    plot(x, y);    
}

EDIT: Mais informações

Estou desenvolvendo um aplicativo de jogo em 3D que demonstra o uso do mecanismo de física Bullet simulando dominós. Em vez de colocar os dominós na cena manualmente, quero usar um pouco de matemática para fazer isso por mim :)

Para quem estiver interessado aqui, está no GitHub .

Descobri isso :) Os dominós agora estão sendo colocados ao longo das coordenadas X e Y geradas pela função.

O código original da pergunta estava traçando uma onda de pontos para fora da posição ou origem central e não era o que eu queria. O que eu precisava era que cada ponto seguisse Archimedean spiralcom um certo espaço entre as espirais.

Inicialmente, usei integervalores para armazenar as coordenadas xe, ymas isso estava causando um erro de precisão ao truncar o floating pointvalor para armazená-lo no integertipo de dados.

O exemplo abaixo gera pontos ao longo da espiral continuamente, em relação ao maxPointsvalor.

float x = 0;
float y = 0;
float angle = 0.0f;

// Space between the spirals
int a = 2, b = 2;

for (int i = 0; i < maxPoints; i++)
{
    angle = 0.1 * i;
    x = (a + b * angle) * cos(angle);
    y = (a + b * angle) * sin(angle);

    plot(x, y);
}

O código do projeto está no GitHub , você precisará de Bullet e freeglut

Essa não é realmente uma resposta direta a essa pergunta (que já tem uma resposta de qualquer maneira), mas pode interessar as pessoas que desejam implementar esse algoritmo em 3D.

Eu tive que tentar implementar esse algoritmo para gerar espirais 3D no blender usando Python (poderia ser facilmente convertido em desenho com PIL ou Matplotlib em 2D). Então, aqui está o algoritmo e o resultado:

import bpy
from math import cos, sin
S = bpy.context.scene

def add_archimedian_spiral( size = 0.1, length = 500, height = 1, name = 'archispiral' ):
    mesh = bpy.data.meshes.new( name = name )

    o = bpy.data.objects.new(name, mesh)
    o.location = (0,0,0) # place at object origin
    S.objects.link( o )

    z     = 0
    verts = []  
    for i in range( length ):
        angle = 0.1 * i
        x     = ( 2 * size * angle ) * cos( angle )
        y     = ( 2 * size * angle ) * sin( angle )
        z    += i / 10000 * height
        verts.append((x,y,z))

    edges = []
    for i in range( len( verts ) ):
        if i == len( verts ) - 1: break
        edges.append((i, i+1))

    mesh.from_pydata( verts, edges, [] )

add_archimedian_spiral( size = 0.2, length = 500, height = 6 )