P=poly(S);

PP=poly2str(P,'s')

PA=polyval(P,S)

PM=polyvalm(P,S)    

S =

     1     1     1     1

     1     2     3     4

     1     3     6    10

     1     4    10    20

PP =

   s^4 - 29 s^3 + 72 s^2 - 29 s + 1

PA =

  1.0e+004 *

    0.0016    0.0016    0.0016    0.0016

    0.0016    0.0015   -0.0140   -0.0563

    0.0016   -0.0140   -0.2549   -1.2089