# some S/Splus  code to generate a random site-swap juggling pattern
# created April, 1995.

nballs <- 4
nsites <- 4


ss <- function(nballs,nsites)
	{
	ok <- 0	
	while( ok == 0 )
	 {
		perm <- sample(1:nsites)
		foo<-sort(c(0,sample(0:nballs,replace=T,
			size=(nsites-1) ),nballs) )
		balls <- foo[2:(nsites+1)] - foo[1:nsites]
		siteswap <- perm - (1:nsites) + nsites*balls
		ok <- ifelse( any(siteswap<0), 0, 1 )
	 }
	return( siteswap )
	}