Hi, Arden,
jnkkk
This is what I have told you last Thursday:  one works, one does not.
When I used 100,
> set.seed(100)
> s32(2,8,1,2,10)
Error in .Fortran("dqrls",: subroutine dqrls: 1 Inf value(s) in argument 5
Dumped

When I used 200,
> set.seed(200)
> s32(2,8,1,2,10)
$"Number of A 1st added in lm":
[1] 4

$"Number of A 1st added in glm":
[1] 3

$"Number of B 1st added in lm":
[1] 3

$"Number of B 1st added in glm":
[1] 2

$"Number of C 1st added in lm":
[1] 3

$"Number of C 1st added in glm":
[1] 5

$"Factor necessary in lm":
[1] 1

$"Factor necessary in glm":
[1] 2

$"1st A necessary lm":
[1] 1

$"1st A necessary glm":
[1] 1

> s32
function(r, n, s1, s2, k, ...)
{
	pa.lm <- numeric(k)
	pb.lm <- numeric(k)
	pc.lm <- numeric(k)
	pa.glm <- numeric(k)
	pb.glm <- numeric(k)
	pc.glm <- numeric(k)
	fact.lm <- numeric(k)
	fact.glm <- numeric(k)
	model.lm <- rep(NA, k)
	model.glm <- rep(NA, k)
	for(i in 1:k) {
		sim <- M3F2L(r, n, s1, s2)
		pa.lm[i] <- round(sim$a, 5)
		pb.lm[i] <- round(sim$b, 5)
		pc.lm[i] <- round(sim$c, 5)
		if(min(sim$a, sim$b, sim$c) < 0.05)
			model.lm[i] <- 1
		else model.lm[i] <- 0
		if(sim$a < sim$b & sim$a < sim$c)
			fact.lm[i] <- "A"
		else if(sim$b < sim$a & sim$b < sim$c)
			fact.lm[i] <- "B"
		else if(sim$c < sim$a & sim$c < sim$b)
			fact.lm[i] <- "C"
		pa.glm[i] <- round(sim$A, 5)
		pb.glm[i] <- round(sim$B, 5)
		pc.glm[i] <- round(sim$C, 5)
		if(min(sim$A, sim$B, sim$C) < 0.05)
			model.glm[i] <- 1
		else model.glm[i] <- 0
		if(sim$A < sim$B & sim$A < sim$C)
			fact.glm[i] <- "A"
		else if(sim$B < sim$A & sim$B < sim$C)
			fact.glm[i] <- "B"
		else if(sim$C < sim$A & sim$C < sim$B)
			fact.glm[i] <- "C"
		ans.lm <- data.frame(cbind(pa.lm, pb.lm, pc.lm, fact.lm, 
			model.lm))
		names(ans.lm) <- c("A", "B", "C", "Factor", "Model")
		ans.glm <- data.frame(cbind(pa.glm, pb.glm, pc.glm, fact.glm, 
			model.glm))
		names(ans.glm) <- c("A", "B", "C", "Factor", "Model")
	}
	list("Number of A 1st added in lm" = sum(ans.lm[, 4] == "A"), 
		"Number of A 1st added in glm" = sum(ans.glm[, 4] == "A"), 
		"Number of B 1st added in lm" = sum(ans.lm[, 4] == "B"), 
		"Number of B 1st added in glm" = sum(ans.glm[, 4] == "B"), 
		"Number of C 1st added in lm" = sum(ans.lm[, 4] == "C"), 
		"Number of C 1st added in glm" = sum(ans.glm[, 4] == "C"), 
		"Factor necessary in lm" = sum(ans.lm[, 5] == 1), 
		"Factor necessary in glm" = sum(ans.glm[, 5] == 1), 
		"1st A necessary lm" = sum(ans.lm[, 4] == "A" & ans.lm[, 5] == 
		1), "1st A necessary glm" = sum(ans.glm[, 4] == "A" & ans.glm[, 
		5] == 1))
}

> q()

Process S+3 finished at Mon Apr 26 16:09:10 1999