Embed Notice

vec3 col(in float x, in float y)
{
vec2 z = vec2(x,y);
int i = 0;
int n = 64;
for (i=0 ; i<n ; ++i)
{
vec2 z5 = cmul(z,cmul(z,cmul(z,cmul(z,z))));
vec2 a = (z5-vec2(1.0, 0.0));
vec2 b = (z5+vec2(1.0, 0.0));
vec2 c = (z5+vec2(pow(0.6625,5.0), 0.0));
z = cdiv(a,cdiv(b,cdiv(c,z)));
if (length(z) > 20.0) break;
}
float t = fract(float(i) / float(n)*16.0);
return hsv2rgb(vec3(0.1, 0.9, t));
}