CALL RUPDT (rup, bup, sup, r, z <*>, b <*>, y <*>, ssq );
See the entry for the RDODT subroutine for details.