Проблема со временем выполнения из-за двойного цикла for

Question

Я создал массив для некоторых дифференциальных уравнений (deSolve) и хочу поместить в массив значения для нескольких предметов (n = 250) с переменной (n = 1000), но теперь я использую два для циклов и это занимает много времени выполнения. Я не знаю, как улучшить мой код. Есть предложения?

array.fruit <- array(NA, dim=c(X+1, length(yini)+1, 250, 1000), 
dimnames = list(NULL, c("time","Glycogen", "FatMass", "LeanMass", "ECF", "BW", "ATh", "BMI"))) #here I created an empty array
for (o in 1:250){
    for (i in 1:1000){  
        array.fruit[,,o,i] <- ode(func=model.fruit, y=yini, parms=pars, times=times1)
    }
}