function(v.scale = 0.00412, gamma = 4.71, peak = 9)
{
#Program compute.upsilon()
#
#Default argumentes based on gamma function fitted to empirical data (Basia).
#
	x <- 1:(aa - 10) - 0.5
	h.scale <- peak/(gamma - 1)
	y <- v.scale * x^(gamma - 1) * exp( - x/h.scale)
	upsilon <- c(rep(0, times = 10), y)
	upsilon <- upsilon[1:aa]	#truncate if necessary
	names(upsilon) <- 1:aa - 0.5
	upsilon <- round(upsilon, 4)
	write.ascii(compute.upsilon)
	upsilon
}