my_model(x[3], a[1])
{
sphere = 81 - x[1]^2 - x[2]^2 - x[3]^2;
my_model = sphere + hfNoiseG(x, 1.8, 0.9, 1.4);
}