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)^4- (xt[2]/3)^4- (xt[3]/8)^4;

xt[1]=xt[1]+7;
tmp = hfTaperZ(xt,-3,3,1,0.5);
el2 = 1-(xt[1]/3.)^4- (xt[2]/2)^4- (xt[3]/8)^4; 

xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];
xt[3]=xt[3]-6;
xt[1]=xt[1]-6;
tmp = hfTaperY(xt,-4,2,1.2,0.4);
el3 = 1-(xt[1]/2)^4- (xt[2]/8)^4- (xt[3]/3)^4; 

xt[1]=x[1];
xt[2]=x[2];
xt[3]=x[3];
xt[1]=xt[1]-2;
xt[2]=xt[2]-6;
tmp = hfTaperX(xt,-2,2,0.4,1.1);
el4 = 1-(xt[1]/7)^4- (xt[2]/2)^4- (xt[3]/3)^4;

my_model = el1 | el2 | el3 | el4;
}