#
# Die Mohr-Coulom'sche Grenzgerade
#
Coh <- 50  # kN/m^2
beta <- 25 # Grad
sigma <- seq(0,200,10)
tau <- Coh + tan(pi*beta/180)*sigma
plot(sigma,tau,type="o",xlim=c(0,1.5*max(sigma1,sigma3)),
       ylim=c(0,1.5*max(sigma1,sigma3)),
       main="Mohr-Coulomb'sche\nGrenzgerade",
       xlab=expression(paste("Normalspannung  ",sigma)),
       ylab=expression(paste("Scherspannung   ",tau)))