如题：

xiao<- function (co, intercept, power, value){
        len.co<-length(co);
        len.int<-length(intercept);
        len.pow<-length(power);
        len.val<-length(value);
       
        n<- len.co*len.int*len.pow*len.val;
        equ.solut<-rep(0, n);
       
        for(i in 1:len.co){
                for (j in 1:len.int){
                        for(k in 1:len.pow) {
                                for(l in 1:len.val){
                                        index<-i*j*k*l;
                                        equ.solut[index]<-(value[l]/co[i]-intercept[j])^(1/power[k])-1
                                        par[index,]<-list(as.numeric(c(equ.solut[index], co[i],intercept[j],power[k],value[l])))
                                }
                        }
                }
        }
       
       
        return(par)
}

co<-c(100)
intercept<-c(0,1)
power<-c(10)
value=c(1000, 800, 600)

xiao(co, intercept, power, value)