Tuesday, May 16, 2023

Write your package

 #========================================================

#-- 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