# NOT RUN {
def <- c(0, 0.0230, 0.0352, 0.0605, 0.1070, 0.1750, 0.2525, 0.3395, 0.4250)
pres <- c(1, 12.5, 25, 50, 100, 200, 400, 800, 1600)
VR <- voidratio(wetsoil = 170.62, drysoil = 134.08, diam.cylinder = 6.95,
height.cylinder = 2.5, dens.particle = 2.61, def)
VR
plot(VR ~ pres, type = "b",
ylab = "Void ratio",
xlab = "Applied stress (kPa)",
main = "Compression curve",
log = "x")
# End (not run)
# }
Run the code above in your browser using DataLab