Фрактальный генерирующий код просто не работает

Я скопировал этот код из интернета, но, похоже, он не работает.

Все, что он делает, это рисует красную линию на одном краю черной картинки, которую делает.

Я играл с ним довольно долго и обнаружил, что когда я печатаю содержимое r1, r2 и r3, он пишет:
(0,0)
(0.866025,0)
(-0.866025,0)
Так что это кажется одной вещью, что он отображает части комплексного числа в неправильном порядке, но почему он отображает r1, чтобы быть нулем, когда это явно нет?

Кроме того, кажется ли это причиной того, что этот код не работает?

GLuint CMyApp::NewtonFractalTexture()
{
const int pic_size = 256;
unsigned char tex[pic_size][pic_size][3];

int MaxCount = 255;
int color_multiplier = 15;
float precision = 0.0001;
std::complex<double> r1 = (1, 0);
std::complex<double> r2 = (-0.5, sin(2 * PI / 3));
std::complex<double> r3 = (-0.5, -sin(2 * PI / 3));

std::cout << r1 << " " << r2 << " " << r3 << std::endl;
std::cout << abs(r1) << " " << abs(r2) << " " << abs(r3) << std::endl;

/*
std::complex<double> roots[birds_num];
for (int i = 0; i < birds_num; ++i){
roots[i] = (bird_positions[i][0], bird_positions[i][2]);
}
*/

for (int i = 0; i < pic_size; ++i){
for (int j = 0; j < pic_size; ++j)
{
//
std::complex<double> z = (i, j);
//
int count = 0;
while (count < MaxCount && abs(z - r1) >= precision && abs(z - r2) >= precision && abs(z - r3) >= precision){
/*
std::complex<double> my_numerator = (z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]);
std::complex<double> my_denominator = (z - roots[0])*(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[6])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[7])*(z - roots[8]) +
(z - roots[0])*(z - roots[1])*(z - roots[2])*(z - roots[3])*(z - roots[4])*(z - roots[5])*(z - roots[6])*(z - roots[7])*(z - roots[8])*(z - roots[8]);
*/
std::complex<double> my_numerator = z*z*z - 1.0;
std::complex<double> my_denominator = z*z * 3.0;
if (abs(z) > 0){
z = z - my_numerator / my_denominator;
}
++count;
}
//
tex[i][j][0] = 0;
tex[i][j][1] = 0;
tex[i][j][2] = 0;
//
if (abs(z - r1) < precision){
tex[i][j][0] = 255 - count * color_multiplier;
}
if (abs(z - r2) <= precision){
tex[i][j][1] = 255 - count * color_multiplier;
}
if (abs(z - r3) <= precision){
tex[i][j][2] = 255 - count * color_multiplier;
}
//
}
}

GLuint tmpID;

// generáljunk egy textúra erőforrás nevet
glGenTextures(1, &tmpID);
// aktiváljuk a most generált nevű textúrát
glBindTexture(GL_TEXTURE_2D, tmpID);
// töltsük fel adatokkal az...
gluBuild2DMipmaps(  GL_TEXTURE_2D,          // aktív 2D textúrát
GL_RGB8,                // a vörös, zöld és kék csatornákat 8-8 biten tárolja a textúra
pic_size, pic_size,     // kép méretének megadása
GL_RGB,                 // a textúra forrása RGB értékeket tárol, ilyen sorrendben
GL_UNSIGNED_BYTE,       // egy-egy színkopmonenst egy unsigned byte-ról kell olvasni
tex);                   // és a textúra adatait a rendszermemória ezen szegletéből töltsük fel
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);   // bilineáris szűrés kicsinyítéskor
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);   // és nagyításkor is
glBindTexture(GL_TEXTURE_2D, 0);

return tmpID;
}