float square(float x, y, minx, maxx, miny, maxy, taper)
{
return  smooth(minx, minx + taper, x) *
                (1 - smooth(maxx, maxx + taper, x)) *
                smooth(miny, miny + taper, y) *
                (1 - smooth(maxy, maxy + taper, y));
}

displace
repeat
        (
        float height = 0.1,
        min_s = 0.1,
        max_s = 0.8,
        min_t = 0.1,
        max_t = 0.8,
        fuzz = 0.1,
        s_repeats = 3,
        t_repeats = 3)
{
float   hump = 0;
vector  n = normalize(N);

float ss = s * s_repeats%1;
float tt = t * t_repeats%1;
float row = floor(t * t_repeats);

if(row%2 == 0) ss = (s + 0.5) * s_repeats%1;

hump =  square(ss, tt, min_s, max_s, min_t, max_t, fuzz);

P = P - n * hump * height;
N = computenormal(P);

}