function S=test31(x)
// x0 = -10
// xo = 8.817
// S(xo) = -8.13

S=(x+10).*(x+6).*(x+5).*(x+1).*(x-7).*(x-10)*1e-4;

endfunction
