P must be a vector. S is required to compute error estimates. S matrix must be n+2-by-n where n = length(p). Zero degrees of freedom implies infinite error bounds.