for j in 1:10
δ= y[j,:]
diff_prob = ODEProblem(Proteo_Deg, u0,t_span,δ)
diff_sol = solve(diff_prob)
end
for
循环创建了一个局部作用域,因此在循环中引入的新变量只存在于该循环中。
diff_sol
。也许你想把它们保存在一个数组中?diff_solutions = Array{ODESolution}(undef, 10)
for j in 1:10
δ = y[j,:]
diff_prob = ODEProblem(Proteo_Deg, u0,t_span,δ)
diff_solutions[j] = solve(diff_prob)
end
现在,diff_solutions
包含了所有10个ODE解,并且也可以在循环外访问。如果这是在一个函数中,那么您可以在循环之后return diff_solutions
。