Dear list, Is there a way to do something like the following pseudo-code - without for loop? complexFn <- function(a,b){ ... return(c) } x[i] = complexFn(x[i-1], x[i-2]) I have found "filter" but would like to do something more complicated than multiplication with x[i-1] and x[i-2]. Thanks in advance! gabor