使用montecarlo从乘法获得分布

Question

来自a = b * c 我需要找到均值，“ a”的P10，P50和P90，从“ b”和“ c”的分布中提取，我试图为此使用Montecarlo，但是我得到的是一个错误。 我知道这有点愚蠢，但是您能帮我吗 请在下面找到我使用的R代码

library(MonteCarlo)
mult<-function(a,b) {sample<-rnorm(b,c)
a<-sample*sample
mean<-mean(a)
return(list("mean"=mean))}
b_grid<-c(50,100,250,500)
c_grid<-seq(1,2,0.2)

param_list=list("b"=b_grid, "c"=c_grid)

MC_result<-MonteCarlo(func=mult, nrep=1000, param_list=param_list)