Function D9upak

Part of:

package slatec
( d9upak < x > < y > < n > )