求救怎样才能通过模拟产生pvalue值～～～～

我现在想通过编程模拟通过已经得到的chi-square值已经各个方格的期望值来求得卡方值相应的p-value值（因为我要求的数据量比较多，如果用spss软件，我累死也输不完呀。所以编程～～～）
举一个例子如下：以下是通过spss处理得到的


或直接打开～～～～http://hiphotos.baidu.com/wubingran/pic/item/c3800113ef72d6caf7039eb9.jpg

上面第一个表格分为gr1 和gr2两组，其中每个格子俩面有两个数，上面的那个数是观测数，下面的那个是理论数，如29是观测数，29.3是理论数。

第二个表格就是卡方检验了
其中，9.372就是卡方值，而0.404就是相应的p-value值，也就是平常说的概率值，我现在就是想知道怎么通过随机种子模拟来求出p-value值

谁来帮忙出个点子，编程一下～～求求这个p值～～～～～