# 位置座標(x,y)に関する(qx,qy)次の多項式基底ベクトルのRスクリプトを自動生成

# 例えばqx=2,qy=2なら,次のRスクリプトを生成

# c(1,x^1,x^2,y^1,y^2,x^1*y^1,x^1*y^2,x^2*y^1,x^2*y^2)

make.basis <- function(x,y,qx,qy){

tmp <- "c(1"

 for(i in 1:qx)tmp<-paste(tmp,",x^",i,sep="")

 for(i in 1:qy)tmp<-paste(tmp,",y^",i,sep="")

 for(i in 1:qx){

   for(j in 1:qy){

     tmp<-paste(tmp,",x^",i,"*y^",j,sep="")

   }

 }

 tmp<-paste(tmp,")",sep="")

 return(eval(parse(text=tmp)))

}