my_model(x[3], a[1]){

e = 3;
sp1 = 1 -((x[1]-5)/e)^4-((x[2]+3)/e)^4-(x[3]/e)^4;
sp2 = 1-((x[1]+3)/e)^4-((x[2]+3)/e)^4-(x[3]/e)^4;
blend1 = hfBlendUni(sp1,sp2,5,2,2);

sp3=16-(x[1]-6)^2-(x[2]-6)^2-(x[3]-6)^2;
blend2=hfBlendInt(sp3,6-x[1],-30,1,1);
my_model = blend1 | blend2;
}