#OK to post
#Rebecca Embar, 2-23-2022, Homework 10

read `C10.txt`:

#1
CountInstabilities:=proc(M,W,pi) local count,n,i1,i2:
	count := 0:
	n := nops(pi):
	for i1 from 1 to n do
		for i2 from 1 to n do
			if not IsStable1(M,W,pi,i1,i2) then
				count := count + 1:
			fi:
		od:
	od:
	count:
end:

#2
GFstab:=proc(n,K,x) local p,co,i,rand:
	p := 0:
	for co from 1 to K do
		rand := RT(n):
		i := CountInstabilities(rand[1],rand[2],randperm(n)):
		p := p + x^i:
	od:
	p:
end: