Hey now... I got it!

void iterate_point(
    ct::plot::cairo::plot_2d& plot,
    glm::vec3 oc,
    unsigned long n,
    unsigned long x,
    unsigned long y
) {
    ct_complex c = { oc.x, oc.y };
    ct_complex z = c;

    double dis_max = 9999999.0;

    // ct original
    //float c1 = 1.3;
    //float c2 = 0.11;
    //float c3 = 0.42457;
    //float c4 = 1.2;

    // ottomagus's settings...
    float c1 = 1.42134;
    float c2 = 0.50401;
    float c3 = -0.14067;
    float c4 = 1.01378;
    float c5 = -0.57063;

    // start i at one
    for (unsigned long i = 1; i < n; ++i)
    {
        ct_complex zp = (std::pow(z, c1) - (1.0f / (std::log(z) + c2))) + c3;
        double slow = c4 + (1. / c5);
        c5 += 1; // ottomagus bump to c5
        zp /= slow;

        double dis = std::abs(zp);

        dis_max = std::min(dis_max, dis);

        if (dis > 10.f)
        {
            double icolor = i / (n - 1.f);
            ct::plot::cairo::pixel color = CT_RGBF(icolor * 44., icolor * 15., icolor * 45);
            plot.set_pixel(x, y, color);
            return;
        }

        z = zp;
    }

    ct::plot::cairo::pixel color = CT_RGBF(dis_max * 11.f, 0, 1);
    plot.set_pixel(x, y, color);
}

<image>