用于模拟反射布朗运动的R代码的问题

Question

我使用的是反射布朗运动的以下定义：

，

有时也称为受控布朗运动。最重要的是，反射的布朗运动应始终位于零以上。

根据上面对反射的定义，我具有以下代码用于通用函数的反射：

reflected <- function(f) {
    g = function(t) {
        if (t==0) return(0)
        else {
            x = optimize(f, lower=0, upper=t, maximum=F)
            return(f(t)-x$objective)
        }
    }
    return(Vectorize(g))
}

但是以某种方式这不能产生正确的过程。我的代码哪里错了？ Brownian motion and reflected Brownian motion (black and blue resp.)