in src/Algorithm/fractals.c [392:405]
void powN(float p, inout vec3 z, float zr0, inout float dr)
{
float zo0 = asin(z.z / zr0);
float zi0 = atan(z.y, z.x);
float zr = pow(zr0, p - 1.0);
float zo = zo0 * p;
float zi = zi0 * p;
float czo = cos(zo);
dr = zr * dr * p + 1.0;
zr *= zr0;
z = zr * vec3(czo * cos(zi), czo * sin(zi), sin(zo));
}