double precision s(mxpart,mxpart) common/sprods/s