R语言处理随机变量函数

想问一下有没有R语言的包可以处理下面的问题：
有两个随机变量x, y, 已知他们的概率分布，想计算(x + y)的概率分布

只要数值解的话
integrate做积分。
联合分布函数
fz=function(z) {integrate(Vectorize(function(x) (fx(x)fy(z-x))),-Inf,Inf)}
未测试，可以是这个思路。