st.andrews <- function(x)
{
	mmax <- 20
	n <- nrow(x)
	m <- ncol(x)
	u <- seq( - pi, pi,  , 100)
	y <- seq(0, 0,  , mmax)
	zmax <- 2*m
	for(k in 1:m) {
		y[k] <- x[1, k]
	}
	plot(u, y[1]/sqrt(2) + y[2] * sin(u) + y[3] * cos(u) + y[4] * sin(2 * u) + y[5] * cos(2 * u) + y[6] * 
		sin(3 * u) + y[7] * cos(3 * u) + y[8] * sin(4 * u) + y[9] * cos(4 * u) + y[10] * sin(5 * u) + y[
		11] * cos(5 * u) + y[12] * sin(6 * u) + y[13] * cos(6 * u) + y[14] * sin(7 * u) + y[15] * cos(7 *
		u) + y[16] * sin(8 * u) + y[17] * cos(8 * u) + y[18] * sin(9 * u) + y[19] * cos(9 * u) + y[20] * 
		sin(10 * u), type = "l", xlab = "", ylab = "", ylim = c( - zmax, zmax))
	title("St. Andrews Curves")
	for(i in 2:n) {
		for(k in 1:m) {
			y[k] <- x[i, k]
		}
		lines(u, y[1]/sqrt(2) + y[2] * sin(u) + y[3] * cos(u) + y[4] * sin(2 * u) + y[5] * cos(2 * u) + 
			y[6] * sin(3 * u) + y[7] * cos(3 * u) + y[8] * sin(4 * u) + y[9] * cos(4 * u) + y[10] * 
			sin(5 * u) + y[11] * cos(5 * u) + y[12] * sin(6 * u) + y[13] * cos(6 * u) + y[14] * sin(
			7 * u) + y[15] * cos(7 * u) + y[16] * sin(8 * u) + y[17] * cos(8 * u) + y[18] * sin(9 * 
			u) + y[19] * cos(9 * u) + y[20] * sin(10 * u), lty = i)
	}
}