T<-seq(0,35,0.5)
es<-6.11*10^((7.48*T/(237+T)))
plot(T,es,xlab="Temperatur in Celsius", ylab="Dampfdruck in hPa",
  xlim=c(0, 35),ylim=c(0, 35), pch=1, lty=1, col="blue", axes=FALSE)

axis(1, seq(0,35,5))
axis(2, seq(0,35,5))
grid (NULL,NULL, lty = 3, col = "grey")
curve(6.11*10^((7.48*x/(237+x))), add = TRUE)