-- by Anton Zubets, Moscow Engineering Physics Institute
-- 3 spins of an atom
my_model3( x[3], a[1] )
{

array xt[3];
array center[3];

center = [5,25,15];
spin0 = hfTorusX(x,center,11,0.3);

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

tmp = hfShift3D(xt,5,25,0);
tmp = hfRotate3DZ(xt,0.5);

center = [0,0,15];
spin1 = hfTorusY(xt,center,11,0.3);

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

tmp = hfShift3D(xt,5,25,0);
tmp = hfRotate3DZ(xt,-0.5);

center = [0,0,15];
spin2 = hfTorusY(xt,center,11,0.3);

my_model3 = spin0 | spin1 | spin2;

}