r=3+0*phi,vector(r,g,b)=vector(sin([20*x+b*cos([b*y+5*sin(10*a*z)])]),1,1),-0.75<sin([20*x+b*cos([b*y+5*sin(10*a*z)])])

a=slider([0,1])

b=slider([0,10])


Graph of the formula

The clipping function produces this fragment program:
!!ARBfp1.0

TEMP R0;
TEMP R1;
TEMP R2;
TEMP R3;

PARAM U0 = program.local[0];

PARAM C0 = {
0.500000000000000000000000000000,
0.750000000000000000000000000000,
20.000000000000000000000000000000,
5.000000000000000000000000000000
};
PARAM C1 = {
10.000000000000000000000000000000,
1.000000000000000000000000000000,
0.000000000000000000000000000000,
0.000000000000000000000000000000
};

MUL R0.x, C0.z, fragment.texcoord[0].x;
MUL R0.y, U0.x, fragment.texcoord[0].y;
MUL R0.z, U0.y, fragment.texcoord[0].z;
MUL R0.w, C1.x, R0.z;
SIN R1.x, R0.w;
MUL R1.y, C0.w, R1.x;
ADD R1.z, R0.y, R1.y;
COS R1.w, R1.z;
MUL R2.x, U0.x, R1.w;
ADD R2.y, R0.x, R2.x;
SIN R2.z, R2.y;
SLT R2.w, -C0.y, R2.z;
SLT R3.x, R2.w, C0.x;
KIL -R3.x;
MOV_SAT R0.x, R2.z;
MOV_SAT R0.y, C1.y;
MUL result.color.x, fragment.color.primary.x, R0.x;
MUL result.color.y, fragment.color.primary.y, R0.y;
MUL result.color.z, fragment.color.primary.z, R0.y;

END