I have an equation to calculate z from multiplying several rasters and numbers as seen below:

z <- a*b*c*d*e

The raster layers and values are as below:

set.seed(234)
r <- rast(nrows=10, ncols=10)

a <- setValues(r,runif(n = 100, min = 4, max=16))
b <- setValues(r,runif(n = 100, min = 0, max=100))
c <- setValues(r,runif(n = 100, min = 68, max=270))
d <- setValues(r,runif(n = 100, min = 0, max=1))
e <- setValues(r,runif(n = 100, min = 3, max=9))

I need to calculate uncertainty based on varying c and e by generating random rasters from within the minimum and maximum values, and by changing the raster values independently. In the end, I need to plot an uncertainty map.