my_model(x[3], a[1]){
array vect[18];
array S[5];
vect = [-8.5, -8.5, 5.0,
         8.0, -9.0, 5.0,
         8.0, -4.0, 5.0,
        -5.0, 7.0, 5.0,
         8.0, 4.0, 8.0,
         8.5, 8.5, 8.5];
S = [1.05, 0.85, 0.8, 0.7, 0.63];

mycurve = hfConvCurve(x,vect,S,0.7);
my_model = mycurve;
}