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

el1 = 1-(xt[1]/2)^2- (xt[2]/6)^2- (xt[3]/2)^2;

tmp = hfShift3D(xt,5.,0.,0.);
el2 = 1-(xt[1]/2)^2- (xt[2]/6)^2- (xt[3]/2)^2;

xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];
tmp = hfShift3D(xt,0.,4.,5.);
el3 = 1-(xt[1]/2)^2- (xt[2]/6)^2- (xt[3]/2)^2;

my_model = el1 | el2 | el3;
}