function dS=gtest20(x)
// Papageorgiou, 1996
//
// x0 = 2
// xo = 0.1429
// S(xo) = [1.0102 1.0102]'

// S=4*(x-0.5).^2+0.5;

dS=8*x - 4;

endfunction
