func_vol <- function(dap, sobr, area_talhao){
altura <- exp(3.7004 - 7.6165/dap)
vol <- 0.0000599458 * dap^1.82083 * altura^0.97645
vol_ha <- vol * 1111 * sobr
vol_talhao <- vol_ha * area_talhao
#cat("Volume no talhão:", vol_talhao, "m³\n")
return(vol_talhao)
}
vol_talhao1 <- func_vol(22, 90/100, 40)
vol_talhao2 <- func_vol(19, 92/100, 35)
vol_talhao3 <- func_vol(25, 91/100, 42)
vol_total <- vol_talhao1 + vol_talhao2 + vol_talhao3
vol_total
dap <- 22
sobr <- 90/100
area_talhao <- 40
altura <- exp(3.7004 - 7.6165/dap)
vol <- 0.0000599458 * dap^1.82083 * altura^0.97645
vol_ha <- vol * 1111 * sobr
vol_talhao <- vol_ha * area_talhao
cat("Volume no talhão:", vol_talhao, "m³")
dap <- 19
sobr <- 92/100
area_talhao <- 35
altura <- exp(3.7004 - 7.6165/dap)
vol <- 0.0000599458 * dap^1.82083 * altura^0.97645
vol_ha <- vol * 1111 * sobr
vol_talhao <- vol_ha * area_talhao
cat("Volume no talhão:", vol_talhao, "m³")
dap <- 25
sobr <- 91/100
area_talhao <- 42
altura <- exp(3.7004 - 7.6165/dap)
vol <- 0.0000599458 * dap^1.82083 * altura^0.97645
vol_ha <- vol * 1111 * sobr
vol_talhao <- vol_ha * area_talhao
cat("Volume no talhão:", vol_talhao, "m³")