x=30;h=10;p=0
for(p in 0:99){
n=h+h*p
f=function(x,n,a,mean,sd){
i=1;j=0;k=0;
while(i<=n){
  w=shapiro.test(rnorm(x,mean,sd))
  u=w$p.value
  if(u>a) j=j+1
   else if(u<a) k=k+1
  i=i+1
}
t=j/n
print(n);print(t);
plot(print(n),print(t))
}
f(x,n,0.05,0,1)
p=p+1
}
做出的是100张散点图，怎么合

用par（）？
http://www.statmethods.net/advgraphs/layout.html
R makes it easy to combine multiple plots into one overall graph, using either the
par( ) or layout( ) function.
With the par( ) function, you can include the option mfrow=c(nrows, ncols) to create a matrix of nrows x ncols plots that are filled in by row. mfcol=c(nrows, ncols) fills in the matrix by columns.
# 4 figures arranged in 2 rows and 2 columns
attach(mtcars)
par(mfrow=c(2,2))
plot(wt,mpg, main="Scatterplot of wt vs. mpg")
plot(wt,disp, main="Scatterplot of wt vs disp")
hist(wt, main="Histogram of wt")
boxplot(wt, main="Boxplot of wt")

100张散点图合在一张太大了，会导致每一个小图太小，不太好看，至少在图形界面是这样的。除非你用PDF文件保存起来，设置图形文件大小。

想过把每次做的结果组成一个矩阵然后提取矩阵中的变量做出图，但是结果没有出来矩阵的形式

x=30;h=10;p=0
for(p in 0:5){
n=h+h*p
f=function(x,n,a,mean,sd){
i=1;j=0;k=0;
while(i<=n){
  w=shapiro.test(rnorm(x,mean,sd))
  u=w$p.value
  if(u>a) j=j+1
   else if(u<a) k=k+1
  i=i+1
}
t=j/n
matrix(c(print(n),print(t)),1,2)
dev.new()
plot(n,t)
}
f(x,n,0.05,0,1)
p=p+1
}