Calculate a matrix operation needed to construct Bayesian confidence intervals
dsms(rkpk.obj)
rkpk.obj
an object returned from calling dsidr
a matrix. See the corresponding Fortran subroutine.