function S=test0(x)
// x0 = 2
// xo = 0.5
// S(xo) = 0.5

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

endfunction
