%QUADFUNGTL
%
function y=quadfunGTL(x)
global s
if s==1
  y=3.2498+11.6854*(x-2)+4.1540*(x-2).*(x-3);
elseif s==2
  y=6.6971+26.0512*(x-2)+9.4751*(x-2).*(x-3);
elseif s==3
  y=11.3311+46.0945*(x-2)+16.9644*(x-2).*(x-3);
elseif s==4
  y=17.1463+71.8137*(x-2)+26.6224*(x-2).*(x-3);
elseif s==5
  y=24.1389+103.2081*(x-2)+38.4493*(x-2).*(x-3);
end
