#========================================================
#-- Multiple functionsctions
#=========================================================
DPGV <- function(x,y,z) {
Bulk_density <- (y/z)
Porosity <- (1-(Bulk_density/2.65))
Gravimetric_water <- ((x-y)/y)
Volumatric_water<- Gravimetric_water*Bulk_density
volumtric_air<-Porosity-Volumatric_water
Void_ratio<- Porosity/(1-Porosity)
out <- table(Bulk_density,Porosity,Gravimetric_water,Volumatric_water,volumtric_air,Void_ratio)
return(out)
}
#----------------------------------------------------------------------------
aggregate_diameter <- function(a,b,x) {
MWD <- ((a+b)/2)*(x/200)
GMD <- ((x/200)*(log(a+b))/(x/200))
out <- table(MWD,GMD)
return(out)
}
#=============================================================
#== Single function
#=============================================================
soil_weight <- function(x,y) {
result <- 100 * 100 * x*y
print(result)
}
##############################################################
No comments:
Post a Comment