function f=hcurve(t);

x=1:1:12;
y=[-31.4920 -19.4311 -3.9216 9.9034 19.0078 22.5756 22.9987 20.3245 8.8065 ...
   -9.1112 -25.9921 -35.7874];
f=interp1(x,y,t,'cubic');
return