rm(list=ls())
#--------------------------------------------------------------------------
#   Berechnung von spezifischen Oberflächen
#--------------------------------------------------------------------------
#
# Funktion zur Berechnung der spezifischen Oberflächen:
spec_surface <- function(rho_s,d) 
{
  spec_surface <- 3/(rho_s*d/2)
  # Angaben in cm
}
#
# Datenvektor mit Korndurchmessern
d <- c(0.2,0.063,0.02,0.0063,0.002)          # mm
d_cm <- d*0.1                                # cm
#
# Datenvektor mit Dichten:
rho_s <- c(1.0,1.5,2.0,2.5,2.65)          # g/cm^3
#
# Erzeugen der Datenmatrix mit den Funktionen "outer" und "spec_surface"
sp_srf <- outer(rho_s,d_cm,FUN="spec_surface")
# 
A <- data.frame(round(sp_srf[,1:5],1),round(rho_s,2)) 
names(A) <- c(paste(d, c(" mm"), sep=""),"rho_s")
cat("--------------------------------------------------\n")
cat(" Spezifische Oberflächen (cm^2/g):\n")
cat("--------------------------------------------------\n")
print(A)
cat("--------------------------------------------------\n")
#write.table(A,"Oberflächen.txt",sep="\t")
#
#   --- Schluss ---
#