my_model(x[3], a[1])
{
array xt[3];
xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];

xt[2]=xt[2]+6; -- shift down y
tmp = hfTwistZ(xt,-6,6,0,3.14);
el1 = 1-(xt[1]/3)^4- (xt[2]/3)^4- (xt[3]/8)^4; 

xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];
tmp = hfTwistY(xt,-6,6,0,-3.14);
el2 = 1-(xt[1]/3)^4- (xt[2]/6)^4- (xt[3]/4)^4; 

xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];

xt[2]=xt[2]-6; -- shift up y
tmp = hfTwistX(xt,-6,6,0,-3.14);
el3 = 1-(xt[1]/8)^4- (xt[2]/3)^4- (xt[3]/3)^4;

my_model = el1 | el2  | el3;
}