In R, if I have a function
myfun<-function(ys, T, N, beta, gamma, sigma) {...}
where beta and gamma are scalars.
If I have fixe
Look at the mapply and Vectorize functions along with expand.grid for the all combinations part.