Thanks a lot, I found the answer in source code. ‘coefficients’ function gives only constant parts of piecewise polynomials. If linear, quadratic and cubic parts of polynomials are needed to get, ‘itp.m’, ‘itp.c’, ‘itp.d’ should be written, respectively.